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CHAPITRE XXIV. 

GENERALITIES SUR LA THEORIE DE LA LUNE. 


311. La iheorie de la Lune doit 6 tre exposde en deux parties; 
dans la premiere partie, on cherche quel serait le mouvement de la 
Lune, si la Lune, le Soleil et la Terre existaient seuls el dtaien, 
reduits a des points materiels; dans la seconde partie, on cherche 
comment ce mouvement est trouble par Fattraction des planetes et 
par l’influence de Faplatissement terrestre. 

La premiere partie n’est done qu’un cas particulier du probleme 
des trois corps, et la difficulty ne provient que de la grandeur rela- 
tivement considerable des perturbations produites. Le rapport de 
la force perturbatrice a Fattraction du corps central est, comme 
nous Favons vu au Chapitre II (p. 537 ), de Fordre de 

m /AC\* 
m, \B0J 9 

m k etant la masse du corps troublant, m 7 celle du corps central 
AC et BC les distances mutuelles des trois corps. Ce rapport est 
le produit de deux facteurs dont Fun est le rapport des masses et 
Fautre le cube du rapport des distances, Dans le cas des planetes, 
le premier facteur est tr£s petit, et le second fini ; dans le cas de la 
Lune, au contraire, le premier facteur est grand et le second tres 
, petit. II en rdsulte que le produit des deux facteurs est petit sans 
-n(a). 


x 
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doute, sans quoi le probleme ne pourrait se resoudre par approxi- 
mations successives, mais beaucoup moins petit que dans le eas 
des planetes, de sorte que Fapproximation est beaucoup plus 
lente. 

La difficult^ et l’importance du probleme ont amene un grand 
nombre de geometres a s’en occuper et l 5 etude de leurs recherch.es 
est extremement interessante au point de vue historique; on en 
lira avec profit l’expose dans le Tome Hide la iMecctnique celeste , 
de Tisserand. Mais aujourd’hui on peut dire qu’il n’y a plus que 
trois methodes qui comptent : celle de Hansen, celle de Delaunay 
et celle de Hill-Brown. 

Celle de Hansen est celle qui a servi a construire les Tables 
actuellement en usage. Ces Tables sont d’une exactitude remar- 
quable et, si elles s’ecartent des observations, les divergences ne 
sont pas dues a un defaut de la methode (au moins tant qu on ne 
considere que le Soleil, la Terre et la Lune), puisque les autres 
methodes, plus satisfaisantes au point de vue theorique, semblent 
devoir conduire aux memes divergences, mais a l’omission de 
quelque terme provenant de Faction des planetes ou a quelque 
cause inconnue. Cette methode est toutefois tres compliqu^e et je 
renonce a l’exposer ici, renvoyant soit aux Ouvrages originaux de 
Hansen, soit au resume qu’on trouvera au Chapitre XVII du 
Tome III de Tisserand. Le succes de Hansen parait du surtout a 
son habilete et a sa patience personnelles, et aussi a ce fait qu’il a 
cherchd directement les valeurs numeriques des coefficients sans 
passer par une expression algdbrique ou les constantes seraient 
representees par des lettres. 

Delaunay a fait tout le contraire, tous ses coefficients sont 
exprimes par des series ou figurent les differentes constantes du 
mouvement de la Lune et dont les coefficients sont des nombres 
rationneis exactement determines. Ces formules sont done appli- 
cables, non seulement a la Lune, mais a un satellite quelconque 
(en le supposant unique). 11 suffirait d’y substituer, au lieu des 
constantes relatives a la Lune, celles qui se rapporteraient a ce 
satellite. II serait ais 6 egalement de voir imm^diatement quelle 
serai t Finfluence d’une correction apportee & Fun des elements de 
la Lune. La determination des nombres rationneis qui servent de 
coefficients a exige un travail dnorme; si M. Andoyer a decouvert 
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quelques ei'reurs, c’est seulement dans les termes d’ordre tres 
eleve, et qui n’entrent pas en ligne de compte avec l’approxima- 
tion habituelle des Tables. Ce travail algebrique est reste long- 
temps inutilise et c’est seulement tout dernierement qu’il a ete 
reduit en nombres. 

Brown a pris une position intermediaire. Ses coefficients ne 
sont ni purement numdriques comme ceux de Hansen, ni pure- 
ment analytiques comme ceux de Delaunay. Us se presentent sous 
forme de series procedant suivant les puissances des divers ele- 
ments, le rapport des moyens mouvements excepte; les coeffi- 
cients de ces series sont calcules numeriquement, mais ces coeffi- 
cients ne sont plus des nombres rationnels, ce sont des fonctions 
du rapport des moyens mouvements, qu’on pourrait egalement 
developper en series, mais dont on se borne a determiner la valeur 
numerique. Comme, d’autre part, la methode de Brown £tait 
beaucoup plus directe que les autres, il a pu pousser l’approxima- 
tion beaucoup plus loin que ses devanciers. 

La methode que nous exposerons ici est celle de Brown avec 
quelques modifications; c’est elle, en effet, qui nous permet le 
mieux d’utiliser les resultats obtenus dans le Tome 1 et de ratta- 
cher ainsi la theorie de la Lune a la thdorie generate du probleme 
des trois corps. 

312. Nous adopterons les notations du Chapitre II, Tome I, et 
nous renverrons en particular au n° 42. Nous d^signerons done 
par 

iU]_ — m2 = m 3 la masse de la Lune, 
m^= — m Q la masse du Soleil, 

in- — m & = m 9 la masse de la Terre; 

par 

a?i, x*, #3 las coordonnees de la Lune, 

x$, xe les coordonnees du Soleil, 
x 7 , Xsj x 9 les coordonnees de la Terre ; 

par A, B, C les positions des trois corps : Lune, Soleil, Terre, 
par D le centre de gravite du systeme Lune, Terre; par 

x' u x t , #3 les trois projections du vecteur AC, 
a?*, » BD. 
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Nous poserons 


tlX a —— lYXet — ■ TiX « — ■ 


m i m 7 


m< — m h = ni R = 


Wr 


W i ( m ! -h />? 7 "1 

nt i -f- WI4. -+- 


/ / 


(zl 

+ 4 

+ 4 ). 

t.= i (4 

+ 4 

- 5 ) 

\ m\ 

m 2 

mj 

- a \ 

% 

^6 / 


de sorte que represente la force vive de la Lune dans son mou- 
vement relatif par rapport a la Terre, en lui attribuant fictivement 
la masse m\, tandis que T 2 represente la force vive du Soleil dans 
son mouvement relatif par rapport au point D, en lui attribuant 
fictivement la masse m', r 

Nous poserons en outre (voi? % t. I, p. 55) 

T . nii fT -f- m 7 ) 

lJl = AC ~ 3 “ BD 

U 3 = m-iTrii (g-g - Xb ) m * m ~ ('BD ~ BO,)’ 

F = Ti -t" T o — f- U 1 -+■ U •> -+■ (J 3 — '+■ 1 ), 

= To -+*• U 2 , in j <^1 = Tj -f- U\ -4- U3. 


Nous avons vu que U 3 est comparable au de et de U] et 
au o ii o W( Fo ^2 et de U 2 , ce qui nous permet de negliger U 3 
devant <$ 0 . 


313. Les equations du mouvement prennent la forme canonique 

dxi __ dy\ ___ d¥ 

dt dy\ ’ dt dx\ 

Si nous faisons i — 4, o ? 6, nous voyons que les d^rivees de 


( l ) Je trouve plus avantageux de poser 
y'=:m\y 

au lieu de 

F = a> 1 >^-m 1 <£ n 
y'-m x y\ 


comme au Tome I. 
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T, + U ( sont nulles, que celles de U 3 sont nEgligeables devant 
celles de <£ 0 , de sorte qu’il reste 


(2) 


dx\ _ d&p dy[ 

dt ciy'i dt 


d<$> o 
dx'i 


( i = 4 , 5 , 6 ), 


ce qui prouve que le mouvement du Soleil par rapport au point D 
peut Etre regarde comme keplerien. 

Si nous faisons i= i, 2, 3, ies dErivEes de <£ 0 sont nulles et il 
reste 


(3) 


dx\ 

dt 


Z l ! 3 

dy L 


dyt 

dt 


— m 


, d& i 
1 dx[ 


(* -i, 2, 3). 


Les seconds membres des equations (3) dependent encore 
de x\, xi, x\ \ mais ces quantity etant determinees par les Equa- 
tions (2) peuvent Etre regardees comme des fonctions connues du 
temps. Si on les remplace alors par leurs valeurs en fonctions du 
temps, les Equations (3) se prEsentent sous la forme canonique, 
mais de telle facon que la fonction caractEristique 4 dEpende 
explicitement du temps (comme au n° 12). 

D’autre part, si Ton veut Eviter la prEsence du petit facteur m 
il suffit de poser comme au n’ 121 (t. I, p. 162) : 

Yi — m\ y\ (* = i, 2 , 3). 


Les Equations (3) restent canoniques et deviennent 


(3 bis ) 


dx'i _ d$ 1 dy"i __ d$ x 

dt dy’i 1 dt dx'i 


31-4. 11 faut maintenant appliquer les principes des Chapitres X 
et XII. Les rEsultats en ont Ete rEsumes en particulier au n° 177. 
On y voit que les ElEments osculateurs peuvent Etre developpes 
suivant les puissances de p. et des expressions 

(4) E*cos«pi, E/ c sinwPA- ? 

suivant les cosinus et les sinus des multiples des arguments w", les 
coefficients des dEveloppements dEpendant encore de rt constantes 
d’intEgration W / (s’il y a n 1 coi'ps). 

Les E sont des constantes d’integration de l’ordre des excentri- 
citEs et des inclinaisons. Les arguments w l et w u varient propor- 
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tionnellement au temps, et I’on a 

= — Y/* ( -+■ w 'i = ; 

les zn et les to 7 sont des constantes d’integration ; les n' et les y' sont 
des constantes developpables (cf. n° 179) suivant les puissances 
de p. et des E-. 

Nous avons vu ensuite au n u J92 que les coordonnees heliocen- 
triques (ici geocentriques), c'est-a-dire x \ , x[ 2 , x' z , sont develop- 
pables de la merae maniere, et qn’il en est encore de meme de y\ , 
y’z * ^ es developpements prennent d’ailleurs la forme (cf. t. 1, 
p. 3a 7 ) 

(5) 2 A E?) dn (2 kW " ^2 piV ) ’ 
les q , les k et les p etant des entiers et Ton a 

( 6 ) 2 />==0 

dans les developpements de x' 3: x ' 6 , r 3 , , 

( 7 ,; ^ = 1 

dans ceux de 3?^ 5 T n ^ 110) 1G.2, 1ST 

et 192, p. 328). 

Nous n’avons que des cosinus dans les developpements de 

x\, r' s , r's- 

Nous n’avons au contraire que des sinus dans ceux de 
^57 y i ? rs, y 45 re 
(cf. n° 190 et aussi n" 192). 

Le nombre des arguments w' est de 4? celui des arguments w u 
de 2 (n° 193). Mais l’un des moyens mouvements y 4 est nul et w\ 
se r^duit a une constante. Si d'ailleurs on prend le plan invariable 
pour plan des x K x 2 , on a E 4 = o, v de sorte que l’argument w\ ne 
figure plus dans les developpements et qu’il ne reste plus que 
cinq arguments : 


(8) 



GENERALITIES SUR LA THEORIE DE LA LUXE. 


7 


Une nouvelle simplification provient de ce fait que la masse de 
la Lane etant tres petite, le mouvement du Soleil par rapport au 
point D peut £tre regards comme keplerien. Le plan invariable 
que nous avons pris pour plan des x K x 2 n’est alors autre chose 
que le plan de Porbite solaire. D’autre part qui n'est autre 
chose, au signe pres, que la longitude du perihelie de Porbite 
solaire keplerienne, se reduit a une constante, et il ne nous reste 
plus que quatre arguments distincts : w \ , w ' 2 , w", w", dont il est 
aise d’apercevoir la signification : w\ est la longitude moyenne de 
la Lune, je ne veux pas dire la longitude moyenne sur Porbite 
osculatrice, mais pour ainsi dire la longitude moyenne moyenne; 

c’est la longitude moyenne du Soleil; — c’est la longitude 
moyenne du perigee lunaire; — w/ 2 , c’est la longitude moyenne du 
noeud. 

De m£me E, est une constante qui joue un role analogue a 
Pexcentricit6 lunaire; E 2 joue le role de Pinclinaison; E 3 joue le 
role de Pexcentricite solaire, et nous pouvons m£me profiter de 
Pindetermination de cette constante E 3 pour supposer qu’elle est 
pr^cisement £gale a cette excentricite solaire. 

Dans les developpements de 

x 'i, y\i y 2 > 

Pexposant de E 2 et le coefficient de w[ 2 seront to uj ours pairs. Ils 
seront toujours impairs, au contraire, dans ceux de 

x y 3 

( cf . n° 191). 

Si E 3 = o, c 7 est-a-dire si Porbite solaire est supposee circulaire, 
nos developpements ne dependent plus de w 3 , mais de l’argument 

k\ w\ •+■ k* w\ -f- p\ -+- P2W21 

oh l’expression 

k 1 — {— /to — p 1 — p 0 

est £gale a z6ro pour les distances mutuelles et pour # 3 , et a 1 
pour x\ et x 2 . 

Si de plus Pinclinaison est nulle, c’est-a-dire si E 2 =o 7 les 
termes dependant de w[ 2 disparaissent; on a done x' :i = o, et dans 
les distances mutuelles des trois corps figure seulement Pargument 

kt w\ -h kt w\ •+• pi w\. 
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OU 

Pi—k\-\-k*. 

Elies dependent alors seulement des deux arguments 

*>'[-+■ w\, 

et elles sont developpables suivant les puissances de 
E! cos(hOo h- w\ ), E! sin ( w\ -+- w\ ). 

C’est le cas du probleme restreint. 

Si nous revenons a l’hypothese E 3 =o, Tintegrale de 

Jacobi a encore lieu. 

315. Nous avons pose plus haut (n° 312) : 

/TZi <f>! = Tj-t-Ui-f- U 3 ; 


nous ferons d’abord observer : 

i° Que U 3 est beaucoup plus petit que T, -h U* ; 

2 ° Que T 4 et U* dependent seulement des masses et m 1 et 
des coordonnees de la Lune, c’est-a-dire des inconnues x \ , 

y 1 1 y 2) y % > 

3° Que U 3 depend en outre de la masse m h du Soleil et des 
coordonnees du Soleil, qui sont des fonctions connues du temps 
et des elements de l’orbite solaire. 


D’autre part, comme AC est beaucoup plus petit que BD, on a 
(cf. t. I, p. 55) : 


( 9 ) 



m\ /rtgzb m-tm'l 

{mi -+- 


P * 


AG" 
BD"-*- 1 ? 


P„ etant une fonction de Tangle des deux directions AC et BD, 
fonction definie Tome I, page 49? 1& serie (p) convergeant 

tres rapidement. 

Nous voyons comment CJ 3 depend des masses etdes elements du 
Soleil. 


i° II est proportionnel a m h ; 

° ^ ne depend < I ue du rapport 


2 { 
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3 ° U 3 depend en outre des elements de Forbite terrestre, a 
sa\ oir : de la longitude du perihelie — qui est une constante ; de 
la longitude moyenne du Soleil — w\, qui varie proportionnelle- 
ment au temps; de l’excentricite solaire, que j’ai designee plus 
haut par E 3 , et enfin du demi grand axe de Forbite solaire que je 
designerai par a ! . 


Voyons comment chacun des termes de la serie (9) depend 
de a'. La distance AC ne depend que des coordonnees de la Lune, 
le facteur P /2 depend d’un angle; il ne dependra done pas de 
mais seulement des autres elements du Soleil. Quant a BD ,i+l , 
e’est la puissance 72 + 1*®"’“ d’une longueur, il sera done propor- 

tionnel a a n+i . Le rapport sera independant de a ! . Nous pou- 
vons done ecrirc 


(ioj 


U 3 



nt\ m'J ± m-m'f 
( m j -t" m- ) !l 


P„AO 


/ a! yn-i 1 

[m) 


et comme tous lesfacteurs sauf le dernier sont independants de a ! , 
nous aurons le developpement de U 3 suivant les puissances 
decroissantes de a! . Le premier terme du developpement (10) est 


— m>+ 


m 1 m- t 
4 - in 7 



a' 6 


? 


de sox'te que nous pouvons ecrire 



Q„ etant independant de m 4 et de a ’ . 

Le facteur ~ est petit, e’est lui qui joue le rdle de p. D’ail- 

leurs a 1 est une constante et notre function perturbatrice se trouve 
developpee en une serie ires convergente procedant suivant les 

puissances de ~ Cette constante ~ pourrait done aussi jouer le 

rdle de ;jl; de sorte qu’en appliquant les principes precedents, 
nous trouverions que nos inconnues peuvent se developper suivant 
les puissances de 


W4 

a 7 * 


et de 


r 
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D’autre part, la troisieme loi de Kepler nous donne 


d’ou 


( m--r- m 4 ') = ci' 3 nl, 



n\ 


i 


111 ; 
ni'+ 


Le second facteur est Line constante connue; elle est si voisine 
de Funite que nous pouvons prendre simplement 

;jl = n\. 

II semble done que nous devions conclure que nos coordonnees 
vont §tre developpables suivant les puissances de n\ et de 

Mais avant d’adopter cette conclusion, il eonvient d’y regarder de 
plus pres. <t> t et U 2 dependent de n 2 de deux manieres : 


i° Directement, a cause du facteur p. = qui affecte U 3 ; 

2 ° Indirectement, parce que U 3 depend en outre de w" a , qui est 
£gal a n 2 1 . 


Pour pouvoir appliquer les principes du Chapitre X, il faut 
regarder n 2 (en tant qu’il entre par w”. 2 ) et p. comme deux variables 
independantes. Dans ces conditions, nos developpements proce- 
deront suivant les puissances de pi, mais les coefficients seront des 
fonctions des difFerentes constantes et en particulier de n 2 . Or 
nous avons vuque les integrations introduisent des petits diviseurs 
de la forme k\ n { -h k 2 rt 2 \ Fun de ces petits diviseurs est pr£cis£- 
ment n 2 . 

Considerons done un terme quelconque du developpement 
dependant de ce petit diviseur qui est n 2 ; soit a Fexposant de p 
et j$ celui du petit diviseur, de telle facon que notre terme con- 
tienne en facteur 

L’expression a — | represente ce que nous avons appel£ la 

clctsse du terme, et nous avons d^montrd au n° 198 que cette 
classe dtait toujours positive ou nulle. 

L’exposant 2 a — J3 de n 2 ne peut done pas etre negatif, mais il 
n’y a pas de raison pour qu’il soit pair. 
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Ce n'est done pas suivant les puissances de n* et de ~ que 
nos developpements procedent , mats suivant celles de n>> et 
de Nous verrons dans la suite, a la fin du Chapitre XXVIII, 

qu’i] peut m&me, mais seulement pour des termes d’ordre tres 
elev6, s’introduire de petits diviseurs en n: 2 ou n\. II en resulte 
que n 2 pourra, dans certains termes du developpement, figurer a 
une puissance negative . 

316. Nous trouvons done finalement 

x' = f{m u m 7 , m; ; n u E 1? E 2 ; u>\, w\ \ a\ E 3 , w\, w' z ). 

Nos coordonn^es dependent en efiet des trois masses, des quatre 
arguments w \ , w'', w \ , ; des deux constantes d’integration E* 

et E 2 introduites plus haut; d’une troisieme constante pour 
laquelle nous pouvons choisir le moyen mouvement n K ; des ele- 
ments de Forbite solaire a\ E 3 et (v' 3 . 

Nous pouvons remplacer m 4 par a lz n\ et introduire une nou- 
velle constante a definie par 

7?2 1 -f- m-i — n\ a 3 . 


Alors m K et m 7 sont des fonctions de — et de ri]a z , de sorte 

7 771 - 1 7 

que nous pouvons ecrire 

X = f{jn-’ a ’ E2 ’ w '*• ( *'*’ a> ’ Es ’ “’»)• 

II faut faire intervenir maintenant des considerations d’homo- 
g^neite : 

a et a f sont des longueurs, de m£me que les x f ; 

— et — sont des temps: 

Jll 7l 2 L 7 

les E sont des nombres, ou du moins on peut profiler de Findeter- 
mination de leur definition pour le supposer; 

— et les w sont des nombres. 

772 7 

Dans ces conditions, pour que la formule soit independante 
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des unites de longueur et de temps, on doit avoir 

, Jn* a\ 
x = a ft — 9 -7 ) • 

J \ni cl J 

Je n’ecris pas explicitement celles des variables qui sont des 
nombres. Nous poserons 

n 2 , a 

777 = 

7l\ CL 

et comme nos expressions sont developpables suivant les puis- 
sances de n 2 et de nous voyons qu’elles le seront suivant les 

puissances de m f et de a; la constante m est le rapport des moyens 
mouvements; a est ce qu’on appelle la parallaxe . 

On aura d’abord, en differentiant par rapport au temps, 


n= 


m\ dx\ __ r dx'i 


m x dt 


1 -b 


mi dt 9 


m-j 


et, par consequent, par raison d’homogeneitd, 

y* — =ntaf(m', a), 

et, d’autre part, 

*1= »!«*/(»', «; § > Jg; e», «-',)■ 


317. II convient encore de faire quelques remarques au sujet 
de la symetrie. Nos ddveloppements procedent suivant les puis- 
sances de 


m , 


„ cos , 
Ei . w\ 
sin 1 


„ cos . 
E 2 . 
sin 2 


e 3 cos w ' 3 

sin 


etl’argument du terme general est 

h w\ -b k 2 w\ -b pi w\ -b p* -b pz vd % . 


On a — 2/> = opour #3 et== 1 P our ^i et x\ v D’ailleurs p 2 

est pair pour x\ et x[ 2 et impair pour x^ ; d’oti il suit que 

£i-b £ 2 — p { — pz 
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est toujours impair. D’autre part, p 3 est toujours de m^me parite 
qae l’exposant de E 3 ; done ko — p { — i est toujours de parite 

opposee a cet exposant. 

Soit maintenant q 0 l’exposant de la paraliaxe a; nous voyons 
que la position du Soleil ne change pas quand on change 

a\ 

en 

— a ' , W% T., W ! 3 TC. 

Dans ce cas, les coordonnees du Soleil ne changeant pas, rien 
ne doit changer puisque, dans nos equations, a f : ne s’in- 

troduisent que par les coordonnees du Soleil. 

Or, dans ces conditions, un tenne quelconque se trouve multi- 
plie par 

( iy/o + ^2-^^3; 

on doit done avoir 

q o -H k 2 -f - pz = o (mod. i). 

D’autre part nous avons trouve 

ki-^r — p z =B i 

pour les trois coordonnees x \ , x'. 2: x :V et par consequent 

*i -f- q Q z= i. 

Pour les distances mutuelles nous aurions trouve 


k t -h k« — p x — p%~ ps — o, p. 2 = o, 
/ti -+- k$ — pi — Pi^ o, 

d’oii 

■+* ^1 •+■ qo = °* 

Un peu plus loin (n <l 320) nous poserons 

x — x\ cos w 2 ~h x\ sin «p 2 , 
y — — x\ sin w s - 4 - x f 2 cos t*> 2 . 


On voit que si k 2 est pair dans le developpement de x\ elx\ 2 , 
il sera impair dans celui de x et y et inversement. On aura done, 
pour x eljK, 

^-t-A'o+p 3 = L 

et pour les termes ind^pendants de a et de E 3 , e’est-a-dire si 
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gr 0= r/> 3 = O, on aura 


k-2 == I . 


318. Nous allons donner une formule importante pour ce qui 
va suivre, et, pour l’etablir, nous reproduirons a peu pres le rai- 
sonnement du n° 120, en nous appuyantsur le theoreme du n°16. 
B’apres ce theoreme, si i’on a un systeme d’equations canoniques 

c lx __ dF dy __ dF 

dt dy 9 dt dx 9 

et si les x et les y sont supposes exprimes en fonctions des con- 
stantes d ? integration a et du temps t , on a l’identite 

dy = d& Aflfa-F dt, 


ou Q est defini par l’equation 


S-'*2-S 


dy 


-'-2 


X 


dF 
dx 9 


et ou les A sont independants du temps et dependent seulement 
des a. 

lei nous avons les equations (3 bis) 

dx'i _ d $> ! rfyj __ 

dt dy'i 7 dt dx'i 

Elies sont bien de la forme canonique, mais <£* depend non 
seulement des inconnues x 1 et y , mais encore du temps, ce qui 
nous oblige & appliquer l’artifice du n° J2. Si 3>* depend expliei- 
tement du temps, e’est par l’intermediaire des coordonnees du 
Soleil, qui elles-memes sont des fonctions periodiques de l’argu- 
ment wl x . Nous introduirons done deux variables auxiliaires u et p. 
Nous pouvons ecrire <$ f sous la forme 

*i ( *>%), 

et poser 

F'= u)-h n 2 v. 

Nous avons alors, si nous voulons que u = w\ 1 = n 2 t - |- 


du dW 

-dt~ n% =lb 9 
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et nous pouvons ecrire les equations canoniques 


/ don' 

dF’ 

df_ 

dF’ 

| Hi 

~ dy"’ 

dt 

dx ' 3 

< 

I du 

dF' 

dv 

_ dF' 

' Hi 

“ dv 9 

dt 

du 


Nous ajoutons arbitrairement la quatrieme equation qui peut £tre 
regardee comme la definition de Nous poserons d’ailleurs, pour 
unifier les notations, 

W\ = w\ , W 2 = W'J, 

(V 3 = w \ , = w 2 , 

et, d’autre part, 

Wi — mt ■+■ , 

ce qui ne change pas la definition de /i M to, etTOo. 

Les equations restent canoniques et F' ne depend plus explicite- 
ment du temps; nous pouvons done ecrire 

(12) 2 #' dy"-+- u dv = dQ. A da — F' dt. 

Mais 

F' dt = (<J>! n 2 v) dt — dt-\- v(du — dm 2 ), 

ce qui nous permet d’^crire 

(1 3 ) dy n = d(Q — da — dt -t- v dxs^ 

Quelles sont nos constantes d’integration a? Nous avons d’abord m ! , 
E, et E 2 , quenous appellerons les constantes j 3 ; nous avons ensuite 
les constantes to qui figurent dans les arguments w. 

Le syst£me (1 1) ^tant du huitieme ordre, il y a une huitieme 
constante, mais nous n’avons pas a nous en inquieter, car elle 
n’entre que dans la variable parasite p. CTest la constante K qui 
figure dans le second membre de Peq nation 

ri2 v -+- #1 = K. 

Quant aux autres constantes, ce sont celles du Soleil, elles doivent 
£tre regard^es comme connues et non comme des constantes d’in- 
i^gration. Parmi les sept constantes que nous conservons, nous 
distinguerons les to et les trois autres E,, E 2 ) que nous 
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appellerons les (3. Je puis ecrire alors, en distinguant les coefficients 
de ces deux sortes de constantes, 


d f - d ( Q — uv )-*- ^ A dm + 2 B 


<l>i dt -h v dm* 


et en posant 


0'= 12 — uv. 


nous arrivons a la formule 


d4) 


dy" = A dm -4- B d$ — <$1 dt -h v dm* 


La fonction Q r se trouve alors definie a une fonction arbitraire pres 
des constantes [3 et try, par l’equation 


(i5) 




Le second membre est une fonction des constantes {3 et des argu- 
ments w, periodique par rapport aux w. 

Designons par H la valeur moyenne de cette fonction periodique, 
ce sera une fonction des [3 seulement. Nous pourrons alors sup- 
poser que 

= H t ~h~ Q ' ■) 

Q" elant une fonction des [3 et des op, periodique par rapport aux w. 
La valeur moyenne de cette fonction periodique peut £tre choisie 
arbitrairement puisque Q' n’est definie qu’a une fonction arbitraire 
pres des constantes, nous la supposerons nulle. Dans ces condi- 
tions Q f/ est fonction seulement des j3 et des w. 

Nous poserons alors comme au n° 129 

os; 2 ■ r ' d y - = 2 w< dwi +■ 2 Cf d?ih 


et nous reconnaitrons que les W et les G sont des fonctions des [3 
et des w, periodiques par rapport aux tv. Nous avons done Fidentit^ 
suivante en rapprochant les equations (i 4 ) et ( 16 ) 

(17) ^ W dw -+• ^ G — d(Hz)= ^ A dt-\-v dm*. 

Cette equation doit devenir une identity quand on j remplace Wi 
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par nit + Alors notre relation (17) peut s’ecrire 

2^ W ( n dt -+- t drt -f~ dm ) -+- C rf ( S — H dt — t d)$ 
= A dm -+- ^2 ® dft — <t> 1 dt -h v dm%, 


d’ou, en idenlifiant les coefficients des diverses differentielles, 


(18) 


i tisH-JWn, 

J Wf = hi 

( Wo = Ao -4- v, 

‘2 w $^ c - b 



Discutons les equations (18) et *appliquons-leur les lemmes des 
n °s X 08 et suivants. Si nous prenons d’abord 1 ’equation W/= A*, 
nous vovons que le premier membre doit etre une fonction des (3 
et des w : periodique par rapport aux op, et que le second doit etre 
une fonction des ft et des to. D’ailJeurs cette equation doit devenir 
une identite quand on y rein place w par nt to. Cela n’est possible 
que si W / et A/ sont feme lions des ft seulement (pour 1 = 1 , 3 ou 4 )« 
Si nous prenons maintenant la derniere equation (18), nous 
verrons que le premier membre se compose d’une fonction C p6- 
riodique des w, et d’une autre fonction periodique des w multi pliee 
par /; et le second membre d’une fonction B des (3 et des to et 
d’une autre fonction des ,3 multipliee par t. L’identite n’est possible 
que si 1’on a separement 

x'l dn d H 

2d^Tft = -dft’ 

G = B = fonction des ft. 

On aura done 

<19 ) ^ Wrfn = ctti, 


et, en rapprochant de la premiere Equation (18), 

( 20) — 2] n dW . 

Nous remarquerons que dans la somme dn, le terme W 2 d/i2 
P. - II (a). 


2 
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ne figure pas, puisque /? 2 est une constante donnee et que par 
consequent dn 2 == o. 

Nous avons dit que W l5 W 3 , W 4 dependent seulement des (3, 
voyons ce que nous pouvons dire de W 2 , nous avons Pequation 

W 2 = A s -t- p, 

ou W 2 est une fonction des (3 et des w periodique par rapport 
aux w, ou 

n 2 v = = K — <i> ]5 

K etant notre huitieme constante, tandis que est une fonction 
des (3 et des w, periodique par rapport aux w ; quant a A 2 , c’est une 
fonction des huit constantes [3, m et K. Nous avons done l’identite 

(21) - 4 - A 2 -H K, 

dont le premier membre est fonction des j3 et des w , periodique 
en (V, et le second membre fonction des J3, de K et des uj. L’iden- 
tite n’est possible que si les membres sont fonctions des [3 seule- 
ment. J’dcrirai alors la premiere equation ( 18 ) sous la forme 

(ti? W 2 *+* ri\ Wi H- /is W 3 -4- 714 W 4 = H, 

ou chaque terme du premier membre et le second membre sont 
fonctions des [3 seulement. 

Nous avons vu plus haut que quelques-unes de nos coordonn^es 
sont des fonctions paires et d’autres des fonctions impaires des 
arguments 

»ii w b 

Le dernier de ces arguments pp! { est une constante que l’on peut 
regarder comme donnee une fois pour toutes; nous pouvons la 
supposer nulle, ce qui revient a prendre pour l’axe des x K le grand 
axe de l’orbite terrestre. Dans ces conditions nos coordonnees 
seront des fonctions paires ou impaires des quatre arguments 

w i = &2= w’b = w' 2 . 

Les fonctions suivantes seront paires : 


y\ ? 

dx\-> 

~dt’ 

dy\ dy\ 

~dt' Uf’ 

da' 
dt ’ 

H, W, A, 

P, 71, 

a/g., X ’dy, 

dt dw 

d£t! 

dw 
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Les suivantes seront impaires : 


rl'r' 

1 

&21 y 1 ? y si 1 


Q', 0", C, B, x'y", 


dx\ dy\ 

dt 3 dt 3 

, </Q' 

* i $ 3 d$ * 


Ce que je voulais faire remarquer, c’est que C devant 4 tre inde- 
pendant des w et en m£me temps fonction impaire des w , devra 
£tre nulle, de sorte que nous aurons simplement 


2 dy" - do!' = 2 w*, 


et a cause des equations (18) et (21) 

( 22 ; 2 x ' d y n ~~~ dQ " = 2 ~ — J ~ » 


ou Ton a pose 

A/ = A/ = W/ (i = 1 , 3, 4)» 

A o = A 2 -4 ^ 

*2 

de telle facon que les A/ dependent seulement des [ 3 . 


319 . De la formule (22), on peut deduire une serie de formules 
importantes, que nous ^crirons sous la foi’me 


t , dy dil" 

A* dp ~‘~df ~°’ 

__ a' 


(i = I, 3, 4), 


= A' 




712 


et, en particulier, pour la constant^/w', 

—— - 

£d, hm ! dm' ~ 


Mais dans Papplication de cette derniere formule, il faut faire 
attention a une chose; nous avons £crit plus haut (p. 12) 

a), y” = niaf^m' ,a), 
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oCi y, = fL. Ilfautalors fairc attention que x' et/' dependent de m', 

CL 

non seulement directement, rnais par l’intermediaire de a qui est 
fonction de m' , et par Tintermediaire de a qui est egal a ^ 7 - 

320. Axes tournants. — On peut avoir un grand avantage a 
rapporter le systeme a des axes tournants auto or de Forigme avec 
une vitesse angulaire n 2 ; ces avantages ressortent deja de ce que 
nous avons dit plus haut au n° 313; nous avons vu en effet que 
lorsque E 3 =o, les equations du mouvement se presentent sous 
une forme particulierement simple en emplojant ce syst&me 
d’axes. 

II serait facile d’etablir les equations du mouvement par des 
procedes elementaires, mais il sera pr^f^rable de rattacher la 
transformation aux principes du Chapxtre I. Je reprends les 
equations ( 11 ) du n° 318; d’autre part, pour me rapprocher des 
notations de MM. Hill et Brown, je designe par x, y, z les coor- 
donnees de la Lune par rapport aux axes tournants, de telle sorte 
que Ton ait 

x = x\ cos u x\ sin u, 
y — — x\ sin u x f . 2 cos w, 

■5= *s- 

L’angle a est Tangle des axes mobiles avec les axes fixes; c’est 
done w 2 = w” 2 = n 2 1 q- rs 2 ; la lettre u a done la meme signification 
qu’au n° 318; je poserai ensuite 

y\ = X cos u — Y sin u, 
yl = X sinw -h Y cosu, 

et enfin 

v r =v — (Xj/* — Yx). 

Dans ces conditions on a 

dx = dx f i cos u -+- dx ' 2 sin u 4- y du , 
dy = — dx j sin u -+- dx\ cos u — x du, 

et, par consequent, 

^j^y" dx'-h p du = X dx -+- Y dy -+- Z dz -f- v r du. 
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Le changement de variables est done canonique et les equa- 
tions ( 1 1 ) ( p . id) deviennent 


( 23 ) 


dx 

dF' 

dX 

dF' 

dt 

“ dX’ 

dt 

~ dx 

du 

dF 1 

dy 

dF’ 

It 

~ do 1 ’ 

dt 

du 


avec les equations qu’on deduit des premieres par permutations 
circulates de x , y , z et de X, Y, Z. 

On a d’ailleurs 


F' = n.p'-f- n 2 (Xy — Yx), 



et les Equations ( 23 ) nous auraient donne 


dx _ 
dt ~ 


■«2 y . 


dy 

~dt = Y n * x i 


dz 

dt 


Z, 


dX 

dt 


_ _ d*\ 
~ y - ~ dx ’ 


dx 


da 

dt 


— Tl 2 , 


en posant 

U t + U 


^ = ~^r l H- n* ( Xy — Y *), *'i = + -+- ^X 2 =F'-»^' 

Remarquons que nous avons 

«s(X>-— Y:r)= ti 2 jk — — n? (a?2 + j/2), 

-(Sr- f •)* ft--^ 

/i 2 {Xjr — Ya?) = — 


<^ 4 ) 


La derniere des formules (24) nous fait connaitre le sens du terme 
complementaire /r 2 (Xjr — Y x); il represente, au facteur constant 
pres n 2 , la constante des aires dans le mouvement absolu. 

Dans le cas ou E 3 = o, F' ne depend plus de x, y, z , X, Y, Z, 
et pas de u. On retrouve done l’integrale de Jacobi qui peut 
s’ecrire 

F # — n 2 f>' = i ^X 2 ^’ ~~ Yar) = const., 
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ou 


bien 


Nous trouvons ensuite 

'j? d y dstf =2^ d x — ( X *X“~ Yar)rfw, 

2 s*’ = 2 X 

£x'dy = 2 l a;dX.-t-(Xy — Y*) da, 
et, en rapprochant de la formule (22) et remarquant que u = <v 2y 
(* 5 ) ^dX-da’^S.dm- *'' dwt 


d’ou 


Al 2 


en posant 


: $>!-+- 7 l 2 (X?' — Ya?) = F 7 — 7 l 2 t>'. 


D’autre part Q ;/ est encore ddfini par 


(26) 


dO? , v 


dt 


-H, 


H etant une conslante choisie de telle fagon que la valeur mojenne 
du second membre soit nulle. II suffit en effet de se reporter a la 
formule (i 5 ) et de remarquer que 


2 


a; 7 


d£ 

dt 


== 2 a7 w Ya7 )* 


321 . Choix des constantes. — II importe de comparer aux nota- 
tions prec^dentes, celles qui ont ete employees par Hansen, par 
Delaunay et par Brown. Delaunay emploie les arguments suivants : 

D s=s w\ — w'i =z — w 2 = distance moyenne Lune-Soleil, 

F = w\ ->r w 2 = Wi -f- = distance moyenne Lune-noeud, 

l = w\ ■+■ w\ = Wi •+- w 3 = anomalie moyenne Lune, 
l' = w\ -h w 3 = anomalie moyenne Soleil; 

ou si <*>3 = o, comme nous l’avons suppose, /' = <*>" = w> 2 = u. 
Hansen prend pour arguments : 

g — HPj-h w 3l ' g’ — w 2 , 

0) SSS fP* -h <P S; C» — fP 4 . 
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Brown adopte les arguments de Delaunay, mais dans ses formules 
figurent plus habituellement les exponentielles imaginaires 

Z = e**> 9 £* = e /F , e u . 

11 faut aussi comparer comment sont notees les constantes qui 
correspondent a celles que nous avons designees plus haul par a , 
m\ E, Eo, E 3 . 

La constante E 3 , excentricite du Soleil, est partout designee 
par e'\ le rapport des deux moyens mouvements que nous avons 
appele ml est designe par Delaunay par m. Au contraire, ce que 
Brown designe par m , c’est le rapport 

/i 2 __ in r _ mouvement sideral du Soleil 

iii — /i 2 _ 1 — in ' mouvement synodique de la Lune 

INous poserons comme lui 

m 

il est aise de voir que m! peut se developper suivant les puissances 
de m et que, par consequent, toutes nos series qui procedent suivant 
les puissances de m 1 peuvent etre develop p^es suivant les puissances 
de m ; la convergence s*en trouve m&me augmentee, on verra plus 
loin pourquoi. 

Les constantes qui correspondent a E* sont designees par 
Delaunay et Brown par a, e ; mais elles sont definies d’une maniere 
differente. Pour Brown a est le coefficient de £ = e* D dans le deve- 
loppement de x iy, et par consequent celui de cosD dans celui 
de x , ou plut6t a est defini de fagon qu : il en soit ainsi, si Ton ne- 
glige E,, E 2 et E 3 ; pour Delaunay, a est defini par l’egalite 

m.\ -+* m<i = n\ a 3 , 

ce qui vaut mieux. 

Delaunay definit e de telle fagon que le terme principal de 
^equation du centre ait meme expression que dans le mouvement 
keplerien; pour Brown, au contraire, e est le coefficient de a sin/ 
dans le developpement de x\ suhv* — «r' 2 cosW|. Ici la difference 
est importante puisque le e de Brown est k peu pres le double de 
celui de Delaunay. 

La constante d’inclinaison qui correspond a E 3 est designee 
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par *v par Delaunay et definie de telle facon que 1 expression du 
terme principal de la latitude soit la meme que dans le mouvement 
elliptique. Elle est designee par E par Brown et definie comme le 
coefficient de ar/ sinF dans le de\ eloppement de >3. 

On voit que toutes les definitions de Delaunay font jouer le pre- 
mier role aux coordonnees polaires, et celles de Brown aux coor- 
donnees rectangulaires. Ges differences n’ont aucune importance, 
et l’on po urrait faire varierles definitions de bien d autres manieres 
sans rien alterer d’essentiel. 


322. II importe de se rendre compte de la grandeur de ces difife- 
rentes constantes. On a sensiblement 


i3 


, I t r 1 

6o 1 20 


e 


— (Brown) ou — (Delaunay), 
io 20 v 


i 



Nos series procedent suivant les puissances de ces diverses quan- 
tites, mais il importe de remarquer que dans le coefficient d’un 
meme cosinus, ou d’un meme sinus, l’exposant de e, e f , K, a est 
toujours de meme parite, de sorte qu’en realite nos series procedent 
suivant les puissances de 


= ou 

IOO 


I 

4oo ’ 




4oo 


36oo 


1 6 o ooo 


Aussi la convergence sera-t-elle beaucoup plus rapide par rap- 
port aux excentricites et aux inclinaisons que par rapport a m qui 
joue le rdle du coefficient p dans la thdorie des plan&tes exposde 
au Tome 1. G’est le contraire de ce qui arrive dans le cas des pla- 
netes. Aussi convient-il de d^velopper, non pas d’abord par rap- 
port a p, puis par rapport aux excentricites, mais au contraire 
d’abord par rapport aux excentricites et ensuite par rapport a m. 
G’est la la difference essentielle entre la theorie de la Dune et celle 
des planetes. 
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323. Nous commencerons par determiner les termes qui sont 
d’ordre zero , par rapport aux excentricit£s, a Finclinaison et a la 
parallaxe, c’est-a-dire par rapport aux E et a a. Ces termes ne 
peuvent dependre de w \ , (v 3 ; ils seront done des termes en 

cos t j 7 

. /^ir/qW’2 . 

sin 7 

Si nous adoptons les variables z du n° 320, nous verrons 

d’abord que s est nul puisque l’inclinaison est supposde nulle; 
d’autre part, d’apres le n° 192, t. I, on a 

k x = -* 2 . 

Ensuite, d’apres le n° 190, t. I, x et Y ne contiennent que des 
cosinus, tandis que y et X ne contiennent que des sinus. 

Enfin, d’apres le n° 317, l’exposant de a etant nul, /c 2 doit etre 
impair. En resume le terme general dans x et dansj^ sera respec- 
tivement en 

cos(a£ -+- i)< w>i— Wi), sin(2A* -+- i)(w\ — w> 2 ), 

ou cos (2 A: -h i)D, sin( 2 A + i)D, en introduisanl l’argument D de 
Delaunay (c/n° 321). 

Pour former les equations du probl&me, il faul reprendre les 
equations generales dtablies au n° 320 et y faire a = E 3 = o. Faire 
a = o, e’est negliger la parallaxe, c’est-a-dire reduire U 3 a son 
premier terme, ce qui donne 
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P 2 etant Je polynome de Legendre 

P — 3 cqs 2 y — i 


On a d’ailleurs, puisque s = o, 

AC 2 =,r 2 -hj' 2 , 


et puisque BD est parallele a 1’axe des x (l’excentricite E 3 <Stant 
nulle et par consequent le mouvement de B an tour de D circulaire 
et uniforme), puisque par consequent y est l’angle de AC avec 
Faxe des x 1 

AC cosy = x t 

d’ou, finalement, 


P 2 AC 2 = 


zx* — y* 
2 


Mais nous avons trouve plus haut au n° 320 

Ui -4- U3 /v v x 

F — = 1 h n z (Xy — Yx). 


Rappelons d’autre part que 


U, 


m.\ 4- m 7 


d’ou 


d’ou enfin 


AC 

U 3 _ — .r 2 

— — n i - , 


a'* 


, X 2 - 4 - Y 2 m x -4- m 1 . 2X Z — y 2 /v v . 
F = - L_J _ ^2 -4-7i 2 (X^ — Yx). 


Comme F ; ne depend pas de a , la variable auxiliaire v f se rdduit 
a une constanle et ne joue aucun role, et les Equations cano- 
niques (23) du Chapitre precedent deviennent 


(1) 


dx 

dt 


= X -4- 7i*y, 


dy 

dt 


= Y 712 X, 


1 dX 
j dt 

| €L 

\ dt 


TYl\ 7W-7 

AC* 


x - 4 - zn\x n- n* Y, 


77l\ -4- /n 7 

ac» y 


ru\y — /z 2 X. 
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MaJs les deux premieres equations nous montrent que 

dx 


*7 


d\ __ d l x dy dY _ d-y 
~dt ~ IF ~ ni ~di' ~dt ~ dt 1 


dt 


dK v </2.r 

-dF- n * Y =dF- 

d\ v 

_ + / , 2 x = -^. 


■ 2 7lo 


dt 

dx 


n\ x, 




de sorte que les dernieres equations (i) deviennent 


(2) 


d-x 

dt 2 


— 2 3 n?, x -+- 

at 


m i -f- ni-i 

\C* X 


= o. 



2 7^2 




111 j -4~ 

AC* 


y = o. 


Ces equations peuvent etre mises sous une autre forme; si nous 
posons 

w'\ — w\ = D = T, 


cela revient a changer l’unite de temps de telle sorte que 1’argu- 
ment de Delaunay joue le role de temps. On a alors, d’apres 


le n° 321, 


dx 

dt 


jii — n 2 = 


71 0 

rn 


Si alors nous posons 


x 


les Equations deviendront 


Til i -4- />? 7 
~1 ~ ^2?’ 


(3) 


d*x 

~d? 


dy 
' dx 


I d*y 
[ dx 2 


dx 
1 dx 


3 171“ X -+- 


XX 



Telle est la forme particulierement simple que prennent les 
Equations du mouvement de la Lune quand on neglige : 

i ° La parallaxe ; 

2 ° L’excentricit^ solaire ; 

3° L’inclinaison. 


Le probleme que nous proposons maintenant est de trouver une 
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solution particuliere de cette equation; a savoir celle qui corres- 
pond au cas oil la constante E< est nulle. Comme nous venons de 
voir que x ety sont developpables suivant les cosinus et les sinus 
des multiples de 2 D, ou de 2 -, nous voyons que cette solution 
particuliere est une solution perioclique. Le probleme a ete 
entierement r^solu par Hill dans un Memoire de tout premier 
ordre dontnous allons exposer les principaux resultats ( American 
Journal of Mathematics, t. 1). 


324. Equations homogenes. — Ainsi que nous l’avons vu au 
Chapitre precedent ces equations admettent l’int^grale de Jacobi, 
qui s’^crit 




r 


C’est d’ailleurs une combinaison immediate des equations (3). 
Ici r = fx*-+-y 2 designe la distance AC. 

Nous avons done trois equations que je puis ecrire en mettant 
dec d 2 x 

x‘ et x ,r pour^— et ce qui n 5 a plus d’inconvenient; nous obte- 
nons ainsi 


( 4 ) 


[ x " — ■ imy ' — 3 m 2 x H- — = o, 


xr 

7T 


< y 4 - 2 mx 


_ c. 

2 2 r 


Entre ces trois equations nous allons eliminer x; nous trouverons 

i xx ff -)-yy"-h ^m^yx' — xy*)->r - (#' 2 h-j >/ 2 ) — - m - ~ = G, 

( 5 ) 2 2 

( yx ,r — xy” — 2 m (yy r 4 - xx' ) — 3 m 2 xy = 0 . 

La premiere s’obtient en multipliant les trois equations (4) par 
x, y et 1 , et la seconde en les multipliant par^, — x et o et ajou- 
tant. 

On voit tout de suite que les premiers membres des Equations 
(5) sont des polynomes homogenes du second degre par rapport 
a x, y et leurs derivees. 
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32o. Equations imaginaires. — Hill met encore ces equations 
sous une autre forme en introduisant les notations 

u = x 4- iy , s — x — 

Elies deviennent alors 

l us" 4 - u" s — 2 mi( us ' — s u j 4- u! s 

(6) j 

f us" — u" s — 2 mi( us' 4- sm') 

ces deux equations peuvent £tre remplacees par une seule 
us' 

(7) us" — imius' h — — (i 5 u-- 3 !- 4- 3 s 2 ) = C, 

a laquelle il convient d’adjoindre son imaginaire conjugu^e 

( 8 ) u"s 4- 2 mis a' 4- — — (1 3 s 2 4- 1 8 us 4- 3 u- ) = G. 

20 

II est aise de verifier en effet que les deux equations (6) ne sont 
autre chose que la somme et la difference des equations (7) et (8) 

326. Les equations (3) forment un systeme du quatrieme ordre, 
elles conviennent pour l’etude des termes de degre zero et pour 
les termes qui dependent seulement de l’excentricite lineaire E*, 
mais sont independants de l’inclinaison E 2 , dela parallaxe a, et de 
l’excentricite solaire E 3 . L’excentricite lunaire E* est l’une de nos 
constantes d’integration ; si nous la supposons nulle, il ne subsistera 
que les termes de degre zero, que nous nous proposons actuelle- 
ment d’etudier. L’ensemble de ces termes de degre zero represente 
done une solution particuliere de nos equations (3); comme ces 
termes dependent d’un seul argument D = t, ils sont periodiques. 

Le probleme revient done a chercher une solution periodique 
des equations (3). 

Si Ton construit la trajectoire T du point x , r par rapport aux 
axes tournants, cette trajectoire sera une courbe fermee, puisque 
x et y sont des fonctions periodiques du temps. Comme x ne eon- 
tient que des cosinus, et y seulement des sinus; comme d’autre 
part les deveioppements de x et de y ne contiennent que des 


~(« -t-s ) 2 = iG, 
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termes dependant des multiples impairs de *7, les deux axes des x 
et des y sei'ont des axes de sym^trie pour cette courbe ferm£e T. 

En diminant x entre les Equations ( 4 ) et en regardant C comme 
une constante arbitraire, nous formons un syst&me ( 5 ) ou (6), ou 
(7) et(8) qui peut £tre regarde comme plus glnlral, puisque ce 
systeme ne change pas quand on change #, y et C en \y, }, 2 C, 
en ddsignant par \ une constante quelconque. Si une courbe 
fermde T satisfait aux Equations ( 3 ), elle satisfera 6galement aux 
Equations (7) et (8); mais ces Equations (7) et (8) admettront en 
outre pour solution toute courbe homothdique a T par rapport a 
Forigine. 

Le probleme a et^ entierement r^solu par Hill, ainsi que je l’ai 
rappels dans mes Methodes nouvelles de la Mtcanique celeste 
(t. I, p. 97). Pour les petites valeurs de /?2, Hill trouve une courbe 
ayant la forme generate d’une ellipse, pour 772 = 1,78, il trouve 
une courbe avec deux points de rebroussement situ^s sur Faxe 
des x. 

S’ilavait pouss£ plus loin, il aurait sans doute lrouv6 une courbe 
avec deux points doubles symdriquement situ^s sur Faxe des x; 
puis ces deux points doubles se seraient confondus en un seul situd 
k Forigine, puis ils auraient disparu et la courbe T serai t redevenue 
une courbe ferm£e sans point double mais parcourue dans le sens 
retrograde. Mais les premieres courbes determin^es par Hill ont 
seules de Fintdr£t pour l’etude du mouvement de noire satellite. 

327 . Calcul des coefficients. — Void k quoi revient la m^thode 
de M. Hill, pour les petites valeurs de m. Aulieu des Equations ( 3 ) 
envisageons les Equations plus g^n^raJes 

I g > 3 4 3 0 xfl? 

I or — % py p*x m*x h — r = 0. 

, . 1 2 2 r * 

{3 bls) j 3 3 xy 

( /+ipx’—-pty+. -m'y- i- -L =0. 

En les traitant comme plus haut on en d&Iuirait, au lieu de (7) 
et (8) les Equations, 

(?bis) us^—ipius'^r — 7 p-us = ^-(i5w*-h 3 s*)h-G, 
240 

(8&iV) 2pisu'^r — 3«*-+-i5s*)-+- C. 

*•*4 o * 
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Nous allons developper la solution suivant les puissances 
croissantes de m 2 ; une fois la solu tion obtenue, nous y ferons p = n ? 
pour relomber sur les equations (3). 

Supposons le probleme r^solu et posons 

u * 

£ = e iz , - = w 0 -+- m 2 Ui ■+• m*9u q -b .. . , 

s 

— = so -h m-Si m+s -2 m*<7s q 

C = Go "+• #? 2 Cj -r~ M** Go TTl-'l Qq -f- . . . . 



II s’agit de determiner par approximations successives, d’a- 
bord u 0l s 0 , Co, puis G|, puis w 2 , s*, C* ; . . . , et ainsi de 

suite. 

On prendra d’abord 


Uq — S Q — — - I ^ 



Quand on aura remplace dans (y bis) et (8 bis ) u , 5 et C par 
leurs developpements ( 9 ), on trouvera dans le premier membre 
des termes de la forme 

m 23H-2p Ka 

et d’autres dune forme analogue ou u a ou sp, outousdeux, ontete 
remplaces par Pune de leurs deux premieres d^rivees. Dans le 
second membre on trouvera des termes de la forme 

m 23H-2p-H2£2 

ou de la forme 
ou encore de la forme 

m 2 aG a . 


Je suppose qu’on ait d£ termini dans les approximations prec£- 
dentes 


«0, 

u u 

. jusqu’a 

u q . 

-li 

SQ t 


. . . , » 

s q- 

-1? 

Go, 

Ci, 

» 

Cq- 

-li 


et qu’on veuille determiner u q: s q: G r Egalons les coefficients 
de m-i. 

Dans le premier membre nous devrons retenir les termes qui 
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contiennent m 2 ? en facteurs, c’est-a-dire ceux ou 

Ces termes ne contiendront que des quantites connues sauf ceux 
ou a = y, jj = o, ou hien encore ceux ou a = o, (3 = q. Dans le 
second membre nous devrons retenir ceux ou 

a -h p = q — i, 

ils ne contiendront que des quantites connues; nous devons encore 
retenir le terme m 2 ? G q qui est inconnu. 

Les termes encore inconnus du premier membre de (7 bis) 
s'ohtiendront done de la facon suivante; prenons par exemple le 
terme 

us\ 

je puis l’ecrire 

us"= u[s^ Zr l Y= U [($?)" — **(*?)'■+■ (*?)]• 

Or, on a, par les developpements (9), 

uXr'- = *£m**n. a . 

Nous devons, d 5 apr£s ce qui precede, retenir les termes tels 
que a-h [3 = nous distinguons parmi eux ceux qui contiennent 
une des inconnues s q ou it q \ c ? est~a-dire les termes tels que a — o, 
|3- = #, ou bien a =<7, j 3 = o. Les termes a conserver sont done 
les suivants (en rappelant que u 0 = 1 , s 0 = 1 , s 0 = o, s' 0 = 0); 

S q 2 IS q S q lL q . 

Nous opererions de la m£me facon sur les autres termes du 
premier membre etil nous resterait comme termes, contenantl’une 
des inconnues, 



Conservons ces termes (10) dans le premier membre, et faisons 
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passer au contraire dans le deuxieme memLre terns les termes qui 
ne contiennentque des quantiles connues; nous obtiendrons une 
equation de la forme 

(i i) A (Sq^ tty) — ^7 *4- G57. 


Dans cette equation, A(^, u q ) n’est autre chose que l’expres- 
sion (10); c’est done une combinaison lineaire a coefficients con- 
stants de s q , Uq et de leurs derivees. La fonction <f> q est Pensemble 
des termes connus qui se trouvaient dans le deuxieme membre ou 
qu'on j a fait passer. C’est une fonction periodique de z. 

En operant de la meme fagon sur (8 bis) on aurait obtenu 

(12) A '(s q , u q ) = & q -\-G q ; 

A l (s q , iig) est V expression 



imaginaire conjuguee de A (s q , u q ) et d&luite par consequent 
de A (s q , u q ) en pennutant s q et u q et changeant i en — i. 

<&' q est l’imaginaire conjugee de $ q ; c’est done une fonction 
connue et periodique de t, d^veloppable suivant les puissances 
positives et negatives de £ 2 . Pour passer de a ii suffit de 
changer £ en et de remplacer les coefficients numeriques par 
leurs imaginaires conjuguees; nous verrons plus loin que ces 
coefficients numeriques sont toujours reels et par suite que cette 
derniere operation est inutile. 

Quoi qu’il en soit, nos fonctions inconnues se trouvent deter- 
minees par les equations (11) et (12) qui forment un systeme 
d’equations lineaires a coefficients constants et a second membre. 

Soient, par exemple, a, a 1 , ? et ^ les coefficients de dans 
u q et Sq] il s’agit de determiner les coefficients inconnus £ et r* 
a 1 ’aide des coefficients connus a et b. 

Pour cela, les equations (1 1) et (12) nous donnent 


P. - II (2). 


4-2jG>4“ 

4- 

3 
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Ces deux equations du premier degre nous donneront £ et tj. 
Le cas de k = o m^rite une mention speciale; d’abord les pre- 
miers membres des deux Equations (i3) deviennent identiques; 
ensuite il fauttenir compte dans les seconds membres du terme C q ; 
les Equations s’ecrivent alors 

04) < ' 4 

I + + = 

Elies ne peuvent £tre satisfaites que si a = oI\ or a et a 1 sont 
imaginaires conjnguees par definition; les Equations ne peuvent 
donc£tre satisfaites que si a est r^el ; nous verrons plus loin qu’il 
en est toujours ainsi. 

Les equations (i4) nous donneront done pour £ -H '0 une valeur • 
r^elle A; on prendra £ = '/j = de telle fagon que £ et */] soient 

imaginaires conjugu^es, et en m&me temps reelles. 

La presence de C q introduit une indetermination dans le pro- 
bl&me. On peut supposer que la constante C est une donnee 
de la question; alors G ne depend pas de m et 1’on doit faire 
dans (i4) C q = o. 

On peut ^galement s’imposer la condition que le coefficient 
de ^ dans u et celui de dans s soient <%aux a i; Ce coefficient 
ne dependant pas de m, on doit avoir i* = 7) = o, de sorte que la 
premiere Equation (i4) se reduit a 

O — CL *4— G qi _ 

ce qui determine G^. 

328. Je dis d’abord que les coefficients de tous les termes du 
d^veloppement de u et de s seront r^els. En efFet, supposons que 
cela soit vrai pour 


HQ) IL\) . . , , Ug—\) 
s Qi • • • 7 £<7—1 j 

je dis que cela sera vrai pour u q et s q . 

En efFet, si cela est vrai pour 

<10; 


( a <7)1 
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cela sera vrai egalement pour 

(16) £u' x , is'v u’v s'i (*<gh 

et par consequent pour les termes connus de(^ bis) qui sonttous 
le produit de deux expressions de la forme (10) ou (16) affectes 
d’un coefficient reel. 

Done, dans le d^veloppement de 0 ><j suivant les puissances de s 2 7 
tous les coefficients sont reels. Done, dans les equations (i 3 ) 
et (i4), les quantites a et a r sont reelles; or les coefficients de ces 
equations du premier degre (i 3 ) et (i 4 ) sont Egalement r^els; 
done nous en deduirons pour £ et r h e’est-a-dire pour les coef- 
ficients de u q et s g , des valeurs reelles. c. q. f. n. 

Je dis maintenant que nous aurons seulement dans u 0 et s 0 des 
termes en £°; dans« 4 et s, des termes en 

S 2 , S- 2 ; 

dans u 2 et s 2 des termes en 

ft 5°, 5" 4 ; 

dans u z et s z des termes en 

C 6 , £~ 2 , S~ 6 ; 

et ainsi de suite. En d’autres termes, je dis que u'C i ~ { et st sont 
developpables suivant les puissances de m 2 £ 2 et/?z 2 s -2 * 

Je dis en effet que, si cela est vrai pour les premieres approxi- 
mations, cela sera vrai aussi pour l’approximation suivante. Je 
suppose done que 

«cti s a (a <q) 

sont des poljnomes homogenes de degr6 a en £ 2 et £“ 2 , et je me 
propose de demontrer que u q , s q seront aussi des poljnomes ho- 
mogenes de degr£ q en £ 2 , D’abord les d^rivees u etc., seront 
comme u ay etc., homogenes de degre a. Considerons maintenant les 
differents termes de nous avons d’abord : 

i° Ceux des termes du premier membre de (7 bis) qu’on a fait 
passer dans le deuxi&me membre parce qu’ils ne contenaient pas 
de quantity inconnue. Ils sont £gaux a un facteur num&nque 
pres k 

(«<£, ?<?; a + P — yh 
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U(L ou s$ p ouvant £tre remplaces par une de leurs derivees ; ce sont 
done des polynomes homogenes de degre q en £ 2 et £~ 2 ; 

2 0 Les termes provenant de llULJL - ils sont de la forme 

? 2 u «Kp (a -+■ p = q— i), 

et par consequent homogenes de degre q en £ 2 , ? 2 ; 

3° Les termes provenant de — - g -~; ils sont de la forme 

(a H- P = ^ — i) ; 

ils sont done aussi homogenes de degre q. 

Done % est un polynome homogene de degr<$ q en CS 
Mais u q et s q se deduisent de par les Equations (i3); les termes 
de u q et s q correspondent a ceux de $> q et contiennent les m£mes 
puissances de done u q et s q sont les polynomes homogenes de 
degrd q en £~ 2 . c * Q- F ‘ D * 

II resulte de la que, si Ton d^veloppe le coefficient de £ 2 * ou 'Q~* k 
suivant les puissances de m : le developpement contiendra seule- 
ment des termes ou l’exposant de m prendra les valeurs 

ik , 2*4-4, 2*4-8, 

Ce developpement proc^dera done non pas suivant les puissances 
de m , mais suivant celles de ttz 4 . C’est ce qui explique la conver- 
gence extremement rapide du precede de M. Hill , chaque ap- 
proximation donnant 4 ou 5 decimales exactes nouvelles. 

329. Je dis maintenant que les coefficients £ et yj calculus d’apres 
les procedds pr^eddents sont des fonctions rationnelles de p. 
Supposons en efiet que cela soit vrai aux approximations ante- 
rieures; je dis que cela est encore vrai a Tapproximation suivante. 
En eflet, cela sera vrai des coefficients de ® q formas par 
multiplication en partant de ceux de u a , s a (<x<<7); oela sera 
d-onc vrai des coefficients a et a 1 qui figurent dans les equations ( 1 3) ; 
les coefficients de ces equations (i3) £tant rationnels en p, il en 
sera de m£me des inconnues £ et 7), c’est-&-dire des coefficients 
de u q et s q . c. q. f. d. 
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Quels seront les facteurs des denominateurs de ces functions 
rationnelles? II est aise de les determiner. En efFet, la resolution 
des equations (i3) introduit au denominateur un facleur qui estle 
determinant de ces equations, lequel est egal a 


~(/? 2 — 4/2 — a -+- < zk 2 ). 

Les facteurs du denominate ur seront done les poljnomes 

p - — 4 p — 2 + 2 /c 2 , 

oil Ton devra donner a k les valeurs paires successives 

1, 4 , 6 , — 

ce qui donneles polynomes 

/ i? 2 - 4/2 - 4 - 6, 

(17) | p 2 — 4/2 4 - 3 o, 

ip 2 — 4/24-70, .... 

Si done on developpe le coefficient de suivant les puis- 

sances de m, le coefficient de chaque puissance est une fonction 
rationnelle de p dont le denominateur est un produit de facteurs 
de la forme (17), chacun de ces facteurs pouvant etre eleve a une 
puissance superieure a 1 . 11 ne faudrait pas en conclure pourtant 
que le coefficient de £ 2A ’ H ? considere comme fonction de pe t de m } 
est une fonction meromorphe de ces quantiles, ou encore qu’il se 
reduit a une fonction meromorphe de m quand on fait p = m. 

329 bis. Quoi qu’il en soit, Papproximation par ces methodesest 
extrdmement rapide. M. Hill a calcule les coefficients ou plutdt 
leurs rapports avec 10 ddcimales; la premiere approximation lui 
donnait gendralement 6 d^cimales exactes, la seconde 1 1 et la 
troisieme i5. D’autre part, la d^croissance des coefficients est aussi 
tres rapide; dans le ddveloppement de u : par exemple, les coeffi- 
cients de C 3 ) • diminuent tres vite, chacun d’eux dtant 

environ la trois-centieme partie du precedent. Le coefficient de £ 
qui est naturellement le plus grand nous donne le terme principal 
de l’orbite nontroublde; viennent ensuite ceux de et de qui 
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correspondent a l’inegalite considerable connue depuis la fin du 
moyen age sous ienom de variation . 

Les valeurs numeriques sont done tres exactement connues; 
mais nous ne possedons encore queles rapports de nos coefficients ; 
les equations (7 bis) et (8 bis) contiennent en effet une indeter- 
minee G, et nous avons profite de cette inde termination pour sup- 
poser que le coefficient de £ dans u ) de meme que celui de 1 
dans est egala 1. Gela revenait a satisfaire aux equations (i 4 ) en 
faisant 

t ■= 'q = O, CL Cq — 0. 

G’est ainsi que nous avons opere a la fin du n° 327 . Nous avons 
ainsi trouve une solution particuli&re des equations (7 bis) 
et (8 bis). 

» ass 4 , (t)» ^ = +1 (C 

caracterisee par ce fait que, A et etant imaginaires conjuguees, 
le coefficient de £ dans A (Q est egal a 1. Alors nous aurons une 
autre solution en faisant 

w = s = ao^i (C)j 

a la condition de changer l’indeterminee C en ci\ C. II reste a de- 
terminer a 0 . 

Pour cela, il faut (apres avoir fait p = m) revenir aux equations 
primitives, qui peuvent s’ecrire 

w h- 2 mi u ~~ - m 2 (u-hs) H — — = o, 

et voir quelle valeur de a 0 correspond k la valeur donnee de x. 
Pour cela soit 

il viendra pour x = o, £ = 1 : 

r = u = s = a 0 ^a k ~ ao + O), 
v! = + etk, u" — — a 0 ^(2/c -H i ) 2 ajc . 

Les coefficients a* etant connus, on connait done facile- 
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ment A (i), <!/(i), et il vient 

(18) «o+ 2 (0[^(0 + aw^'(i)~ 3 m^(i)] = x, 

d’ou l’on tirera sans peine a Q . 

330 . On voit facilement que la resolution de liquation (18) 
entraine F extraction d’uneracine cubique. Nous devons done nous 
attendre a trouver que, dans le voisinage de eertaines valeurs sin- 
gulieres m 0? P expression de #0 et par consequent celle des autres 
coefficients de u qui sont a Q ou, . contiennent en fac- 

—i 2 

teur ( tyi — m 0 ) 3 ou (m — m 0 ) 3 , P ar exemple. 

C’est en effet ce qui arrive pour m 0 = 1 . 

Supposons done qu’on veuille developper a 0 suivant les puis- 
sances de m\ alors la convergence sera comparable k celle d’une 
progression geometrique de raison 

m 

m 0 etant le point singulier le plus rapproche; si ce point est i r 

puisque m = *^> cette raison sera 

m _ 1 
m {) ~~ 12 

Si Ton avait d^veloppe suivant les puissances de 

. Tio m 

m = — = > 

ni i-h tn 

comme Fa fait Delaunay, la raison aurait 

m! __ 2 _ 
m 0 """ 1 3 ? 

done notablement plus grande. 

On arrive au m<*me resultat en ce qui concerne le ddveloppement 
de a t ; le point singulier le plus rapprochd provient avant tout de 
la presence du facteur p*— bp + § dans nos denominateurs 
(c/. n° 328 in fine) ; dans ce facteur on fait p — m, de sorte que 
l’on a a peu pres 


m .\ — 4 ni o *4-6 = 0, 



40 CHAPITRE XXV. 

d’ou 

[ m 0 1 = v/6 

et, pour la raison de la progression, 

12 /6* 

A 

= «^IL. 

12 /6 

On voit combien il est plus avantageux de choisir m comme pa~ 
rametre au lieu de m f . Naturellement la meme difference se 
retrouve dans le calcul des termes suivants, puisque les nouveaux 
d^veloppements que Ton obtient pour ces termes se deduisent de 
ceux que nous venons d’obtenir pour les termes d’ordre z£ro. 

330 bis. Une solution particuliere remarquable s’obtient en fai- 
sant x = o ; on trouve alors 

u = s = [3£ 4 '-l- 

a, [3 et k etant lies par les relations 

a ^ — X 2 — 2 pk — \p*) — ~ P — °? 
p X 2 -4- X — ^ ~ a /n 2 = o, 

d’ou 



ou, en faisant />=/??, 

X 4 — m 2 X 2 = o. 

Pour que cette solution nous convienne, il faut que k soitentier 
impair, et par consequent que l’on ait 


m 


avec on aurait eu 




I mo ! 


I “f“ T)Xq 


et, pour la raison, 


m = i , 


m = 3 


m = 5 , 
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La courbe decrite par le point y se reduit alors a une ellipse. 
Mais cela correspond a x = o, c’est-a-dire, en se reportant a 
liquation (18), a 

4 '*( I ) [<K'(0 -+- 'zrni'b'(i) — 3 m*b( i)] — o. 

Si alors on donne a x une valeur donnee finie, il faut faire a n — oo. 
Done, quand m se rapproclie de la valeur i, la trajectoire T du 
point x, y tend a £tre de plus en plus semblable a une ellipse, mais 
les dimensions absolues augmentent an dela de toute limite ; e’est 
de cette facon que s’introduisent dans a 0 des facteurs tels que 

_ i 

(m — i) 3 , 

ainsi que nous Pavons explique plus.haut. 

L’ellipse a laquelle on parvient en faisant m = 3 , m = 5 , etc., 
ne fait pas partie de la serie de solutions que nous envisageons ici, 

puisque la periode n’est plus arc, mais 
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331. Nous allons maintenant determiner les termes du premier 
degre par rapport a l’inclinaison E 2 ; nous negligerons done comme 
dans le Chapitre precedent la parallaxe et Texcentricite solaireE s ; 
seulement nous ne supposerons plus z — o. 

Dans ces conditions, nous avons encore 


^ =-^3 4 P 2 AG*, 

m' /y/3 


mais on a 
d’ou 

D’ailleurs 

d’ou 


p 2 — ^-P °- S — (. f, AG.cosY = a7; 


AC 2 sa r 2 = x * -+- y* -h z*, 

U s _ 2 — y 1 — Z 1 

m\ ~~ n 2 ± 


Ui 


nx\ -4- m-i 


F'= /ijp'-j- 


X*H-Y 2 -{-Z 2 mi -h nii 


•7z 2 (X r - Yar) 


v . , ‘2a? 2 — r 2 — 

— Y x) — n\ 


Cela pose, nous retrouvons les equations (r) du Chapitre prece- 
dent, avec cette difference que AC = r ne represente plus \]x 2 
mais z 2 . 

Nous trouvons ensuile 

dz _ d¥' _ 

~dt ~ ~dL ~ 

d*z dZ dF f mi -+- m 1 

~di* - dt ~~ n * z ^ *• 
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Si nous posons, comme au Chapitre precedent, 


T = fl% ) t. 


n 2 


Til — /i 2 


' ” dz* 


m t -f- m 1 
(ni — nif ’ 


liquation pr^cedente devient 

(i) = o. 


avec 


e = /??*-+- — • 

,.3 


Si nous cherchons les termes de l’ordre de Pinclinaison, nous 
negligeons le carre de Pinclinaison et par consequent s 2 ; nous 
pouvons alors prendre 

r = \J x t -t- jk 2 , 

de sorte que nous retombons sur les equations (i) du Chapitre 
precedent, sans changement. Et comme nous devons faire E< = o, 
puisque les termes que nous voulons calculer sont independants 
de Pexcentricite lunaire E <? nous devons prendre pour x et y la 
solution particuliere etudiee dans le Chapitre precedent. Done x 
et y sont des fonctions connues du temps, developpables suivant 
les puissances impaires positives et negatives de £ = e ix . 

^ nous est ensuite donne par Finlegrale de Jacobi 

x, __ a?' 2 -hj y' 2 3 mix’ 1 ^ 
r ~ 2 2 Jy 


d’oh Fon deduit ais^ment et 0. 

r 8 

Done 0 est une fonction connue du temps d£veloppable suivant 
les puissances paires, positives et negatives de 

Nous avons vu que, quand on introduisait les deuxparametres/? 
et m qui figurent aux Equations ( 3 bis), (7 bis) et (8 bis) du Cha- 
pitre precedent, 

«C _1 = iy), 

= C ( rr — iy) 

sont developpables suivant les puissances de p , m 2 £“ 2 ; il 
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en est done de meme de 0, de sorte qu’apres qu’on aura fait p = m , 
0 sera d^veloppable suivant les puissances de 

772 , t 2 , m 2 t~ 2 . 

De plus, quand on change t en — t, ou £ en « se change 
en s, de sorte que r et 0 ne changent pas. 

Si done nous posons 

8 

on aura 

©A = 3 

et l’on pourra ecrire 

@ S3 © 0 4- 2 01 COS 2 Z -+- 2 02 COS 4 T ■+• 

0 ^tant developpable suivant les puissances de m , /n 2 ^ 2 , /n 2 ^” 2 , 
le coefficient 0* contient en facteur m 2k , de sorte que les coeffi- 
cients 0 decroissent tres rapidement. 

332. II s’agit done d’int^grer liquation (i). 

Cette equation est de m£me forme que celle que nous avons exa- 
minee au Chapitre XVII des Methodes nowelles de La Meca - 
nigue cileste. Nous savons done : 

i° Que cette equation admet deux solutions particulieres de la 
forme 

(a) z e= 3 = t) , 

^(t) etant une fonction periodique de la forme 

ou k prend toutes les valeurs entieres positives et negatives. 

Nous devons remarquer en effet que liquation (i) ne change pas 
quand on change t en — ?; l’existence de la premiere des solu- 
tions ( 2 ) entraine done celle de la seconde. 

2 0 Nous aurons ensuite les deux solutions 


^ = F('u) = + — t), 

* “/('O = XjV*rc^(T) — 
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dont la premiere est une fonction paire et la deuxieme une fonc- 
tion impaire de v. 

Nous pouvons acheverde determiner (qui n’est encore de- 
termine qu’a un facteur constant pres) et la constante a, de telle 
sorte que 

F (o ) = i, F'(o) = o, 
f(o) = o, /(o)=i, 

c’est-a-dire 

3° En vertu d’un th^oreme ddmontre dans mon Memoire Sur 
les groupes des equations lineaires ( Acta mathematica , t. IV, 
p. 212 ) et rappele dans Les methodes nouvelles de la Mecanique 
celeste (t. II, Chap. XVII, p. 23o), la fonction F(?) consideree 
comme fonction des 0* est une fonction entiere. 

Nous trouverons alors 

F(T) 

k variant de — cc a -+• oc. 

La fonction 'i('r) est periodique de periode ”, de sorte que 
Ton a 

(3) 6(x) = ^(°) = F(ie) = COS^-TT. 

333. Liquation F(tc) = cos^t: a une grande importance; c’est 
elle en effet qui va nous permettre de determiner g-, c’est-a-dire le 
mouvement du noeud. 

En effet, la solution gdnerale de liquation (i) peut s’e'crire 
z = Ej F ( *? -+- Wit ), 

E 2 et tsj 4 <5tant deux constantes d’integration, la premiere reprdsen- 
tantla constante d’inclinaison, et la deuxieme la constante tu 4 qui 
figure dans la formule du Chapitre XXIV, 

w k =«»J=n 4 i-r w*. 

Si nous faisons en effet g = ^ on voitquela formule pr£- 
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cedente peut s’ecrire 
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Z =^E 2 &*COS(c*^-!- Wi-h 2/:Wi — 

ou nous reconnaissons la forme des formules du Chapitre XXIV . (II 

IT * 

suffiraitde changer laconstante arbitraire enu 4 +-> pour avoir 
un sinus au lieu d’un cosinus.) 

Cela pose, les passages au noeud s’obtiendront en faisant z — o ; 
le terme principal du developpement Stant 

E 2 6()Cos«>4, 

les passages au noeud ascendant auront lieu sensiblement pour 

v 

Wi H- WP4 = gt -H TU4 = - • 

Dans l’intervalle de deux passages au noeud consScutifs, la diffe- 
rence des longitudes moyennes de la Lune et du Soleil aura aug- 
ments sensiblement de — ; si 1’on considere n ■+■ 1 passages conse- 
secutifs au noeud ascendant, cette difference aura augments 
sensiblement de et cela d’autant plus exactement que n sera 
plus grand. Done g mesure le moyen mouvement du noeud. 


334. Pour calculer F(t), voici comment nous allons procSder. 
Posons 


©o= 2^ 


etecrivons l’Squation ( 1 ) sous la forme 

( 4 ; jBr+q*z=s(e 0 — e)*. 


puis cherchons a dSvelopper F(t) suivant les puissances de 0 4 , 
etc.; le dSveloppement sera tres convergent, puisque nous avons 
vu que ces coefficients sont tres petits, et d’autre part que F(t) cst 
une fonction entiere de ces coefficients. 

ReprSsentons alors par z v , . . ., z k les termes de F(?) qui 
sont de degrS 0 , 1 , . . . , k par rapport aux coefficients 0< , © 2 , . . . ; 
nous aurons, pour dSterminer successivement 


•® 0 > ^17 « • • j 


* • r 
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le systeme d’equations 

*o-*-7**o = o, 
z]-h q^Zi = (e 0 — Q)z 0 , 



*£4- 7 2 **= (0 o — 0)-sa-i* 

Ce sont des Equations lineaires a coefficients constants et a 
second membre, imm^diateraent intdgrables; il faul, pour achever 
de definir z 0 , z<, . . . , se donner les conditions initiales qui 
seront 

*0=1, ^ 0 = 0 , Z\ = z\ = o, z k =z' k = o. 

Nous trouvons d’abord 



- *0 = cos q t, 

et pour Zk une expression ou figurent des termes en 

( 6 ) cos + ^F^smO/ -h ?)t, t^cos (2/ h- 2)^, 

j et [jl etant entiers et 2 pi i etant au plus £gal a k . 

En effet, il est a.ls6 de verifier que, si cela est vrai pour cela 
sera vrai ^galement pour 

Ainsi le terme general du d^veloppement de F(t) sera de la 
forme (6); il est aise d’en apercevoir la raison : nous avons trouve 

F(x) =^^yCos(^ + 

D’autre part, g est developpable suivant les puissances de 0<, 
0 2 , ... et se reduit a q pour 0 { = 0 2 = . . . = o. 

Je puis done dcrire 

F(t) =2*y' cos r(? + V> + {g — q )t] 

et d^velopper suivant les puissances de g — q\ soit alors 

cos, = ^ sin , =2 

les coefficients otp, et (3^ ayant les valeurs numeriques bien coDnues; 
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| F(-w) = ^bj^g - <7)*f^cos(<7 + *j)t 

I — ^ pjite — sin(? -h V>* 

On peut ensuite developper les bj et les puissances d eg q 
suivant les puissances de @i, ©o •• - 5 et ^ on devra retrouver le 
m£me developpement que par le mojen des equations ( 5 ). On voit 
que le terme general est bien de la forme (6). 

333. Voyons comment on peut se servir de ces developpements 
pour determiner g et les coefficients bj. 

Supposons que Ton ait determine 

- 0 , • • • 7 Z k 9 

on a done, a des quantity pres de l’ordre de /ft 2 ** 2 , 

F (t ) = 

et F(t) se presente sous la forme d’un developpement (7) dont 
tous les termes sont de la forme (6)* 

Dans ce developpement, conservons seulement les termes p£rio- 
diques purs en cos (q -f- 2 j) t, ceux ou t ne sort pas des termes tri- 
gonom6triques ; les coefficients de ces termes seront precis^ment 
les coefficients cherches, au m£me degr£ d’approximation, e’est- 
a-dire a des quantiles pres de Pordre de m 2k+2 ; et en effet a 0 est 
< 5 gal a 1 . 

Pour determiner g nous calculerons F (tc) en faisant % = dans 
le developpement (7) et nous aurons ensuite g par liquation 

F (tz) = COS 71 . 

En faisant cette substitution on trouve 

F(ir) = Hcos^ir -H H'singrrc, 

H et H 7 etant des fonctions entieres de ©< , © 2 , 

D’autre part, les coefficients de ces fonctions entires seront des 
fonctions rationnelles de q \ car, si les coefficients de Zk-\ dependent 
rationnellementde q , il ensera de m 4 me, en vertu des Equations ( 5 ), 
des coefficients de z&. 
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JLes premiers termes sont 

r '] 

* L 2 (! — ^ 2 ) 2 J 

r — 7T@f (l5q'+— 35<7 2 -r- 8) 770ri 

" Sia?TC L 4 ^( 7 ^) + 6 W-/)H 4 -^ |- 

Quand on change Tenx+^ 0 <, 0 3 , 0 5 , ... changent de signe ; 

la valeur de cos^-tt ne doit pas changer. Done, dans le develop- 
pement de F(tc), 0 ( , 0 3 , 0 5 , . . . entreront avec un exposant pair. 
Supposons un instant 

o = ©i = ©3 = ©5 = . . . ; 

alors 0 admettra la periode ^ ; alors, quand on changera xenx+^ 

0 2 , 0 6 , © 10 , 

changeront de signe; done 0 2 , 0 O , 0 1O , ... ne pourront figurer 
qu’a un degrt§ pair, a moins d’etre multiplies par 

@17 © 3 > **M ©1©37 ©1© 5 , .**• 

De meme 

©4, ©is, ©20 j ••• 

ne pourront figurer qu’a un degre pair, a moins d’etre multiplies 
par 

©17 ©37 ©1 ©3? ©27 ©63 *•*? ©2 ©6 j •••• 

336 . Method© de Hill. — Hill ramene le probleme a la resolu- 
tion d’une infinite d’^quations du premier degr£ a une infinite 
d’inconnues. On pourrait se demander si cette m^thode est suffi- 
samment rigoureuse; elle peutetre completement justifiee, maisje 
renverrai pour cette justification aux Methodes nouvelles de la 
Mecanique celeste (t. II, Chap. XVII, p. 260 et suiv.). 

Faisons dans l’equation (1) 

z = e^ x ^(t) 

il viendra 

— 2 b f ,( -ik -+- gy xyk+s 0 y b^v-^t+g = o 


p. — 11 (•?) 


4 
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ou, eu egalant ]es coefficients de K 2 ^, 

(8) [e»— (a/ -»-<?•)«] bj-h^Qj-kbk- o- 

On obtient ainsi une infinite d’equations lindaires entre les in- 
connues bj; pour les determiner, considdrons les expressions 
lineaires 

k 

Ces expressions doivent s’annuler pour \ = g- Formons leur 
determinant et appelons-le □ (5). 

Je suppose pourcela que l’on donne ay et a k toutes lesvaleurs 
entieres positives et negatives depuis — p jusqu’a H- p inclusive- 
ment. Nous aurons alors un determinant fini de 2p-M lignes et 
de i colonnes. Alors □ (5) sera par definition la limite dece 
determinant quand p croit indefiniment*. 

Ce determinant converge; nous pouvons remarquer en effet que 
les elements de la diagonale principale sont egaux a i. 

Dans ces conditions, on a 

I A I <e**\ 

en designant par a les divers elements du determinant autres que 
ceux de la diagonale principale. Le determinant convergera done 

avec^|a|; e’estee que demontreraient les considerations expo- 
sees au Tome II des Methodes nouvelles ( loc . cil .). Or ici 

Or cette serie converge absolument et uniformement, sauf si 
Qo = ( V ■+■ £) 2 

ou, puisque 0 O = si 

S = - V =*= 2 - 

Done □ (5) = lim A existe ; e’est une fonction de £ qui ne change 
pas quand on change ? en £ -p 2 ou en — 5 ; cela revient en effet a 
changer Pordre des lignes et des colonnes de notre determinant 
d’ordre infini, soit en faisant avancer toutes les lignes et toutes les 
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colonnes d’un rang (si l’on change £ en £ 2 ), soit en retournant 

le determinant (si Ton change £ en — £); on aura done 

□ 

□ (£) est done une fonction meromorphe de £, qui a pour poles 
simples 

£= — 2 jdriq. 


Je dis pour pdles simples; en effet, les elements de la f° me ligne 
seuls deviennent infinis pour £ = — 2 j ±q et infinis du premier 
■ordre. 

Or un terme quelconque de notre determinant ne peut contenir 
•en facteur qu’un seul element de la j^ me ligne. 

Quand la partie imaginaire de £ tend vers 1’infini, tous les ele- 
ments en dehors de la diagonale principale tendent vers o; done 


Soit 


limfH(£) = 1 . 


F(£) = D(0 


cost:£ — cos Tzq 

COSTCg — COS7C g* 


Cette fonction est meromorphe; elle ne pourrait devenir infinie 
que pour 

5 = -* 


qui rend □(£) infini, mais qui annule egalement cos7t£ — cost iq, 
ou pour 

£ V 


qui annule le denominateur costt£ — costt g, mais qui annule ega- 
lement □ (£) puisque nous avons vu que 

□ ($0 = o 

et par consequent 

□ (dz^) = D(- 2 / ±*) = o. 

La fonction F(£) est done entiere, mais elle tend vers 1 quand 
la partie imaginaire de £ tend vers Pinfini ; elle se r£duil done a 
i’unite, de sorte qu ? on a 

□ (£)= COS7Z \~ 

x cost:? — costt^ 
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On peut en deduire plusieurs manieres de calculer g; par 
exemple, a l’aide de liquation 

(9) cosit^* = I — 2 sin* ^n(o). 

337. U est ais£ de voir quelle est la forme de notre determinant. 
Supposons qu’un terme du developpement contienne comme fac- 
teur Pel^ment de la a iemc ligne et de la b [ °'" 0 colonne et, par con- 
sequent, il devra contenir un element de la b xhme ligne, soit 

celui de la c'* mo colonne et, par consequent, il devra contenir 

ensuite Telement de la c' k '" G ligne etdelad^ me colonne, c’est-a-dire 
® c _d, etainside suite jusqu’ace qu’on retombe surla a [bme colonne, 
ce qui arrivera, par exemple, si Ton trouve un element de 
la d i6 '" e ligne et de la a lhn>e colonne dependant de 

Done le terme envisage contiendra en facteur le produit 

— b &b—c ® c—d ® d—a 1 

ou un produit analogue, ou plusieurs produits analogues. 

Dans to us les cas la somme algebrique des indices doit 6tre 
nulle . 

Si Ton observe que 

e y = 0-j, 

et que par suite on ne fasse pas attention au signe des indices, nous 
dirons que les indices, toils regardes comme positifs, doivent 
pou voir se diviser en deux classes de telle fagon que la somme des 
indices de la premiere classe soit egale a la somme des indices de 
la deuxi&me classe. On pourra avoir par exemple des termes en 

Le coefficient de chaque terme se presente sous la forme d’une 
serie infinie, mais toutes ces series peuvent £tre somm^es k l’aide 
de la formule 

fio) ^ — - — = T:cota7Tc. 

v 1 JUx^r n 

Considdrons, en effet, un terme du developpement,* ce sera le 
produit de certains elements du determinant; parmi ces elements, 
les uns appartiendronta la diagonale principale et nous n’avons pas 
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a nous en occuper, puisqu’ils sont egaux a i ; les autres n’appartien- 
dront pas a cette diagonale, et ils seront en nombre fini (sans quoi 
notre tcrme contenant en facteur m a une puissance infinie serait 
infiniment petit). 

Appelons [f, £*] V element de la f i6me ligne et de la k i{ime colonne, 
et supposons par exemple que le terme envisage soit le produitdes 
elements 


[<*, b] [b, c] [c, d] [d, a] [a', b'] [ b c'][c\ a']. 

II contiendra alors en facteur, comme nous venons de le voir, le 
produit * 

(i i) Q a -b&b— c Q c —tf(dd— a 

et le coefficient de ce produit (n) sera une fonction rationnelle 
de £ que j’appelle R(£). 

Considerons maintenant le terme 

[a-h/i, b n].. .[d-+- n, a + n][a'+ n, b' r f- n ]. . ,[d' -+- n, a'-h n \ ; 

il contiendra egalement en facteur le produit (n), mais avec le 
coefficient R(£ 4- zn). 

Le coefficient definilif sera done 

y 272 ’)* 

Pour evaluer cette sorame, decomposons la fonction ration- 
nelle R(£) en Elements simples : 

d’ou 

2a<s+.»>- f 22 n 4= J - 

La sommation par rapport a n peut se faire par le moyen de 
Tequation (io). Si dans la decomposition de R(£) figurait un 
terme 

A 

(£-«)/>’ 
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on obliendrait la sommation par rapport a n de 


i 


A 

(£ -h 2 n — a )i> ’ 


en difierentiant l’equation (to). 

En fait notre fonction. R(£) n’admettra pour p61es, comme nous- 
Pavons vu plus haut, que £ = — 2 / :± <7 et ces poles seront simples m r 
Papplication de la formule (io) introduira done seulement 


cot 7r = ±: cot ~(t±.cj) 


Cela s’applique au cas ou nous n’aurions en facteurs qu’un seul 
produit de la forme 

®b~c @r— d ®d—a* 


Dans le cas ou l’on aurait en facteur deux (on plusieurs) produils 
de cette forme, comme il arrive, par exemple, dans le cas du pro- 
duit (i i), les choses seraient un peu plus compliqu^es. 

Supposons que nous ayons a envisager le produit des elements 

[n, n h- a][n -+- a, + n] [n\ n! -+■ a'] [n'~h a\ /i'], 

et que nous fassions varier n el n\ en supposant que a, b , a' ne 
varient pas. 

Nous aurons alors partout en facteur le produit 

a © — b a’i 


qui ne dependra ni de n , ni de n'; ce produit sera multiple par un 
coefficient qui dependra de £, de. n et de n et qui sera de la forme 

R<6 + n)R'«-hn'), 

R et R' etant rationnels. Nous aurons done a calculer 


en donnant a n et a n f toutes les combinaisons de valeurs possibles, 
en excluant seulement celles pour lesquelles la difference n — n f 
prend certaines valeurs parliculieres; car une des quantiles n* 
n + a, /i H- ne doit pas £tre egale a une des quantity n\ n/-{~ o! , 
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Supposons que a,, a 2 , . . . , cl p soient les differentes valeurs que 
ne doit pas prendre n r — /? ; 1’expression a evaluer, 

^R(5 + n)R'(« + A 

pourra s’ecrire 

2R(?-t-«)2 R '^- + - n )-2 P(?+ra ’ *«)— •— 2 P(?4 ' re ’ **>» 


ou 

P(5, «/) = R(E)^ + 4 


Chacune des sommes qui figurent dans cette expression est de 
la forme ^R(£ + n) et peut s'evaluer comme nous venons de 
l’expliquer. 

On peut voir que □ (?) est de la forme 


(12) 


□ ({) = A-+-Bcot^(?-Hgr)-i-Ccot^(5— q) 
-l-D cot2(|-i-gr)cot^(5 — q), 


A, B, C, D etant developpables suivant les puissances de 
0 2 , . . , , et les coefficients de ce d^veloppement £tant des fonctions 
rationnelles de q. 

Nous avons trouv^, d’autre part, 

□ (£)_= COS Ttg — COSTCO 

^ cost:* — costt q 
et 

F(tc) = cosTzg = H cosyn H- H' sin^Tr, 

H et H ; ^tant comme A, B, C, D des series procedant suivant les 
puissances de 0*, © 2 > . . . ; et avec des coefficients rationnels en q . 
On a done 

_ />v ( 1 — H ) cosctt: — H' sind77t 

□ (£) = iH- — r 1 

v ' cost:* — cos izq 

II suffit pour retomber sur le developpement ( 12 ) de remplacer 
(cositJ — costt #), . qr > 


2 sm 




par 
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de facon a tout exprimer en fonction des lignes trigonometriques 
des deux angles “(? — #)• 

338. II reste a determiner les coefficients 6*. Pour cela reprenons 
Je determinant □(?) et remplacons-y dans la ligne de rang zero 

0_ 2 0-1 0 i 

par des indeterminees 

1 > ^ 0 ? 2 ? 1 * • * * 5 

le determinant continuera a converger pourvu que les indetermi- 
nees soient limitees (cf. Methodes nouvelles, t. II, Chap. XVII, 
p. 260). II sera de la forme 

A = . . .-+- B_! Bo^oH- hl^l -r- B 2 5?2-f“* • •, 

les B etant des fonctions de £, qui admettent les m^mes p61es 
simples que □ (£) sauf £ = ± <7; de telle sorte que 

By(C0S7r5 — COSTt^) 

est une fonction entire. 

Je dis que, pour £ = g*, on a B* = b&; il suffit pour cela de de~ 
montrer que les satisfont aux equations (8), c’est-a-dire que A 
s’annule quand ony remplace xj par 1 et Xk{j^k) par 

®j-k 

®o — ( a /-*-*) 1 * 

C’est ce qui arrive car, pour y ^ 0, on obtient ainsi un determinant 
ayant deux lignes identiques; et, pour j = o, on obtient le deter- 
minant □ (g) qui est nul. c. q. f. d. 

339. Nous pouvons nous rendre compte de l’ordre de grandeur 
des coefficients b *. Pour cela reprenons les equations (5) et cher- 
chons a former a l’aide de ces equations non plus la solution F(t), 
mais la solution 

z = e 1 '?* ^( t ) =2 

Nous observerons alors que @ 0 — €> est developpable suivantles 
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puissances de /?z, m 2 ^ 2 , m 2 ; et je dis qu’il en sera de m£me 
de 

Si l’on integrait les equations (a) par approximations successives, 
on verrait, par une analyse toute pareille a celle du n° 334, que ^ 
se presente sous la forme suivante, 

2 &/C.V.-V- &■**, 


de sorte que z^~? est developpable suivant les puissances de m , 
de t, de £ a etde ^ _2 ; je me propose d’etablir que est develop- 
pable suivant les puissances de m, t, m 2 ^ 2 , cela peut se 

demontrer par recurrence, car les Equations (5) montrent que, si 
cela est vrai pour Zk ~\ , cela le sera egalement pour z&. 

De cela il resulte que bk contient en facteur ?n* k , ce qui montre 
que les coefficients bk doivent decroitre rapidement. 

340. Nous avons vu que g definit le mouvement moyen du 
noeud, mais nous n’avons qu’une premiere approximation. 

En eflet nous avons neglige E,, E 8 , a et le carre de E 2 . 

Le veritable mouvement du noeud peut, d’apres le Cha- 
pitre XXIV, se developper suivant les puissances de 

Ef, EL EL a 2 , 

et g(n t — n 2 ) ne nous donne que les termes de degre zero de ce 
developpement. 
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341. Je suppose qu’on ait un systeme quelconque d'equations 
diflferentielles ; pour simplifier, je supposerai deux equations seu- 
lement et deux inconnues x el y, et j’dcrirai ces deux equations 
sous la forme 

( f (^, y, y« y\ a i ? a s) = o, 

j Fi(x, y , y, x",y n , *i» *2) = o, 

ou a, el a 2 represented deux param£tres tres petits. 

Je suppose qu’on sache trouper une solution particuliere S des 
equations ( 1 ) en supposanl a, = ou = o; je me propose d’en deduire, 
pour les petites valeurs de a, et de a 2 , toutes les solutions peu dif- 
ferentes de ]a solution S (il y en aura evidemment si a { et ou sont 
tres petits). 

Je me propose de developper ces solutions suivant les puissances 
des parametres a et de certaines constantes d’int^gration {3, , [3^ 
p 3 , [3 4 qui s’annulent pour S. 

Nous connaissons dej& les termes de degre zero de ce develop- 
pement et nous voulons determiner successivement les termes de 
degr£ 1 , de degr6 2 , etc. 

Soient x — x 0j y = y 0 la solution S et posons 

x “ ^ 0 F(a701 ^ y«» ys, o, o), 

cette d^rivee dtant calculee en regardant x 0i y 0 , x' 0 , y' 0 , x" o: y\ 
comme des variables ind^pendantes, 

dF v , dF 



MOUVEMENT DU PERIGEE. 5 $ 

Supposons qu’on ait trouve le developpement exact jusqu’aux 
termes da k'* me ordre inclusivement et soit 
x = x k , y = y k 

ce developpement; si nous substituons ce developpement a la place 
de x et de y dans les equations (i), les premiers membres devront 
s’annuler aux quantitds pres du ( k H- i) ieme ordre. Done 

F i(x k ,y k ) 

seront des fonctions connues dont le developpement commencera 
par des termes d 7 ordre k -r- i . 

Soient maintenant 

x = x k + &r, y = y k - h ty, 

et supposons que nous voulions calculer ox et oy jusqu’aux termes 
du (k + i) lkme ordre inclusivement et en negligeant ceux du 
(k H- 2 )’ eme ordre. 

F(a?, y) = F («■*-*- S x.yk^r^y) 

pourra se d^velopper par la formule de Taylor, sous la forme 


Les termes en ox 2 (de m£me que ceux en §x oy r : oxox r , etc.) 
sont negligeables, car &x- esl d’ordre zk + 2. 

Dans le coefficient de 8#, nous pouvons faire 

x-Xq, y~yo, a 1 =a 2 =o. 


d\* 


Cela donne une erreur du premier ordre dans le coefficient 
et, par consequent, une erreur d’ordre k ~h 2 dans le produit 


dF 


^ ox , puisque %x est d’ordre ^ -hi. Cela revient a remplacer 


par 


dF dF dF 
dx’ dy’ dx ' ’ 

X, Y, X', .... 


II reste done (avec ce degre d’approximation) 
F{x,y) = F (a?*, 
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en posant 

^ X = X 8a? -+- Y o^-h X' ox ' -+- Y' o /+X' fcp'-t- Y' Sy", 
de sorle que les Equations (i) deviennent 

j 2 X s^=- F 

(2) ^ 

I Xi 8a? = — F, (a?*, yjc) \ 

les seconds membres sont connus de m£me que les X, de sorteque 
les equations (2) sont des equations lineaires a second membre. 
Les premiers membres demeurent les mimes a toutes les ap- 
proximations . 

Supposons en particular qu’on veuille determiner les termes de 
degre 1, en supposant a* = a2=o; les seconds membres de- 
viennent alors (pour la premiere equation, par exemple) 

F(.r 0 ,^ 0 , ^' 0 ,ro» *2), 


ou , puisque a 4 = a 2 = o, 

F(# 0 , Jo, x f 0f y<n x" 0 ,yl o, o). 

Mais cette expression est nulle, puisque 
^ = ^0, y = JK 0 

est une solution des equations (i) pour a< — a 2 = o. 

Les equations (2) deviennent done 

[ Vx 3 * = o, 

(3) ^ 

( 2 Xiaa7==o * 

Ce sont des equations lineaires sans second membre . 

Mais on sait que, pour integrer des equations lineaires avec 
second membre, il suffit de savoir integrer les equations sans 
second membre; on n’a plus a effectuer ensuite que de simples 
quadratures. 

Done, quand on saura determiner les termes de degre o 
et ceux de degre 1, pour a< = a 2 =o [ e’est-a-dire quand on 
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saura integrer les equations (3)], on saura determiner par qua- 
dratures les termes de degre superdeur quels que soient a, 
et a 2 . 

Les equations (3) ont regu le nom d’ equations aux variations 
(cf. Methodes notice lies } t. I, Chap. IV). 

Dans le cas qui nous occupe, ce sont la parallaxe a et l’excentri- 
cite solaire E 3 qui jouent le role des parametres a, et a 2 ; le r6le 
des constantes ( S est jou6 par 

(4) £i e izs t, E 2 e iT *2, E 2 e-«*». 

Nous avons deja determine au Chapitre XXV les termes de 
degre o; il nous reste done a determiner les termes de degre i par 
rapport aux constantes (4), e’est-a-dire par rapport aE, el a E 2 , en 
supposant 

a = E 3 =s= o ; 

on n’aura plus ensuite a eflectuer que de simples quadratures. Au 
Chapitre XXVI, nous avons calcule les termes de degr£ i par rap- 
port a E 2 ; nous avons maintenant a calculer les termes de degre i 
par rapport a Ei . 

Tel est le principe qui va nous servir a calculer les differents 
termes de nos d^veloppements ; nous verrons dans les Chapitres 
suivants quelles modifications de detail il convient d’apporter a ce 
principe pour qu’il s’adapte parfaitement a notre objet. 

342. Voulant calculer les termes de degre i par rapport a E 4 , 
nous devons faire 

Eo = E 3 = a = o ; 

nous retombons done sur les equations du Chapitre XXV ; nous 
avons d'abord z = o, puis les equations (4) du n° 324 : 

x n — T.my 1 — 3 m^x ^ = 0 , 

| y n -\- nmx r = o. 

Nous avons trouve plus haut au Chapitre XXV une solution par- 
ticuliere de ces equations (5), soit 



x~x 0 , y=y 0 ; * 



CHVP1TRB XXVII. 


62 

il s'agitd’en trouver une solution plus approchee 
x = a? 0 -H jr=js 0 -h oy, 

en negligeant le carre de E t et par consequent celui de ox et 8y. 
Formons done les equations aux variations des equations (5); il 
viendra 

j ox" — 'i m 0 y — 3 ox -+- A Sa? •+* B By = o, 

<6) I Zy+imSx' ■+■ B Sx •+. C 8jk = O, 


en posant 



et en rempla$ant bien entendu dans A, B ? G, apres differentiation, 
x ety par x 0 et jv 

Les equations (6) sont des equations diffdrentielles lineaires 
en 8x et 8 y, puisque A, B, C sont des fonctions connues. Le 
svsteme forind de deux equations du deuxieme ordre est du qua- 
trieme ordre; il admettra done quatre solutions 'independanles. 
Deux de ces quatre solutions sont deja connues. En effet, les solu- 
tions du Ghapitre XXV, 


x~x 0 , y —yo, 

dependent en rdalite de deux constantes arbitraires; nous avons 
trouvd en effet pour x 0 une solution de la forme 

# 0 = ?(*c, rn) 


et de m^me pourjK 0 - Or 

•c = — 


m = 


n 2 

Tl\ — Tl 2 


La solution continuera a convenir si 1’on change t en t s, de 
sorte qu’il resle 

3\>=<? («!—«!) — 

L n l — n 2j 

avec deux constantes arbitraires s et /z*, ou bien 
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Nos equations aux variations (6) admeltront done comme solutions 
particulieres 


dx q 

<s 

dy Q 

ds * 

? y = 

ds ’ 

dx 0 


dy 0 

dm ’ 

oy — 

dm 3 


en prenant pour variables ind^pendantes t , e et m, ou bien en 
revenant aux variables independantes t et m, 


< 7 ) 


;i 2 

dx 0 


lift dy 0 


ni 

dz 9 

oy = — 

m dz 9 



dx Q dx 0 

<x 

~ d )' 0 

dy* 

m 

dz dm 

°y = - 

— — ...... — _j_ 

m dz 

dm 


Inutile d’aj outer que, dans la premiere ligne de (7), on pourrait 
supprimer le faeteur constant en n 2 et m, puisque les equations 
sont lin^aires. 

Connaissant deux solutions parLiculieres d’un systeme du qua- 
trieme ordre, on sait qu’on peut ramener ce systeme au deuxieme 
ordre; mais il vaut mieux operer autrement, parce que la seconde 
solution (7) n’est pas periodique. 

Partons done de l’int^grale de Jacobi 

2 2 7* * 


et formons-en Pequation aux variations. 

Si nous supposons 3C=o, de fagon a slimmer la seconde 
equation (7), il viendra 

x^ 0#' H-j ' l 0 By — 3 m i x 0 ox -h Y ^—‘ ox -+- 

r o r o 

ou, puisque ety 0 satisfont aux equations (5), 

8 ) x f 0 Bx'+y 0 S y—{x" Q — zmy 0 )Bx~ (yl -+■ 2 mx^By = 0. 

On verifiera sans peine que la premiere solution (7) saLisfait bien 

a (8). 

Combinons alors la premiere equation ( 6 ) avec ( 8 ) ; nous aurons 
un systeme ( 9 ) qui sera du troisieme ordre, etqui admettracomme 
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solution particuli&re 



°r 


_ dy_ o __ 


d-z 


y\ b- 


Grace a la connaissance de cette solution particuliere nous pou- 
vons ramener ie systeme au deuxieme ordre. Posons en effet 

( lx -h i 8jr ) = ( $ H- ir \ ) (x' 0 ■+■ y''o )» 


avec son imaginaire conjuguee, et prenons pour nouvelles va- 
riables £ et r t . Alors, nos Equations admettent comme solution 
particuliere 

ox = x' 0y oy=/o> 


d’ou 


l-i-iv = $- = i ? ? = = °; 


done | satisfait a une equation lineaire du troisieme ordre, qui 
admet pour solution particuliere i , et 7j a une equation du deuxieme 
ordre de la forme 


( i o ) t/h- H 7} # — f- Ktq — o; 

pour ramener cette equation a la forme canonique, posons 

p etant noire nouvelle inconnue, etep une fonction pdriodique de t, 
telle que 


Alors p satisfait a une equation du deuxieme ordre de la forme 
(II) p ff -t-0p = O, 

ou 0 est une fonction connue de t. 


343. II nous faut maintenant montrer que cette equation (n) 
est de m£me forme que l’equation (i) du Chapitre precedent, 
e’est-a-dire que 0 est une fonction paire de t, periodique de 
p&riode tc el toujours finie. 

Dans les Equations (6), les coefficients A, B, C sont des fonc- 
tions periodiques, de sorte que ces equations ne changent pas 
quand on change t en t -4- rc. Si done nous designons pour un 
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instant par 8,#, S^, , 7), ce que deviennent ox , or, 5, 7j quand 

on change z en w -f- tc, et si 8#, 5 y est une solution, il en sera de 
meme de 8^, Mais nous avons 

lx -rio/ = (5 + ir^ix'o -+- *y 0 ), 

avec sa conjuguee. Si dans celte Equation je change z en z -f- tz, 

Zx, s, •*), &o, y 0 

se changeront en 

Sir? -nii — a?' 0 , — y'o, 

d’ou 

Oja?H- JOjj' = — (?i-H iT],) (a?Q iy f 0 ), 


ce qui montre que, si q, 7] est une solution, il en est de meme 
de — , — r Ji et par consequent de $ ce qui veut dire que 

les equations en 5 et en 7) ne changent pas quand on change z 
en z + -re. Done, dans l’equation (to), H et K sont periodiques. 

Si 1 } on change t en — z et oy en — oy les equations (6) ne 
changent pas, car A, B, C se changent en A, — B et C. D’autre 
part, x f 0 et y' 0 se changent en — x' 0 et y' 0 . Done 


se change en 


Sh -”1 = 

ox — i 8 y 

K — Wo 


lx -+• ily 

*o + Wo 
»-£-+■ /r 4 . 


ce qui veut dire que £ change de signe et que r, ne change pas. 
L’equation (io) ne doit done pas changer, ce qui veut dire que K 
est une fonction paire et H une fonction impaire de z. 

Si nous reprenons liquation 


nous verrons que 


Hfp TT 

— -f- H = o, 


, = «-*/" d \ 


Or H est une fonction p^riodique et sa valeur mojenne est nulle 
puisque e’est une fonction impaire; done j* Hdz , et par conse- 
quent, cp est une fonction pdriodique et paire. On en conclura 
p. — ii(2). r> 
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que0 est une fonction periodique et paire. Nous verrons d’ailleurs 
bientot la facon de determiner completement cp. 

11 reste a montrer que 0 est toujours finie. Pour cela nous 
remarquerons que ox et 3 y sont finis ; done 

„ . 3 X-r-i%Y 

% -+• i-r t — 4 

r 0 -+- ly o 

ne pourrait devenir infini, t etant reel, que si 1’on avait a la fois 

ce qui ne peut arriver, les trajectoires fermees de la Lune (rappor- 
tees aux axes tournants), etudiees au Chapitre XXV, ne prdsen- 

tant de point de rebroussement que pour ~ = i , 78. 

II faudrait faire voir maintenant que cp et par consequent 0 sont 
finis; pour cela nous allons determiner tp, mais il est necessaire de 
reprendre la chose d’un peu plus haut. 


344 . II serait facile de determiner cp en faisant le calcul tout au 
long, mais il est plus instructif de proceder autrement. Les equa- 
tions (6) admettent quatre integrates lineairement independantes ; 
nous connaissons deja deux d’entre elles qui sont les integrates (7). 

Les deux autres, d’apres les proprietes generates des equations 
lineaires a coefficients periodiques, seront de la forme 


(12) 


| 3# -H i 3 y =2 ^ 1 , 

I 3a? — ioy =2 ctS*****^ 


et, ea changeant t en — r et 5 y en — 3 y, ce qui ne change pas les 
equations, 


(ia bis) 


| 8 x — i &y =2 b k Xr c ~ lk - 1 , 
I 3a? *+- i 8/ = c > Xr c ~- k ~ 1 . 


En additionnant ces solutions, on trouverait une solution reelle 

dx (£/,-*- C&) COS[> s -+- (^ + {2^ + 1)^], 

dy =2 ( b k —Ck ) sin[w 3 -h 
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les coefficients bk et Ck sont reels. Le nombre c nous fait connaitre 
le mouvement du perigee de meme que le nombre g nous faisait 
connaitre celui du nceud, la remarque du n° 340 restant appli- 
cable. On pose 

-+- W Z = CZ -f- s, 

e etant une constante arbitraire el w$ representant la quantite 

= n 3 t h- m 3 

d^finie au Chapitre XXIV avec les conditions 


-+- ri\ 

n i — IX 


e = m„, 


'h 


It i • — 71 2 


Entre les quatre solutions (7), (12) et (12 bis) existent certaines 
relations bilin&ires dont nous allons indiquer l’origine. 
Gonsiderons un systeme d’equations canoniques 

dx _ dF dy dF 
^ ' dt dy ’ dt ~ dx’ 


d’apres le theoreme du n° 16 , rappele au n° 318 , il existe une fonc- 
tion Q telle qu’on ait 


^ a? dy — dil ^ h.d% — Fdt, 


les A dependant seulement des constantes d’integration a; on aura 
done, par exemple, 




dQ 

d*i 

dH 

d%i 


At, 

A,. 


Si nous diff&rentions la premiere par rapport a a», la seconde 
par rapport a a, , il viendra 


^ dx dy 
JmddQLi d%\ 

2 dx dy 
dx 1 dat$ 



d i y _ dQ 

dx± dx% dx\ dx% 

d ± y _ dQ 

dxi d%% ~ dxi dx % 




+ 


dX\ 

<57 ’ 

d Aj 

dT t ’ 
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@u, en retranchant, 

/ dx dy dx dr \ _ d\ j d\ 2 

Jmd' dl.2 d%i dni d%*j dz 2 d%\ 

Comme les A ne dependent que des constantes d’integration, le 
second membre se reduit a une constante. 

Supposons qu’on forme les equations aux variations des equa- 
tions ( 1 3 ) ; nous obtiendrons deux solutions particulieres de ces 
equations en faisant 


^ dx 


=^. 

OX = — T 

CIOL i 

ojk = 

di\ 

„ dx 


dy 

OX = -7—7 

a* 2 

or = 

di % 9 


et nous obtiendrons ainsi toutes les solutions ind^pendantes de ces 
equations aux variations. Si done 

dx = jj, dy = 7j. dx = $*, dy = ij* 

sont deux solutions quelconques de ces equations, on aura entre 
elles la relation bilineaire 

2^*-^= const. 

Nous pouvons appliquer ces principes aux equations qui nous 
occupent, puisque les equations (4) et (6) du Chapitre XXV 
derivent directement des Equations canoniques (i). Les variables 
conjuguees sont alors 

(x,X), (y, Y) 

et leurs variations sont 

dx, dX = d — /i 2 dy = ( — n 2 ) ( dx f — m dy), 

dy, dY = d ^ -+- n 2 dx =(/ii — n i )(dy'-+- m dx). 

Soil alors 

dx = di X , dX = Si X, dy = Oiy, SY = S, Y 

une solution particuli&re des Equations aux variations, et repre- 
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sentons de m^rne par des indices 2 une seconde solution particu- 
liere; on aura 

(8i#8 2 X — 8iX8 2 ^)*+-(^i^^2 Y — 81 Y 827*)= const., 
ou, enremplacant oX, ol r par leurs valeurs et divisantpar n s — /z 2 , 

( (8j# 8 2 #' — 8 2 #8i x') 

4 ( -+-(SiyB 2 y r — Ziyoiy' l-tri m(diyS 2 07 — dzyoi#) = const. 

La constante du second ruembre est nulle si les deux solutions 
5 , et o 2 sont identiques ; elle est nulle encore si Ton combine une 
des solutions (7) avec une des solutions (12) ou (12 bis). Si Ton 
combine en effet la seconde solution (7) avec (12), le premier 
xnembre de (i4) ne contiendra que des termes en 

^rc-i-Ti 0ll £c+n ( n etant entier). 

Aucun de ces termes ne peut etre constant a moins d’etre nul. 

D’autre part, si nous considerons une solution S l7 la constante 
du second membre ne peut etre nulle, quelle que soitla solution o 2 ; 
sans cela on aurait entre les quatre quantiles 

3j;r, 8 i y. Si#', Si 7 ' 

quatre relations homogenes du premier degre dont le determinant 
n’est pas nul, et ces quatre quantites devraient s’annuler a la fois. 

Nous devrons done conclure que la constante du second membre 
n’est pas nulle quand on combine entre elles les deux solutions (7) 
ou les deux solutions (12), (12 bis). 

345 . Que devient la relation (14) quand on fait subir a nos 
equations les transformations du n° 342 ? II est clair que nous 
allons avoir une relation bilin^aire entre deux solutions quel- 
conques des equations transformees 

rlu 

**12 > Is? ^2* 

Cette relation bilineaire’ dev ra etre ^videmment satisfaite, la con- 
stante du second membre etant nulle quand l’une de ces deux 
solutions corresponds a la premiere solution (7), c’est-&-dire quand 
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ou bien 


& = *n = Vi = 

?2=i, £' 2 = r i2 = r/ 2 = o. 


Nous en devons conclure que £< et ne doivent pas figurer dans 
notre relation bilineaire; nous aurons done une relation bilineaire 
entre 

Si» tfn 

$2) r i2 * 


Si, dans la relation (i 4 ), nous prenons pour la solution 3 * la 
premiere solution (7), nous retomberons sur Fequation (8) deduite 
plus haul de Pintegrale de Jacobi. Transformons cette equation (8) 
en faisant 

8x i 8y = (? -h i7])(a?'-h *>'); 


elle deviendra 

(i5) 


il _ ^o/u- ri» 
2T, a?i s -*-yo* 


Si, en partant de cette equation (10), nous remplacons £' et £' 2 en 
fonctions de r H et r l2 , ii restera une simple relation bilineaire entre 


'nu r/ u r i2j TQi, 


qui devra £tre de la forme 


<K T )Ohil*-- r iaVi)= const., 


etant une fonction connue de t; il est d’ailleurs ais£ de voir 

que 


cette fonction <p ^tant celle du n° 342 . 

346 . Reprenons les Equations (6) et passons aux variables £ 
et nous aurons deux Equations lin^aires entre 

r, v, w u' f 71, 

ou £ ne figurera pas puisqu’elles doivent 6tre satisfaites pour 5 = 1 , 
71 “ o. 
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Multiplions-les par — v 0 eix r 0 et ajoutons de facon a eliminer q" ; 


11 restera 


remplacons ensuite en fonction de r, a l’aide de l’equation (i 5 ); 


on trouvera am si 


(16) =°* 

M etant une fonction connue de t; si nous comparons a l’equa- 
tion (10), nous trouverons 


H _ 2 f® £Vji2j!Z« 


vVo 2 +y 0 2 


Nous en conclurons d’abord que cp ne peut devenir ni nul ni infini, 
par consequent que p reste fini, et enfin que 0 (qu’il esl d’ailleurs 
aise de former) est toujours fini, car, si @ devenait infini, l’une des 
deux integrates de l’equation(i 1) devrait devenir infinie. 

Done l'equation (1 1) est de meme forme que [’equation (1) du 
Chapitre XXVI, et tout ce que nous avons dit dans ce Ghapitre 
devient applicable. On peutse servir en particulier du determinant 
de Hill pour calculerle mouvement du perigee. La seule difference 
e’est que 0j est notablement plus grand, et il en resulte deux 
choses : d’abord la convergence du developpement est moins 
rapide que pour le mouvement du nceud et e’est ce qui explique 
les circonstances qui avaient tant e tonne les mathematiciens du 
xvm e siecle; ensuite certaines in^galites ont un coefficient notable. 
C’est ainsi qu’a c6t6 des termes en b 0 et qui represented les 
termes principaux de F equation da centre, nous avons les termes 
en b_ K et c { qui nous donnent la grande inegalite connue sous le 
nom d ’evection. 


347 . On peut obtenir immediatement un systeme du second 
ordre auquel satisfont S# et 8y, je veux dire les solutions (12) 
et (12 bis) qui nous interessent. Reprenons en effet la relation 
bilindaire (i 4 )> et imaginons que SaJK (que nous designerons 
simplement par Zy en supprimant Findice 2) represented 
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Pune des solutions (12) ou (12 bis) et que o K x, o K y represented 
une des solutions (7); la constante du second membre sera nulle, 
mais les solutions (7) peuvent &tre regardees comme connues, cela 
va done nous donner une relation lineaire entre ox, 8 y, So/,* 3 y ! ; 
comme nous avons deux solutions (7), cela va nous donner un 
systeme de deux equations difFerentielles du premier ordre entre 
ox et By, 

Prenons d’abord la premiere solution (7), 


> ~ 0 _ rJ 

OlX X n , 


=/o> 


0^ = ^, gi y=y 0 ; 


nous trouverons 


(16) #' 0 8 x'-hy ' 0 l/— x ’ 0 dx —f 8 y -4- 2 m(y Q 8 x — x ' 0 By) = o. 

Prenons ensuite la seconde solution (7), 

^ t , dx 0 * x . dy {) 

Oir = __^ 0+ ._ 5 0 iy = -- y 0 + - JL , 

d’ou 

a,* — 2.4 + P—Zk, s,y=--yi-hp-^.. 

Nous trouverons j^en tenant compte dela relation (16), en vertu 
de laquelle les termes en se detruisentj 


I P Sx’+psy- p lx -Ply 
1 dm dm J dm dm 


ox h-xo °y , 


II s’agit d’integrer le systeme (16), (17). On pourrait chercher a 
le ramener a la forme canonique, ou bien en faciliter Pintegration 
par Pemploi de Partifice du n° 327. Reprenons les Equations (3 bis) 
de ce n° 327, que j’dcris 


(5 bis) 


x ff — 2 py* — - p*x • 


5 _ Y.X 

- m 2 x h 7=0, 

2 r-* 7 


y* % px r — - p-y — - m- y -+~ = o. 
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Nous pourrons remarquer que ces equations peuvent Etre de- 
duites d’equalions de forme canonique; il suffitde prendre comme 
variables conjuguEes 

X, Y, ar, y 

avec 


F = 


X 2 -r Y 2 771 ! m- 


-h n 2 (Xy — Yx)- 


■ 3 h* 


a? 2 — y 2 


de former les equations canoniques 

dx __ d¥ dX _ dF 

dt ~ dt ~~ dx 9 


d’eliminer X et Y, et de poser 


dz — ( /? j — n z )dt, jiz=zp(n x — /i 2 ), 

771 1 -f- 7777 
A ~ (ft I 7^2 ) 2 * 


k = 772 ( 71j — Tig 


Recommencons sur ces equations (5 to) les memes raisonnements 
que sur les equations (5). 

Les equations (5 to) admettront le systeme de solutions perio- 
diques formees au Chapitre XXV et que nous ecrirons 

x 0 ~q(t, p , 771), y 0 = o t (z. p, m). 

Nous formerons les equations aux variations analogues aux Equa- 
tions (6) en posant 

x = x 0 -h8x, y=y 0 +8y- 

Ces Equations admettraient des solutions analogues aux Equa- 
tions (7) qu’on obtiendrait en diSErentiant, par rapport aux deux 
constantes d’intEgration z et n i} l’expression 


d’ou les deux solutions particulieres 


CLX Q / v &Xq 

0^7= _ =(71i— 71 S ) — 


Sx : 


dz 

dxp dx 0 


n i dx 0 h dx 0 

(ni — /i 2 ) 2 dp (/lx — n*)* dm 1 
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ou, en supprimant des facteurs constants, 


dx o dx o 

~-£ +p -df 


dx o 
dm 


qui remplacent les solutions (7). 

Dans ces conditions, (16) et (17) deviennent 


(16 bis) 2^0 5;r ~“‘ r o 5 7)= °» 



On va ensuite chercher a developper suivant les puissances 
de m 2 , et il arrivera alors la meme circonstance signalee a la fin 
du n° 328 que le developpement du coefficient d 7 un meme terme 
procedera suivant les puissances lion de m 2 , mais de m k . 

Mais, pour m = o, on a simplement 

#0 = A cost, jKo= A sinT, 


A etant une fonction de p facile a former; il vient done 

A = A. cost, 

,f^= Re ;n 


x' 0 — — A sinx, 
dx 0 


p~4= Bcost ’ ^^ =BsinT ’ 


en posant pour abr^ger 
Soient alors 


B =/> 


dk 

dp 


F(3* 0 , /<>, Fi(^ 0 ,/o, m) 


les premiers membres des Equations (16 bis) et (17 bis ) que j’ecris 
sans mettre en evidence ni /?, ni les inconnues 8#, 8^ et leurs 
derives. Nous designerons de m6me par 

F( A cost, A sinT, o) 


ce que devient F(# 0 , y 0: m) quand on j remplace m par zero, x 0 
et y§ par leurs valeurs approchees A cost, Asm?, et par cons&- 
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/ ' dx^ dy o 

quent Xq , y 0 , 


• •• par les valeurs correspond antes. Alors 


F (£r 0 j yoj m) — F (A cost, A sin t, o) = — IV1, 
F i ( Xq , y 0 , m ) — Fi(A cost, A sin?, o) = — Mj 


contiendront m 2 en facteur. Les equations (16 bis) et (17 bis) 
pourront alors s : 6 crire 

( F (A cost, Asiiit, o)= M, 

' \ Fj( A cost, A sinT, o}= M 1} 


et elles pourront s’integrer par approximations successives; on 
negligera d’abord m 2 , de sorte que les seconds membres serontnuls. 
Ensuite on remplacera dans les seconds membres 8# et 3 y par 
leurs premieres valeurs approchees, de sorte que les seconds 
membres seront connus et que les equations (18) se presenteront 
sous la forme d’equations lineaires a second membre, et Ton 
obtiendra a l’aide de ces equations de nouvelles valeurs approchees, 
et ainsi de suite. 

On voit aisement que les premiers membres sont de la forme 

f A( — sinT cost G( cost -f- sinT By), 

\JQj | B( cost 3a?' -t- sinT 8y f )-h D( — sinT 8x -+- cost By'); 

A et B (d£ja d^finis) de meme que C et D sont des coefficients 
numeriques dependant de p et faciles a former. 

Les equations (18) etant des equations linbaires a second 
membre, il faudrait savoir integrer les equations linbaires sans 
second membre, c’est-a-dire les equations obtenues en egalant a 
■ zero les expressions (19). Or, si nous posons 

£ = f-i Bu = ZrHBx-hiBy), 
n — K $s =Z {Bx—iBy), 

ces Equations sans second membre prendront la forme 

6 f -4- = o, 

V+tE-*-* 1 ! - °i 

a, { 3 , y, 3 etant des coefficients constants; les Equations (18) 
prendront la forme 

H- a£ -h {2nq “ P, 

— Qi 
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P et Q etant des fonctions connues, equations qui s’integrent 
aisement. 

348 . On peut se servir du systeme (i 6 bis), (17 bis) et du pro- 
cede d’approximations successives que nous venons d’exposer pour 
la determination de c et des coefficients bk et ch . II suffit d’appli- 
quer les principes du n° 335 . 

Determinons, par exemple, par le procede precedent, la solution 
particuliere des Equations (16 bis), (17 ’bis), qui, pour x = o, donne 
comme valeurs initiales 

037 = i, By = o. 

On aura alors 

So? ~ i By = ^ £a *£ <m " 2Am " 1 
avec la condition 

I 

(que nous pouvons supposer remplie, puisque les rapports des 
coefficients bk et Ck sont seuls determines). 

La valeur de Zx pour t = tz sera alors 

COS CTZ, 

ce qui determine c. 
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349. La determination des termes d’ordre superieur se fera en 
appliquant les principes da n°341. Je suppose qu’on ait deter- 
mine nos coordonnees jusqu’aux termes du k ikme ordre inclusive- 
ment, et soil par exemple x = Xk la vaieur approchee ainsi obtenue ; 
posons x — Xk- 1-8.2; et cherchons a determiner ox jusqu’aux 
termes du (£-bi)^ me ordre inclusivement; nous serons amenes a 
former des Equations analogues aux equations (2) du n° 341; ce 
sont des equations lineaires a second membre. Les premiers 
membres sont toujours les memes, quel que soil k. Nous avons 
appris a integrer les equations sans second membre aux Cha- 
pitres XXVI et XXVII, en formant les termes du premier ordre. 
L’integration des equations a second membre, et par consequent la 
determination des termes d’ordre superieur, peut done s’operer par 
de simples quadratures. 

Toutefois une complication se presente : nous n’avons pas seu- 
lement trois fonctionsinconnues qui sont nos coordonnees x, y, z ; 
nous avons encore deux constantes inconnues g et c d’ou depen- 
dent les mouvements du nceud et du perigee. 

Ces deux constantes sont developpables suivant les puissances 
de 


Ef, Ef, Ef, «*. 


Nous n’avons determine jusqu’ici aux Chapitres XXVI et XXVII, 
ainsi que nous l’avons remarque au n° 340, que les premiers termes 
de ces developpements, ceux qui sont independants des Ef et 
de a 2 . 

Supposons alors qu’on ait determine g et c jusqu’aux termes 
du (k — i) l4me ordre inclusivement, etsoient^-i et ces valeurs 
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approchees; soient ensuite gk~\ -h 8 g, C/ui + Be des valeurs plus 
approchees jusqu’aux termes du k ihme ordre inclusivement; en 
meme temps que 3#, ly, 82 , 1 ] nous faut determiner og et 3 c. Nous 
ferons cette determination, comme on le verra dans la suite, de 
facon a faire disparaitre les termes seculaires. 

3o0. Rappelons en quoi consiste la m^thode de Lagrange pour 
l’integralion des Equations a second membre. Considerons pour 
fixer les id^es un sjsteme du troisieme ordre, forme de trois equa- 
tions du premier ordre. Soient done X, Y, Z trois combinaisons 
lineaires de y, z; soient A, B, C trois fonctions connues, x\ 
y’ , z r les d^rivees de x, y, z\ nos trois Equations pourront s’ecrire 

(1) 2/4- X = A> /+Y=B, -s'+Z= G. 

Supposons qu’on ait int^gre les equations sans second membre 

( 2 ) X =y-h Y = z'-+- Z = o, 

et soient 

x — xt, y~yu z~zi (t = 1 , 2 , 3 ) 
trois solutions ind^pendantes du sjsteme ( 2 ). Onposera 
x = X*a7 2 -+- 

y-^hyt, 3 

etl’on cherchera a determiner les nouvelles fonctions inconnues A/. 
Le sjsteme (i) deviendra 

( 3 ) 2 X ^= A > 2 ^= b > 2 v,3, ' =g ’ 

et Ton en deduira les X'- sous la forme 

V/= «/A 4 - 4 - y iC (i = i,2,3), 

les a*-, les et les y* dtant des fonctions connues de de sorte 
que la determination des X/ etpar consequent celle de j', s sont 
ramenees a des quadratures. Les fonctions a/, j3/, peuvent s 5 ob- 
tenir par Fintegration des equations du premier degre (3). 

Tel est le procede classique, mats il y a des cas ou quelques- 
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simplifications sont possibles. Supposons d’abord un systeme du 
second ordre au lieu du troisieme, de sorte que les equations (3) 
s’^crivent 


(3 bis ) 


a', x x -f- )4 x 2 = A, 
^yi-^^y* = B. 


Supposons qu’on ait entre les solutions du systeme sans second 
membre une relation bilineaire 


— JKi^a = H. 

H etant une fonction connue de t . On aura 
X't H = A.yt — B 


d’ou 

et de meme 


ai -ir 




$2 = 


Xi 


* Q2 “ H' 

Supposons par exemple qu’on ait une equation de la forme 
x* 0 x = B, 

ce qu’on peut remplacer par le systeme 

x'—y = o, y+6x= B; 

on aura alors 

x 1 y 2 —x 2 y i = i, 

d’ou 

x = Xjari -+- X 2 a?j, 

Xi = — x 2 B, V 2 = x t B, 


351. Supposons un systeme d’^quations canoniques 


(4) 


dx _ dF 

Hi “ 5x’ 


^X 

dt 


dF 

dx* 


avec les variables conjugu^es 


y, 

X, Y, Z. 


Soit X 0? ... une solution particuliere de ces Equations; 
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posons 

X=X 0 -~$X, 

et, negligeant les carres de ox, . .., formons les equations a ux 
variations des equations (4)5 (#)? (X), ... etant des fonctions 
lineaires des six variables x , X, . . . ; soienl 

(5) o#'-r- (#) = o, 3X'-h('X;=o, 

ces equations aux variations. Soient 

lx~x h o X = X* ... r/=*i,a.3, 4, 5,6) 

six solutions independantes des equations (5 ). Considdrons main- 
tenant les equations a second membre 


(6) OX -r* — X, 

oX'-i- (X) = A*, 

011 A, A* sont des fonctions connues. Posons 


X =2>>/X, ( 

d 5 o 11 


^ViXi- k, 

2Jv,X,= A.*, 

et 



V £ = a £ A-b B -f- y*C ■+• a? A* -4- £3? B * -+- y* G*. 

II s’agit de determiner les fonctions a, 

A cet effet, rappelons-nous qu’on a, d’apres le n° 344, 

(x £ X A — x/cAt) -4- (yiYk — -+* (s 4 -Za- — «aZ/) = const., 
ce que j’dcrirai simplement 
( 7) — x A Xi) = const. 

On peut choisir les solutions particulieres . . . du systeme (5) 
de telle fa§on que la constante du second membre soit egale k i 
pour 

i = i, k = 2 ; i = 3, £ = 4 ; i — 5 , A- = 6, 

et a z^ro dans tous les autres cas. 
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En se servant alors des equations ( 7 ), on trouve aisement 

a, = X 2 , Y 2 . Tl = Z 2 , 

aj = — x 2 , = — y-2, 7 l= — -2,* 

a 2 = — Xi, , 

*2 = X U , 

et de merae 


0t3 — X4, 0t4 — X3 , Otg — X C , OCjj — X 5 . 


352. Reprenons les equations canoniques (23) du n° 320: 

dx __ d¥^ dX ___ 

dt dX 1 dt dx 

Supposons qu’on aittrou\e une solution exacte jusqu’aux tertnes 
du A ,,dme ordre inclusivement (par rapport aux E eta a). Soit 

x — x ttl y=y/< , - = 2 /,, X = X*, 

celte solution. II convient, eomme je Fexpliquais au debut de ce 
Chapiti'e, de tenir compte des constantes c et g . Je suppose que 
nous possedions des valeurs approchees de ces constantes, 

C = Cj : —\ , g — g A— 1 , 

exactes jusqu’aux termes du (A* — i) ieme ordre inclusivement. Je dis 
d’abord que l’erreur commise sur les deux membres des equa- 
tions ( 8 ) est du (A* + i)^ me ordre. Cela est evident pour les seconds 
membres, puisque nous y substituons a la place des inconnues des 
valeurs exactes jusqu’au A ,1 * me ordre exclusivement. Pour les pre- 
miers membres, ou figurent les constantes c et£*, cela exige un peu 
plus d’attention. 

En efiet x, par exemple, est une fonction periodique des quatre 
arguments Wi =. n L t -f- vsi. On aura done 

dx _ ^ dx m 
dt Jmd 1 dwi 3 

mais, eomme 

^3 = ( c — 1 — 7/1) ( /ii — n % ), ( g—i — m )( n 1 ^ n x ), 

dx 

-gj depend de c et de g . Quelle est Ferreur commise si nous sub- 
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stituons xu , Cfc-i el gk~\ a la place de x, c et «*? Si nous rempla- 
90ns x par a?*, nous commetlons une erreur du (A" -f- i) ien,e ordre; 
si nous reinplacons ensuite c par nous coinmettons une nou- 

velle erreur 

, w > dx k 

(«i— n 2 )(c — 


Or c — Ck~\ est du k'* me ordre ; je dis que est du premier 

ordre, car x & — x 0 est du premier ordre, et est nul puisquea? 0 

ne depend que de — <r 2 . L’erreur est done du (k - 4 - i) ieme ordre:. 
et il en est de meme quand on remplace g par gh~\. 

c. Q. F. D. 

Soil done 

(9) A, A*, B, B% G, C* 

ce que deviennent les differences 

_dx dF dF 

cit dX’ cLt dx' dt ' H dX ' ‘ ’ 

quand on y fait cette substitution. 

Les expressions (9) seront des fonctions conrtues qui seront 
du ( k 4- i) u * n,e ordre. 

Posons alors 


x — Xk -+- ox, X = X/e- f- oX, c ss C/t—i -h or. 


et proposons-nous de pousser Papproximation jusqu’au (A*-h i)‘ 4me 
ordre pour jusqu’au A' i, * me ordre pour c. Nous pourrons n^gliger 
les carres de 5 #, oc, . . . , et nos equations prendront la forme 


Vio) 


^ dx ^ dF . „ d\ ^ dF 

0 yr — 0 ~pr = A , 0 _ _i_ 0 — = A*, 

dt dX dt dx 


Les seconds membres sont des fonctions connues; quant aux 
premiers membres, ce sont des expressions iineaires par rapport 
aux inconnues et a leurs derivees. On a par exemple 

„ dF ' d* F\ d'-F* „ v 
° dx dx- °' r dxdX^ u 


ou, dans les deri\ees secondes de F, les inconnues doivent £tre 
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remplacees par leurs valeurs approchees x = x 0l X = X 0 , ... 
(en inline temps que les E et a par zero), ainsi qu’on l'a explique 
au n" 341. D’autre part, 


^ dx d ox N dx 

0 clt dwi Jflt dwi 


Dans le premier termed ng nous pouvons remplacer n$ 
et n % par leurs valeurs approchees 


(c 0 — i — m)( /ii — Tii), (g 0 — I — m)(/ii— n 2 ), 


deduites de Tanalyse des Chapitres XXVII et XXVI. L’erreur 
commise ainsi sur /i/est du premier ordre (et m£me du second), 

et, comme ox est du (A*-f-i) ieme ordre, l’erreur sur sera 

du (/c -f- 2 )'* ,IHJ ordre au moins. 

Dans le second termed 3/2/ ou figurent 

on 3 = (/it — n t ) o/i 4 = (/ii — n t ) 8g, 


nous pouvons remplacer x par x K ; I’erreur ainsi commise sera du 
second ordre, et, comme 3 m est du k' htte ordre, I’erreur sur le 
produit sera du (/r -H 2 )’* me ordre. 

Si nous faisons passer le termed 3 m dans le second membre, 
les equations (io) deviennent 




2 d ox * dF 1 . ^5. dx j 

n ‘d^ - 6 dK =A -l Sn ^’ 

doX ^ dF' . . dXi 

2* ni dW+°iz =k -Z oni d^- 


Les premiers membres restent les m£mes a toutes les approxi- 
mations. 

Dans le calcul des lermes du premier degr£, et en sup- 
posant a = E 3 = o (et aussi3c = o^= oetpar consequent S/z/= o, 
puisque a cette approximation c et g se reduisent k c 0 et g 0 ), les 
seconds membres sont nuls; les Equations (i i) doivent alors se 
reduire k celles que nous avons int4gr<5es aux Chapitres XXVI 
el XXVII, Chapitres dans lesquels nous avons pr£cis<$ment d6ter- 
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mine ces lermes da premier degre; c’est dire que les premiers 
membresdes equations (i i) se reduisent a 


(12) 


I d 


oX — n- 2 oy. 


dt 

( nil -+- in*) o ~ — 2 n\ ox — BY, 
dt ' r* 

— 3 Y -+- 71^ $x ; 

-+• (m i-h m 1 ) o ■— -h 7120X, 

dt 
dlZ 
dt 


~sz, 

-4- ( — /i 2 ) 0 OZ 


[c/. equ. (r), Chap. XXV, n" 323; equ. ( 1 ), Chap. XX\ 1, n° 33i ; 
equ. (6), Chap. XXVII, n M 342]. 

Ohservons d’ailleurs qu’on aurait 

7 . 3 ~ = A ox -+- B 8 y, 

/.o — = B ox ■+■ C 5/, 

^3 

A, B, C ajant merae signification que dans les equations (6) du 
Chapitre XXVII; nous ecrivons aussi pour abreger 

d d 

Jt = 2i n ‘d^’ 

les rti £tant remplaces par leurs valeurs approchees 

n u n i} (co— 1 — m)(n x — ti 2 ), (go — 1 — m) (n x — n 2 ). 


353. Si nous supposons pour un instant que les eonstantes oc 
et Zg (et par consequent les 3 m) sont donn^es, les seconds mem- 
bres sont connus; les premiers membres sont les expressions ( 12 ) 
et nous savons integrer les equations sans second membre; le pro- 
blemepeut done £treconsid£re comme r^solu par 1’application du 
proedd^ classique du n° 350. Nous devons toutefois faire les re- 
marques suivantes : 

i° Les seconds membres se presentent sous la forme de foncrions 
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p^riodiques connues des quatre arguments pp,*; mais dans les pre- 
miers membres figurent non pas les d^rivees 

d d 

dt ~~ 2^ Ul dwi 

mais les derivees 



oii les n\ sont les valeurs approchees 

= /ii, n\ = /i 2 , 

nl = (cQ— i — m)(n l — n 2 \ — — m)(n l — n 2 ). 

Cela ne change rien d’ailleurs au principe du calcul; seulement 
il faut, avant l’integration, remplacer dans les seconds membres 
les Wi non pas par ntt to/, mais bien par n? £-j- to/. 

Supposons done qu’on ait forme les fonctions que nous avons 
appelees a'. au n° 350, 

XJ =a/A+,..; 


ce sont des fonctions p^riodiques donn^es des w; on doit ecrire 
alors non pas 

d'ki vi dlf 

-di = Z nk d^ = aiX+ -' 


mais bien 
Si done 

( 13 ) 


d'ki V' o dli . 

ai A ■+■ . . . = ^ h e ^ 


les k a 4tant entiers, on en deduira non pas 

h c 2/k' a w a 

— — — , 

he^^oL^a. 


mais bien 


*,-2 




Cette analyse suppose que le second membre de ( 1 3 ) ne contien t pas 
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de terme constant. Onverra plus loin, au n° 3 o 6 , comment onpeut 
s’arranger pour qu’il en soit ainsi. 

3 o 4 . On remarquera que les quatre premieres equations (11) 
forment un systeme qui ne depend que de 8#, 8 y, 8X, 3 Y, tandis 
que les deux derni£res forment un systeme qui ne depend que 
de 5 5, 8Z. Ces deux systemes peuvent done 4 tre traites separe- 
ment. 

Nous remarquerons ensuite que nous nous trouvons dans le cas 
ou le procede dun" 351 est applicable, puisque nos Equations sans 
second membre sont les equations aux variations d’equations ca- 
noniques. 

Consid^rons done ies quatre solutions particulieres (7), (12) 
et (12 bis) du Chapilre XXVII; soient 

ox = Ei> oy = * 11 ; ox = ^2? oy — 7j 2 

les deux solutions (7), ou ces deux merries solutions muitipliees 
par un coefficient que nous pouvons choisir arbitrairement ; 

oy = r ti ; Zx = q/ = r j4 

les deux solutions (12) et(i2 bis ), ou ces memes solutions muiti- 
pliees par un coefficient que nous pouvons choisir arbitrairement. 
Soient 

ax = av =: tj? 

les valeurs de 8X et 3 Y qui correspondent k 8«r = £/, 8 y = r u -; on 
aura la relation bilindaire 

( V U ) H- ( vni — T u V.) = const., 

•qui n ? est autre chose que la relation (14) du Chapitre precedent. 

Nous savons que la constante du second membre est nulle, 
:sauf dans le cas ou i = 1, k = 2 et dans celui ou i = 3 , k = 4 - 
Nous pourrons choisir les coefficients arbitraires dontnous venons 
de parler [et par lesquels les solutions (7), (12) et (12 bis) sont 
muitipliees] de telle fagon que dans ces deux cas la constante du 
second membre soit egale a 1. C’est ce qui est le plus commode 
pour Pexposition; mais dans le calcul on pourra faire un autre 
choix; par exemple on pourra avoir avantage pour i = 3 , /r = 4 a 
supposer la constante ^gale a sf— 1. 
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Soient alors 


L, L% M, M*, N, N* 

les seconds inembres des six equations (1 1), de telle facon que 

*- i - 2 »*£- 

Appliquons le proc^de du n" 3 ol ; il viendra 

* Xes 2 M; ’ 

avec les conditions 


U 4 ) 


~ =-5TL + £,L*-^M+r tl M*, 


dAi 


r =-«!*+ 6,L*-r I IM+r 11 M*. 


3 oo. Le meme procede s’applique au second systeme forme des 
deux dernieres equations (it). Soient 

Z = Z = ^6 

les deux solutions (a) de 1 ’equation (1) du Cliapitre XXVI, multi- 
pliees au besoin par un facteur constant arbitrairement choisi. 
Nos equations sans second membre [formees avec oz coinme cette 
Equation (1) avec z] admettront. les deux solutions 

§* = ? Sj sz = a, 

a* = : 6 , ez = 

entre lesquelles nous aurons la relation bilineaire 

const. 

II ne faut naturellement pas confondre les & avec s = 

Nous pouvons supposer que la constante du second membre est 
£gale a 1. Nous aurons alors 
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avec 


1 15) 


1 ^X5 
dt 
d \ 6 
dt 


= tfN-C.N*, 
= — «N-hC«N\ 


3o6. II faut, comme nous l’avons vu a la fin du n° 353, que 
les seconds membres des Equations (i4) el ( 1 5), qui sont des fone- 
tions periodiques des cv, ne contiennent pas de terme constant. 

£?X 

Examinons-les successivement el consid^rons d’abord Je dis 

que cette expression est une fonction impaire des w et ne contient 
pas de terme constant. II est aise, en effet, de constater que 

Si» r iu M 

sont des fonctions paires, tandis que 

61, ^ L, M* 

sont des fonctions impaires, ce qui demontre la proposition 
enoncee; il ne peut done pas s’introduire dans X 2 de terme secu- 
laire. 

Passons a ; nous savons que £, est egal, a un facteur constant 
pres, a et que est ^gal, a un facteur constant pres, a 

t da ? 0 dar Q . 

m dz dm 3 


nous pouvons done choisir les facteurs constants de telle sorte que 
Ton ait 

5s ” ^1 +($s)i £2 = (?2 )> 

i)s = <i)i-t- (*)*)» 11?== <1)1 H-(* 1?), 


rp, £*, 71% (£2), (via)? (i*)? (*n*) etant des fonctions periodiques 
de Posons alors 


il viendra 


Xi = — fX 2 -f- fxi ; 


d’oii 


. <^Xa ^ dpi 

dt ~ dt dt 9 
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Le second membre est encore une fonction periodique des 
paire cette fois; mais dans ce second membre figure A 2 qui n’a 
ete determine que par une integration, e’est-a-dire a une constante 
pres. Nous pouvons disposer de cette constante de telle facon que 
le terme constant du second membre disparaisse. Dans ces con- 
ditions p, ne conti endra pas de terme seculaire. 

D’ailleurs les deux premiers termes de ox 

XiJi-4-A2$2 

se reduisent a 

f*i$i -*-**(&). 

de sorte qu’on a 

X 2 (b) -+- X 3 ?3-4- X*£ 4 . 


357. Examinons maintenant et J e dis que nous pouvons 

choisir oc de facon a faire disparaitre a la fois le terme constant 
dans ces deux expressions. Nous avons 




Comme o/Ji eto/i 2 sontnuis et que 


dx i 
dw w 


o, puisque x t ne depend 


que de t — — w 2 et de tv ;} , le dernier terme du second membre 

se reduit a 


^ dx 1 • * , , dx 1 

— on 3 — T — = — oc ( ni — n > ) 

dw % ' dw z 


D’ailleurs on a 


*i=E t«3+ 50- 


Nous voyons done que 

peut se diviser en deux parties et qu’on peut ^crire 

— <p 3 — — n 2 ), 

ou est ce que devient l’expression de ~~ quand on y remplace 


L f h*, M, M 
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A, A*, B, B*. 

et i 8 ce que devient cette expression quand on y remplace ces 
quantites par 

da ?i dX i dy\ d Y, 
dZ? dw 3 ’ dw z ' ~d\v z 


On aura de meme 

= ?4— '{'V Sc(rt, — /ij), 

« 4 et 'i 4 elant formes avec ^ comme » s et avec Je vais dis- 

poser de Sc de facon a annuler le terme constant de je dis que 
la \aleur de Sc qui annule ce terme est reelle et qu’elle annule en 
meme temps le terme constant de °-~- 
Nous avons suppose 

•+■ «is = A ^ , 

(.-•VAjetW 


les coefficients 6 a et c* sont reels, et h cst un coefficient constant 
dont nous avons dispose de fagon a reduire une certaine constante 
a i. 

En effet, quand je change les <v en — (*>, x qui est une fonction 

• • dx dx 

paire des w ne change pas; au contraire, -jp X, y , 

changent de signe. 


dx i 
dw z ' 


B% 


dw 3 


changent de signe; 


dX , 

3 »,’ B ’ 




ne changent pas. 
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^ et se changent en ^ et tandis que ^ et ^ se changent 


uau^^m, v^xx c. c y“ j kauuis *-|ixc Cl j SC 

Done ~ se change en — y~ en 


xi 2* 
h ' h 


en — Ji’ ~fi en — X* ^ onc terme constant de par exemple, 
est egal a celui de — ^ • 

Si maintenant nous changeons i en — L ~ -?•* Yj A se per- 

° ’ h h ft h 1 

mutent avec jL, les quantites A, A*, qui sont 

reelles, ne changent pas. 

Done jp -jr* -jr se changent encore en — — -jr * — > 

— Done le terme constant de est imaginaire conjugue de 
celui de — 

h 

Si le terme constant de ? 3 est d’une part egal a celui de — 

d’autre part imaginaire conjugue de celui de — c 7 est que ces 

deux termes sont egaux et reels. Et il en est de m£me en ce qui 
concerne et i 4 . 

Soient done 


a, (3, — a, — 


les termes constants de 


?* 4*3 ?4 

h 7 X* X 5 


ils seront r^els, et il suffira de prendre 


Hi) 

pour annuler a la fois le terme constant de et celui de Nous 
n’aurons done de terme seculaire ni dans X 3 ni dans ), 4 . 

C. Q. F. D. 


On demontrerait de la m^me maniere qu’on peut choisir S^* de 
fa$on qu’il n’y ait de terme seculaire ni dans A$ ni dans X*. En 



CHAPITRH XXVIII. 


92 

eflet, le second membre de la cinquieme equation (1 1) s^crit 


n = g-2 


on i 


dwi 3 


on a ici 


^1 ■> dzi * , \ dzi 

2* oni asr, = ora > ^ m ^ ■ 


Le reste du raisonnement s’acheve comme pour X 3 et X 4 , 
og jouant le role de Sc, IN et N* celui de L et L*, £5 et C G celui 
de i 3 et i 4 , etc. 


358. L’integration, comme nous l’avons vu a la fin du n° 353, 
introdiiit le petit diviseur 


ou les k a sont des entiers; comme 


/i§ = Co(rti— no) 

et 

nl = g 0 (ni— n t ) 

sont d6veloppables suivant les puissances de m 2 , ce petit diviseur 
sera lui-meme developpable suivant les puissances de m-. 

Les premiers termes du developpement sont 

/i 2 = n% = m(n t — n z ), rii = n\ = (1 ■+■ m) {n x — n 2 ), 

«3= (*1 — .), 

nj = (m-iit) 

puisque 

3 3 

c = h-/?2 + 7 /?z- + ..., g = i ni — 7 . . , 

4 4 

d’ou 

a/l ^ ~ ^‘l -r- (^1 -f* -r - ( X' 3 — £4) ttt 2 -+- .... 


11 sera done divisible par m si k i — o; il sera divisible par m- 
si k\ = = o, c J est-a-dire si le terme correspondant ne depend que 

des longitudes du perigee et du nceud; dans ces deux cas il sera ce 
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que j’appellerai un petit diviseur analydque . Enfin il sera divi- 
sible par m* si k { = k 2 = o, /r 3 = k :t , c’est-a-dire s’il depend seu- 
lement de la somme des longitudes du perigee et. du noeud. Ce 
sera alors un tres petit diviseur ancdytique. 

Mais il arrive ici que les termes en m z ont de tres grands coef- 
ficients, de sorte que c Q — i — m 7 au lieu d’etre a peu pres egal en 

3 

valeur absol ue a - m 2 , et par consequent a g 0 — i — /??, est a peu 

pres deux fois plus grand. Il en resulte que les tres petits diviseurs 
analytiques, quoique di visibles par m ; \ sont numeriquement de 
l’ordre de m‘ 2 . Si, au contraire, on a /r 4 = = o, 2 /r 3 = A* s , le 

diviseur, quoique non divisible par sera numeriquement de 
I’ordre de m z . Ce sera un tres petit diviseur numerique (ine- 
galite de Laplace, cf. Tisserajvb, t. Ill, p. 1 58 ). 

Dans le calcul numerique des coefficients, ce sont les tres petits 
diviseurs numeriques qui importent Au contraire, si Ton se pro- 
pose, comrae le faisait Delaunay, de developper ces coefficients 
suivant les puissances de /?*, il faut s’inquieter des tres petits divi- 
seurs analjtiques. 

Les tres petits diviseurs analjtiques se presentent pour la pre- 
miere fois dans les termes en EjEijE* ou en EjE 2 E.', et les tres 
petits diviseurs numeriques dans les termes en Ei;E : *a, EiEgEjJ, 
EiE 2 E?,a ou E^EijEij. 

339. Tel est le principe de la methode de Brown. 

Je n’insisterai pas sur les perfectionnements de detail qu’il y a 
apportes et dont les principaux sont les suivants : 

Au lieu de quatre equations lin^aires du premier ordre 
a second membre, il emploie deux equations lin^aires du 
deuxieme ordre a second membre (obtenues par l^liminalion 

de X et Y); il en resulte que les expressions des — sont un peu 

inodifiees et se presentent sous la forme d’une somme de deux 
termes seulement et non de quatre. 

a 0 Au lieu de §# et o ly, il prend comme inconnues les quantiles 
imaginaires 

(lu — ox -h i cjy, os = £x — iSy. 

3° Au lieu d’employer des lignes trigonomelriques des multiples 
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des w, il simplifie les notations en introduisant la variable 

£ = e 1 ’. 

4° Pour les approximations d’ordre eleve, il a recours a Fartifice 
par lequel Hill etait passe des equations ( 2 ) aux equations (5) du 
Chapitre XXV. Les calculs de substitution s’en trouvent un peu 
simplifies. 

Je me bornerai a dire que la methode est applicable aussi bien 
au calcul des termes du premier ordre en a et E 3 qu’a celui des 
termes d’ordre superieur. 

Pour plus de details, je ren\errai a son Ouvrage original {Me- 
moirs of the Royal Astronomical Society } t. LI1I, L1V et LV1I).. 

On pourrait se demander si, a cause des petits diviseurs divisibles 
par m 2 ou on n’arrivera pas dans la suite des calculs a des 
termes contenant en facteur une puissance negative de m. On peut 
demontrer que cela ne peut arriver que quand interviendront les 
tres petits diviseurs anaiytiques . 11 enresulte, d’apres le numero 
precedent, que cela ne peut arriver que pour des termes d’ordre 
tres eleve; cela ne peut arriver si l’on suppose E 2 =o, on 
bien E 3 2 = o, puisque dans ce cas il ne peut j avoir de tres petits 
di\iseurs analvtiques, mais tout au plus de petits diviseurs anaiy- 
tiques divisibles seulement par/n-. Pour la demonstration, je ren- 
verrai au Bulletin astronomique, t. XXV, p. 32 1 . 
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360. La seconde methode que nous allons exposer presente sur- 
tout des a vantages quand on veut obtenir non seulement la valeur 
numerique des coefficients, mais leur developpement analytique 
en fontion de m, comme le faisait Delaunay. 

Nous aurons avantage a employer, an lieu des arguments tv/, les 
suivants : 


T = W»i C*> 2 , Ti = iVi -f- T* = W x -+- ^3 — W 2 . 

De cette facon, dans le coefficient d’un terme dependant du sinus 
ou du cosinus de 

pZ -4- pi Ti -+- p 2 Z» -I- 

f’exposant de E/ sera au moins egal a | pi [. 

Nous partirons de la formule (25) du n° 320 : 

Cette formule peut recevoir utilement di verses modifications : 
d’abord nous pouvons passer des variables aux variables t; nous 
pouvons ensuite remarquer que ( I>, se reduit a une constante ( in- 
tegrate de Jacobi) lorsque E : j= o. Gela nous permet d’ecrire 

<!>' 

= K-+-E 3 <l\ 

II 2 

K etant une constante et E ;i 4> etant divisible par E s . On a en eflel 

d>' 

Si =6 h-E 3 9, 

6 et Q etant des fonctions de y, X, Y, Z, de t 3 =w> 3 et 
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des constantes E 3 et a. D’ailleurs 'i est ind^pendant de E 3 et t 3 . 
Soit ensuite x* ce que decent le developpement de x quand on y 
fait E 3 = o. 

Soit ce que devient i quand on y remplace x, . . . par x * 7 . . . ; 
alors A* sera une constante, A — A* sera divisible par E 3? etnous 
pourrons poser 

= K, -t- E 3 e = E 3 4>. 

On a done 

Af dwi — ■ - = B dz -f- Bi dz\ — h B 2 dz% -+- B,j dz 3 — E 3 dz$, 

d’ou 

dX — dQ” = ^Bdz — 

Les B sont des constantes, de merae que les A) et K; je veux dire 
par la qu’ils dependent seulement de 

Ej. E 2 , E.j, oc, m 

et sonl independants des t. 

Posons 

S = Q" — x Q X — y 0 Y — z\ Z, 

Xo et jK 0 etant les termes de degr^ z 6 ro calcules au Ghapitre XXV ; 
z { represente Fensemble des termes determines au Chapitre XXVI 
(comme z ne contient pas de terme de degre zero, on aura z 0 = o) ; 
il \iendra 


0) 


dS = 2 ^ — ^dX—^Xdxo 

-+" ( -3 — ) dZ — Z dz 1 — ^ B dz -+- E 3 ^ dz% . 


Dans cette formule (1), on doit regarder E,, E.>, m et les t 
comme des variables; au contraire, E 3 et a sont des constantes 
donndes, de sorte qu’on aura 

c?E 3 = dv. = o. 


J’ajoute que dans cette formule ^(37 — x 0 ) dX. et^^^o re- 
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presentent simplement 

(or — Xq) dX -+- (y — y&) Xdx^^dy^ 

361. Ccla pose, supposons qu’on ait determine les termes 
d’ordre k et d’ordre inferieur de x et de y, de X et Y, de S et B, 
les termes d’ordre k — i dc z et Z, et qu’on ait par consequent 

( 2 ) x = x/ l ., y = y it & — za-i , S = S*, 

Je pose 

(3) x = Xk- J r OX, y =y A -h by , 

et je me propose de calculer ox, 3)', oS jusqu’aux termes 
du (k i)' ime ordre inclusivement, 5s jusqu’aux termes du 
k iimc ordre. 

Substituons d’abord dans ( 1 ) k la place de toutes nos variables 
leurs valeurs approch^es ( 2 ), la difference des deux membres sera 
d u (/ f iy«erae or( j re • nous pourrons la mettre sous la forme 

u et v etant des fonctions connues. 

Substituons maintenant a la place de ces variables leurs valeurs 
approchees (3) ennegligeant les puissances supSrieures de ox, . . 
ce qui cst permis puisque nous n^gligeons les termes d’ordre k -4- 2 ; 
il viendra 

d8S =^8xdX -hj?(x — x 0 )d8X. — dx Q -h oz dZ 

-^(z — z^doZ — SZ dzi — ^SB dz -f- E 3 8<* dv. 

Le second membre est susceptible des simplifications suivantes : 
consid^rons-en d’abord la premiere ligne ; comme ox et 3X sont 
d’ordre k -+- 1 , nous pouvons remplacer x et X par et X 0 . De 
m^me dans la deuxieme ligne, comme 3:s et 5Z sont d’ordre k , nous 
pouvons n^gliger dans leur coefficient ■3o qui est de deuxieme 
ordre et remplacer z et Z par s { et Zj . 

Enfin 8<I> est d'ordre A* -H 1 , car 



P. -II (a). 
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Or ox et 3X sont d’ordre k + 1 , 02 et 3Z sont d’ordre Ar, 
mais et sont divisibles par E 2 et par consequent du premier 

ordre. Comme d’ailleurs E 3 est du premier ordre, nous pourrons 
negliger E 3 3<I> et il reslera 

d SS =^ 5 # dX 0 - ^ SX dx 0 

-4- 0-3 — SZ eft, — SB d~ -f- u dv. 



362. Prenons xnaintenant la formule (26) du n° 320, 


dor 

dt 


►s+2‘ 


eft 


— H, 


ou ^ = F' — n 2 v r et oil H est une constante choisie de telle fagon 
que Q rr soit periodique. Nous deduirons 


( 6 ) 


dS 

dt 


= + *0) 

dZ 


dX 

dt 


H-(*~ 


flfZ ^ e£r 0 

} eft 2d dt 


— z 


"ST 


- H. 


Substituons dans (6) les valeurs approchees (2) a la place des 
variables et soit G la difference des deux membres; G sera une 
fonction connue du (A* -4- i) i6,ne ordre. Substituons maintenant les 
valeurs ap proch^es (3) et negligeons tout ce qui est du (A* + 2 ) ihme 
ordre ; il viendra 


; dS 

dt 


>«f-2 




SX 


eft 




D’autre part, 




'^SX- 

I Ct A 


d$>\ y 

!u QZ ' 


d*\ 

dZ 


SZ. 


Mais, en vertu des Equations du mouvement, on a 


d’ou 


dx 

dt 

idX 


dF f 
dX : 


d®[ 
dX j 

dx , 


^\dX* Vi dx e/Z . dz 
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et par consequent 


3^ = G- 
at 




% dZ 

(S - Zl) t-, 


(dz dz, 
1 dt 


dt 


)* 


Z. 


Dans la premiere ligne on peut remplacer j; par et, dans la 
seconde, s par pour les memes raisons qu’au numero prece- 
dent, ce qui nous permet d’ecrire 

( 7 ) zf t =G-m. 


Qu’est-ce maintenant que 8 On a 


d’ou 


dS v-i dS 
dt ni du>i 


^ dS 
o 


dS ^ do S N 

dt=2* ni d^+Z ani 


dS 

1 dwi 


Dans le premier terme du second membre nous pouvons rem- 
placer m par de sorte qu’il se r^duira (en reprenant les nota- 
tions du Chapitre precedent) a 

"V o __ d oS 

2mi ni dw t ~~ dt 

Dans le second terme figurent deux constantes indeterminees, 
8/13=30(71! — /i 2 ), on^=:Bg(n 1 — /z 2 ); 


la premiere est d’ordre A**, car nous supposerons que c a ete deter- 
mine jusqu’aux termes d’ordre k — i inclusivement ; la seconde 
sera d’ordre k — i , car nous supposerons que g a ete determine 
jusqu’aux termes d’ordre k — 2 . Posons 

S = Sq-t- S t -+- Sj , 

So, S 4 , S 2 , ... representant respectivement les termes d’ordre o, 
i, 2 , Nous pourrons alors, dans le coefficient de S/a 3 , rem- 

placer S par Sq — {— S* ou par S f , puisque S 0 ne depend pas de w 9 , 
k et, dans le coefficient de 8« 4 , remplacer Spar S 0 -f- S 4 4- S 2 ou 
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par S 2 , puisque S 0 et S< ne dependent pas de II vient done 
(8) - = G OH — 0,* 3 — - Ort; 


ilw 4 


363. Reportons-nous aux notations da n" 318; nous avons 
trouve dans ce nuinero les form ales suivantes : 

<t>, == FI —V W/t, W i =A,= A' ( (i = i, 3,4), 

/ijWj -j- 4*1 — /i> A 2 -f~ K. = n% A 2 , 

= 2 W = — 2 » rfW, 

Ja leltre K. aj'ant le mdme sens qu’au n° 318, et nous en tirerons 
(en nous rappelant que dn 2 = o) 

H 2 71 ' = °* 

Mais il vaudra mieux revenir aux arguments 7 que nous avons 
introduits au debut de ce Chapitre ; soit done 

dz = v dt , = v* c/t, 

de telle sorle que 

V ~ Tl j — rt 2 , Vi = /2j -h /l 3 = V£, v 2 = /2i-h /l* = V^, V 3 =/i 2 . 

Nous aurons 

Kvj-i-^Bv =2 A ' rt ’ Krfv,+2 B <A/ =^AV/i, 

la lettre K ayant le m^me sens qu'au numero precedent, et par 
consequent 

(9) Hss^JBv-hKvj, dR ^viB+v 3 c/K = o. 

Toutes ces quantiles H, B, v ... dependent des constantes /», 
E et a, et ne dependent pas des arguments 7. Supposons que dans 
les equations ( 9 ) on substitue d’abord les premieres valeurs ap- 
prochees («) de nos inconnues, et soient 

< 10 ; p. 2 Qe/? ’ 
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les differences des deux membres; P, Q, q seront des fonctions 
connues et les expressions (10) seront d’ailleurs du (/" + i) ,ern< * 
ordre. 

Cela pose, substituons dans les equations (9) les valeurs plus 
approchees (3) et negligeons les termes d'ordre k -4- 2; il viendra 

| m =2 B 3v + 8B 4 - P ■+■ v 8 oK, 

^i) 

I doli d ov dv -hJ^Qdq. 


Remarquons que 3H, <^3H et SB sont d’ordre k -r 1 ; que 
ov = ov 3 =o; que 3v< est d’ordre k et Sv 2 d’ordre k — 1; que 
d 3 v/ est du meme ordre que 3 v/ ; que Bi est divisible par E, et B 2 
par E“ et sont par consequent du second ordre, ce qui permet 
de negliger, par exemple, B,3v { ; alors nous ecrirons 


( 12} 


| oil = B 2 ov a SB -i- P h- v 3 3 K, 

I d oil = B* d 3 v* +2 SB dv -4-J? Q dq. 


364. Tous les termes de nos d^veloppements contiennent en 
lacteur un certain monome 

l± = ofl° E'f 1 E^ 8 ’ , 

que Brown appelle leur carcicteristique ; la somme des exposants 

q$- f- qz 

est le degre du terme. Dans les calculs precedents, nous pouvons 
supposer que 3S, par exemple, ou 5x, au lieu de representer tous 
les termes de degre k + 1, par exemple, represente seulement 
l’ensemble de tous les termes ayant une caract^ristique donn^ep. 
Mais, comme le degre d’approximation n’est pas le rn^me, par 
exemple, pour Sz et ex, il est n^cessaire que je precise. Je con- 
viendrai done que 

$ar, oy, 8X, oY, $S, 8B, 3K, 8H 

represen tent 1’ensemble des termes de caracteristique p; que 3s, 
dZ represented 1’ensemble des termes de caracteristique jjr; que 
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oc comprend les termes de caracteristique et og ceux de carac- 
teristique • 

En ce qui concerne les valeurs approchees ( 2 ), je supposerai que 
xk, yio S j { , ... 


comprennent toils les termes dont la caracteristique est un diviseur 
de p. ( 4 x lui-m£me etant exclu) et que de meme 


Zfc—l j C^u gJc-1 

comprennent les termes ayant respectivement pour caracteristique 
urx diviseur de 

[X JJL JJL 

K’ Wy’ K| * 


365. Cela pose, nous devons distinguer trois cas : 


i° u n'est pas divisible ni par E, ni par E 2 . 

Dans ce cas ni*r,, ni (on ce quirevientau meme ni w z ni pp> 4 ) 
ne figurent dans nos developpements. La formule (8) se reduit 
done a 


6t 


= G-SH; 


G est connu ; on disposera de oH de telle fa$on que le terme con- 
stant du second membre soit nul, et l’on aura SS par une simple 
quadrature. La fonction SS est ainsi determinee a une constante 
pres, mais cette constante doit etre nulle puisque SS doit £tre une 
fonction impaire. 

Passons maintenant a la formule (5), et remarquons que Z* , . . . 
sont nuls, et de plus que nous n’avons que trois variables par 
rapport auxquelles nous puissions differentier et qui sont m, t 
et t 3 ; il vient done 


U3) 






dX a 


5X 


dm 

dx 0 
dz 


dv 

' dm } 


SB- 


2 dv 


et, puisque x„ et X 0 ne dependent pas de T a , 


(t i 


dSS 


V 1 „ dv 
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De liquation (r4) on deduit 


-St-*]- 


en designant par lc ternie constant de u ; et cela deter- 

mine oB, v 

On remarquera d’autre part que dans la seconde equation ( 12 ) 
on a B 2 = o, puisque B 2 doit etre divisible par E; et que, u n’etant 
pas divisible par E 2 , nous negligeons E 2 ; il reste done 

del H = 

ou, puisque dv 3 = o, que dv { et ofy 2 n’interviennent pas, 
dm dm dm 

ce qui determine SB, puisque 3H, Q et q sont connus et que 


Nous avons enfin 

j>y ? ^ •> s.x’- > dy -> 

o\ — 0 — “ /io o v ) Ol 0 — r~ — {— /Zj ox 

W f fli 

Comme ici nos d^veloppements ne contiennent ni w$ ni «\ 5 , et que 
Ton a d’ailleurs 


nous poumons eenre 


h —— - _ 37 % dV _ d ojk # 

dt ~~ dt ’ ° dt dt * 


mais, comme nous retrouverons les memes equations un peu plus 
loin, j’aime mieux trailer tout de suite la question d’une facon un 
peu plus gendrale. Reprenons done 1’dquation 


dx 

— -X-n t jr = o, 


et soil A ce que devient le premier membre quand on j substitue 
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les valours approchees (a). Nous aurons, en prenant les valeurs 
plus approchees (3), 


D'ailleurs 


^ dx 
0 dt 


~ 8X — rii oy = A. 


c&c XT' cfcc „ cfo? XT' c? Zx rs cfo 

dt ^d nt dwt ° dt Mmi 11 dwi jZd. 1 1 dwt 


Nous pouvons remplacer dans le premier terme n/ par n° £ puisque 
ox est d’ordre A* 4- i, de sorte que ce terme se reduit a -^r> dans 
le second on a 



dx 

dwt 


= V oc 


dx 

d =7 


“T* v o g 


dx 
dx 2 


Or, oc et o** etant respectivement d’ordre k el /c — r, nous pou- 
vons, dans le coefficient de oc, remplacer x par x { et, dans celui 
de og , remplacer x par x» (oil x 0 , x i} x» representent les trois 
premieres approximations de x). Nous pourrons done ecrire les 
Equations suivantes : 


( i5) 


o Zx 

“5T 

dZy 

dt 

,dX o 
dm 


— ax- 


* ^ ciur I J uu , 2 

n s o r = A -voc^- -vdsr-gr> 


3Y -h n 2 ox = A' — v 8c — v Zg-^~> 


dm dm 


12 ** 

2 ^- 2 * x ^ 


dh 

’ dm 


3Z 


i ^ dv 
1 ~~2d U d^i 


dm 


d 5S 
dx 


dx dx Jrni dx 


Les termes en 8s et 3Z sont nuls dans le cas qui nous occupe, de 
meme que A, A' et les termes en oc, Sxr ; mais je prefere completer 
tout de suite les equations (i5), afin de pouvoir encore m’en 
servir dans les deux numeros suivants. 

3S et SB sont connus; les arguments et t 2 n’intervenant pas, 
nous devons regarder 8c et 3 g comme nuls; les seconds membres 
sont done connus; nous avons done a integrer un systeme d’equa- 
lions lineaires a second membre, Les equations lineaires sans 
second membre ne sont autre chose que celles que nous avons 
integr^es aun° 347; seulementnotre systeme d’^quations dififeren- 
tieiles est du deuxi£me ordre au lieu du quatrieme. 
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Nous n’aurons done qu’a appliquer les methodes des n us 349 
et suivants ; il n’y aurait de difficulty que si le determinant 

dx o dy n dx 0 dy 0 
dm d-z drz dm 

pouvait s’annuler, ce qui n'a pas lieu. 

11 ne s'introduira pas de terme seculaire; cela ne serait possible 
que si l’argument de l’un des termes du second membre etait le 
rneme que celui d ! un des termes de ee que nous appelions dans le 
Chapitre precedent ou e’est-a-dire 

~ i -+- ( 'X k -+- 1 ) ~ . 

Or cela est impossible, puisque nos seconds membres sont inde- 
pendants de t, et 


366. Passons maintenant au second cas : 


2 ° u est di\ islble par E 2 , mais pas par E, . 

Alors nos fonctions dependront de (e’est-a-dire de mais 
pas de cr 3 (e’est-a-dire de T|), et l’equation (8) s’ecrira 


(8 bis ) 



8H — 


dS* 


Comme le second membre ne doit pas avoir de terme constant, et 
que — — ■ n’en a pas, SH ne sera autre chose que le terme constant 
de G; SH etant ainsi connu, les equations ( 12 ) donnent 


d $H _ £Sv 2 
dE> ~ 4 dE 2 




dq 


Or v ne depend que de m, v* n’intervient pas et v 3 est constant; 
on a done 

A-a— 

d’ou 


rf8H 
dh* 2 


• B* 


d 8v 2 

"517 


?R dWt 

oB! 5e; 


-2«afc 


Nous prendrons oB 2 arbitrairement (en le prenant toutefois nul si 
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les exposants q de la caracteristique ne sont pas toils pairs). Tout 
sera connu, nous en tirerons done cette quantite el par 

consequent 

Sv 2 = 


d 8v 2 


<72 — 2 r/E 2 


puisque Sv 2 est homogene d’ordre q « — 2 en E 2 . 

Connaissant ovo = o/i 4 , nous connaitrons le second membre de 
Fequation (8 bis) et par consequent oS a une constante pr&s qui 
est nulle, puisque 3 S est une fonction impaire. 

Cela pose, venons aux equations ( 5 ); elles nous donnent 


(16) 


/ dl S 
\ dE, 
j d3S 
[ dz 2 


dz. 



-2 

8B s +2 


dv 

“ 3 e 7 

dv 

7/ * 


w, 5,, Z| sont des fonctions connues, SB 2 a ete choisi arbitraire- 
ment, on vient de determiner SS; nous pourrons done determiner 
sans integration les deux inconnues restantes oz et 8 Z a Faide 
des equations du premier degre (16). Le determinant de ces equa~ 
tions, 

d r L\ dz>\ dz\ dLi 
dE* dz 2 dEf> ciz-2 


ne peul s’annuler, car il se r^duit a une constante; on n’aura done 
pas, pour resoudre les equations (16), a effectuer de division. 

Les equations ( 5 ) nous donnent ensuite 



ce qui montre que oB 3 est egal au terme constant de ^ ll ~dT* Les 
equations (12) donnent 

am dm dm dm dm 

Tout etant connu excepte SB, cela determine SB. 

Venons enfin aux equations (10). Dans ces equations, t 4 n’inter- 

venant pas, tout se passe comme si Sc etait nul; Sg= ^ est connu; 
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on vient done de determiner oS et3B; tout est done connu, sauf 

Sr, S>% SX, SY. 

Ces quantites se determineront done facilement par Integration 
des equations (id); on demontrerait comme au numero precedent 
qu’il ne peut pas s’introduire de termes seculaires. 


367. Passons au troisieme cas : 

3° jx est divisible par E, et par E^. On peut alors eviter une inte- 
gration. 

Les equations (3) nous donnent 

d SS __ vi dv 
~2d u dK [ ’ 


car# 0 , X 0 , z i9 Z, ne dependent pas de E, ; si alors SS et v sont 
homogenes de degres q K et k en E,, on en tire 


ce qui determine SS. 
On a ensuite 


?.ss=2 


huv, 


dZ S 
a?E 2 
dZ S 
dz 2 


— 06 — 77 T - 

fl^E 2 g?E 2 


-2- 


dv 

d\ii* 


* dZ t dz\ ^ dv 

— oz SZ 8B 2 -+* 

dZr, dz •-> uTj 


>7- 




ce qui determine 82 et 3Z, la constante 8B 2 pouvant £tre choisie 
arbitrairement. Puis 


dob ^ dv dob V 1 

1=7 = - oB 3 -’-2“^' =- oBi+ 21“^r' 

ce qui montre que 8B 3 et SB, sont egaux aux termes constants de 

2 dv V71 dv 
“d 7s’ 2d U dT l ' 


II reste a determiner SB et 8^, d J ou depend 8v 2 ; pour cela nous 
nous servirons des equations (ta), qui nous donnent 



o8 


CIIAPlTKK XXIX. 


Dans le coefficient de SB, nous pouvons remplacer les v par leurs 
valeurs approchees, v = v, = c 0 v, . dans ces condi- 

tions, les v ne dependent pas de E, et il reste 


cU H 
dH t 


m _ d 3v, , aq 

£ t 2 dEi ‘ Zd ^ d\ii 


dq 


et Ton aurai t de meme 

d 3H _ d Sv 2 

~ H 


d\ 


dE 2 


2«afc' 


et Fon en deduit 


#1 3H =s B 2 ^i 3v 2 


q-i oH = B 2 ( q * — 2 ) Sv 2 4- 


q etant suppose homogene de degres k et Id tant en E, qu’en Eo. 
De ces deux equations on tirerait 8H et 3 v 2 (d’ou og ). 

Le procede deviendrait illusoire pour q 2 = o, mais dans ce cas 
Fargument tt 2 et par consequent Sv 2 ndnterviennent pas. 

On trou ve ensuite 


d oH ^ g d Sv 2 
dm ~~ 2 </m 







2« 


6 ?^ 


doii l’on tire SB. 

On determinera enfin So?, Sy, 3 X, 3 Y par le moyen des equa- 
tions (10); dans les seconds membres tout est connu, a Fexception 
de la constante 8c. On disposera de celte constante de fagon a faire 
disparaitre les termes seculaires, qui cette fois ne sont pas nuls 
d’eux-m£mes. La determination de toutes nos inconnues est done 
achevee. 


B68. Cette methode a ete exposee, mais sous une forme et avec 
des notations differentes, dans le Tome XVII du Bulletin astrono- 
mique, p. 87 et 167. Quand on veut Fexpression analytique des 
coefficients, elle presente Favantage de rendre plus rapide le travail 
de substitution (puisqu’on n 5 a qu’a faire les substitutions dans <£' , 
au lieu de les faire dans les trois d^rivees de cette fonction) et 
d’amener a Fint^gration d’un systeme du deuxieme ordre au lieu 
du quatrieme. EUe est susceptible de plusieurs variantes * 
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1 “ Nous pouvons employer an artifice analogue a celui des 
347 et 347. en envisageant un probleme plus general que le 
probleme propose. Nous avons 


SX' 2 

'X 


7>2| ■ 


n. 2 { Xr — Ya*) 



*2 


/l l 0 , 


le dernier lerme — representant fensemble des termes conte 
nant en facleur a ou E 3 . Prenons la fonnule plus generale 


SX 2 


m 1 -+- mi 


-f- n 2 1 Xy — Yx) 


x 2 -hy 2 


■ Z 2 X- 

— 3 A 2 — 


-r 2 


-A 2 e, 


qui se reduit a la premiere pour h = /zo. Xu n" 347, nous a\on> 
pose 

no = /? v, /i = /;?v; 


nous poserons cette fois, ce qui revient au rneine, 

72-2 = 771V, Jl = S 7?l v . 

Nous appliquerons ensuite la methode en developpant, non plus 
seulement suivant les puissances de a et des E, mais suivant celles 
de j3, de a et des E, de telle facon que .r 0 , par exemple, represente 
fensemble des termes independants de a, des E eC cle k 3. Le 
nombre des termes a calculer se trouve un peu augmente ; en 
revanche, # 0 ,jr o , X 0 , Y 0 , -i, Z, se reduisent a un seul terme en 
cos t, sin t, cos t-j ou sin t 2 . 

2 ° On pourrait, au contraire, achever d’abord le developpement 
par rapport a a en regardant les E comme nuls, puis developper 
ensuite par rapport aux E, en regardant x 0 par exemple comme 
fensemble des termes de degr£ zero par rapport aux E seulement, 
mais de degre quelconque par rapport a a. Ou bien developper 
d’abord par rapport a E* et E 2 , et ensuite par rapport a a et E ;l . 
Cela entraine dans la methode quelques petites modifications sur 
lesquelles nous n’insisterons pas. 

3° Au lieu de prendre les coordonnees rectangulaires y y s et 
leurs variables conjugates X, Y, Z 7 on peut prendre les coordon- 
nees polaires et leurs variables conjugu^es. La methode fondee 
uniquement sur les proprietes des equations canoniques restera 
applicable, sauf quelques modifications de detail. 
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Th6qremes d’ Ad axis. 


369. Reprenons l’equation du n° 362, 

dt ~ 9i ■ Zi dt ’ 

ou cette fois^#X signifie xX-hyY + s Z. Posons 

ccX 


V= Q"- 


il viendra 


2 X A 
~ 9 


d\ i ^ dX i dx 

Mais 

dx_dF dX _ d¥' _ d*>\ 

dt <a£X ~ dX 7 dt ~~~ dx ~~ dx 

Done 

< EL = v, iy~g*l IVxgj H. 

aff 1 2 dx 2 ^ d'L 

Soit 

o* repr^sentant Fensemble des termes homogenes de degre k en 
ic, jK? X, Y, Z; on aura, en vertu du th^oreme des fonctions 
homogenes, 

2*2*2*£-2*- 

d’ou 

S-2(— t)**-"- 

V ^tant une fonction periodique, H ne sera autre chose que le 
terme constant de <p*. 

Mais, si Fon neglige la parallaxe, on a 


_ , 2X 2 77ii-i-m 7 Q ^ _ 

= — f-n 2 (\7* — Yor) — /if PjAG 2 




{c/. n° 315). 

Tous les termes sont homogenes de degr£ 2 , sauf le terme en 
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qui est de degre — i . Done 

2(-!)»-. !a r 

Done j si la par allaxe est nulle, H nlest autre chose que le 
terme constant de y au facteur constant pres 

3 , 

- (/«i h- m : ). 

370 . Cela pos£, prenons 1’equation (9) du n° 363 , 
dH = 2 B d ' ) - 

Comme v et v 3 ne dependent ni de E, ni de E 2 , elle nous donne 

I'd ± 

M dEt ’ 

dv 


d'J dv 

dEt* 2 dEl 


H = H 0 +H j +H 4 + ..., 

H* representant l’ensemble des lermes de degre k par rapport a E 4 
et Eo, et de degre quelconque en a et E 3 . Le theoreme des fonc- 
tions homogenes nous donne 

„ dH _ dH „ 

El dEl * E * dEl = 2Hl+ 4 H * + - • " 

Done 2H2 repr^sente l’ensemble des termes de degrd 2 dans le 
second membre de (18). Or B { et JB 2 sont di visibles respectivement 

par EJ et E* ; d’autre part, E { + E 2 s’annule avec E 4 et E 2 . 
Done les termes du second degre sont nuls, done 

H# = 0. 


<17) 

j dll dv 1 . D rfv* 

1 dE[ ~ 1 rfEi ' 2 rfE, 

I dH dti d ~> j 

( dE\ ~ 1 dEl ^ 2 dEt 

d’ou 


{18) 


Soit 



Done les coefficients de E* et de E* sont nuls quels que soient 
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a et E 3 dans le developpement de H; ils sont done nuls, si a = o 
et quel que soit E 3 dans le developpement du terme constant de - • 


371. Soit main tenant 


Soicnt 


H 4 = aE\ - 4 - 2 &Ef E| - 4 - cEf. 

v i — Xt •+■ (AjEf h- pi Ef +■..., 

V 2 = Xo * 4 “ [i.o Ef -+* JJ -2 ^2 "+"*••? 

B t =p E? 4 -..., B 2 =* Y E|-h... 


les premiers termes des d^veloppements de v,, v 2 , [3, y suivant les 
puissances de E, et de E 2 ; ces coefficients a, 6 , c, \ p, (3, y sont 
eux-memes des fo notions de E ;l ou de a, on de E 3 seulement si nous 
supposons a = 0 . 

Les equations ( 17 ) nous donnent alors 


4 & Ef ~ 1 4 b Ej E| -4- , . . = 2 E i -4- 2 Y(J-2 Ei K| -4- . . . , 
4 6 E? E 2 -h 4 cE| -4- . . . = 2 pfXi EJ E 2 -+- 2 Y Pi E| . ., 

d’ou 

2 Cl=P[A h 2 d = 3 ll\ =yn 2l 2 C = V [2 2 , 

ou 

a __ b b _ c 

r ^2 ~~ [4 * 


G’est la une relation entre les coefficients du developpement 
de et de v 2 , et par consequent de c et de g d’une part, et ceux du 
developpement de H, et par consequent (en supposant ol — o) du 

terme constant de -• 
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ACTION DES PLANETES. 


372. Pour etudier Paction d’une planete troublante sur le sys- 
t&me forme par le Soleil et une planete troublee, on commence 
par former les equations du mouvement de ce systeme comine si 
la planete troublante n’existait pas. Ge mouvement est alors keple- 
rien et, en seconde approximation, on eLudie les perturbations 
de ce mouvement keplerien par la planete troublante; pour cela 
on applique la methode de la variation des constantes. 

Nous opererons absolument de la meme maniere pour etudier le 
mouvement du systeme quadruple forme par le Soleil, la Terre, 
la Lune et une planete. Comme premiere approximation, nous 
integrerons les equations du mouvement du systeme triple : Soleil, 
Terre, Lune; c’est ce que nous avons fait dans les Chapitres pre- 
cedents; nous avons obtenu amsi les coordonnees des trois astres 
de ce systeme en fonction du temps et d’un certain nombre de 
constantes d’integration G. Nous devons ensuite Etudier les pertur- 
bations de ce mouvement par la planete, e’est-a-dire determiner les 
petites variations des constantes G dues a Taction de cette planete. 
Cette fa^on d’appliquer la methode de la variation des constantes 
a £te proposee et mise en oeuvre par M. Newcomb. 

Reprenons les notations des n os 42 (t. I, Chap. 11} et 312 
(Chap. XXIV). Soient 

A la Lune, 

B le Soleil, 

C la Terre, 

P la planete, 

D le centre de gravity du systeme Terre, Lune, 

G celui du systeme Terre, Lune, Soleil. 

P. — II (a). 


8 
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Soient 







#i. 

#2, 

#3, 

mi ~ 

772 2 

= 7 ? 2 3 les 

coordonnees et la masse de 

A; 

#;, 

# 6 , 

# 6 , 

m w — 

7 ? 2 S 

= ^6 

» 

B; 

#77 

# 8 , 

#9, 

mi = 

772 * 

== 772 9 

» 

C; 

# 10 , 

# 11 , 

# 12 , 

/?l 10 = 

THU 

= ™12 

» 

P; 


Soienl 


7i = rrii ■ 


#17 

#21 

#3 

# 4 , 

#i, 

#; 

#107 

#117 

#i 2 


= t? 2 ; 

ass 771; 


= m' 5 

= m\ 

™io 

= ™ii 

= 772, 


t t / Tn \ ^7 

m* = 771 f, 838 m 3 = : 

1 2 3 m! -f- / 7 Z 7 

, _ Jn 4 (m,H- 7 n 7 ) . 

77 ? r — Til t; " c “ — - » 

* 5 6 /TZi-r- t? 2 4 4- m 7 

/ 7 nio(/?ii-H /n 4 4 - m 7 ) t 

10 11 mi 4 - 4 - 772.7 4 - m 10 

, ,dx'i 

Soient 

t __ i yz! = iy/! 

2 ^ m 2 rri 

Fenergie cinetique et 

r T _ ( m,i ( m l m 1 77247227^ ( m i0 t ' 

U — \ ~AB 1 AG 1 BC - / — \ ~ PA 1 PB~ 


Fdnergie potentielle; F = T -J- U Fenergie totale; nous aurons les> 
equations canoniques 

dx’j _ <fF __ a?F 

~ 7 * 

Nous diviserons T et U en plusieurs parties; nous poserons 
T = Tj 4 - Tj4- T 3 , 

U = U j 4- Ui 4- U3 4- U 4 4- U 5 4- U 5 , 
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et 



(*' = !> 2 , 3), 

(* = 4, 5, 6), 

(i = 10, II, 12), 


rT «i m- , mi) 

= AC~ ’ Us = BD ’ 

j t m i 0 (m t -+- /?Zt -f- m-) 

u * = pg ’ 

U$ = ui\ mm ^ pQ pg ^ ™io( H— nii ) ^ ~~ ■ pjQ ^ > 

U 6 = 7?? i m io ^ pg J ■+" m ~ m io ^ pg pc) * 

On voit que T\ et U* dependent seulement des coordonnees de 
la Lane (et des y\ correspondants); T 2 et U 2 des coordonnees du 
Soleil; T 3 et U 4 de cedes de la planete; U 3 de cedes de la Lune et 
du Soleil; U 5 de cedes de la planete et du Soleil, et enfin U 0 de 
cedes de la Lune, du Soleil et de la planete. 

373. 11 faut maintenant voir quel est Fordre de grandeur de ces 
differentes quantitds. Je supposerai que Fon ait pris une unitd de 
longueur de Fordre de BD, de telle sorte que AC soit de Fordre 
de la parallaxe a, ce que j’ecrirai 

BD ^ i } AG a ; 

je prendrai de m£me une unit^ de masse de Fordre de la masse 
du Soleil m 4 , de sorte que 

i * 

Pour les planetes inferieures, on aura 

m io /n 7l PD rsj i. 

Pour les grosses plan&tes, m i& sera beaucoup plus grand que m 7 ; 
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mais en revanche PA, PB, PC seront beaucoup plus grands que i , 
et cela fera une sorte de compensation; nous admettrons done 
dans tous les cas 

7?l 10 ~ ^ 7 , PD^I. 

Nous trouvons ainsi 

rr tt tn x m n 

U 1 ~ > 

a 

To ~ L *2 ~ ^7) 

U3 ~ a 2 mj, 

T3 ~ U4 ~ m-, 

U 5 ^ 77if, u 6 /^ 

Nous poserons maintenant 

F = F'-+-F% 

FJ = T2+T3+U2+U4, F' = FJ •+• U„ 

f = t 2 4-t 3 +u 2 +u,+ u 5 , 

F ff 0 « T, -h U, -+• U 3 , F ff = F" + U 6 . 


Nous observons : 


i° Que F 7 ne depend pas des coordonnees de laLune; 

2 0 Que U3 et U 0 sont n^gligeables devant F 7 ; 

3 ° Que T< et U 4 ne dependent que des coordonnees de laLune. 


II en r^sulte qu’en ce qui concerne les coordonnees du Soleil et 
de la planete, nous pouvons nous contenter des equations 


daft dF dri 

dt ~~ dy r f dt 


dF' 

dx'i 


( i = 4, 5, 6, 10, ii, 12). 


Pour la determination des coordonnees de la Lune, nos equa- 
tions se reduisent a 


(2) 


dx^ 

dt 


_ dF " 


dfi 

dt 


dF" 

dx'i 


(t = i, 2, 3). 


374 *. En premiere approximation, nous negligerons U 5 devant 
F' 0 , et U c devant F'' . Dans ces conditions, nos equations se reduisent 
a 


(1 bis) 


dx'i 

d F£ 

dy'i _ 

dFjl 

dt 

" dy'i ’ 

dt 

dx’i 

dx'i 

d F'i 

dy'i _ 

_^n 

dt 

“ dy'i ’ 

dt 

dx'i 


(2 bis) 


(t= 1,2,3), 
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Or on voit que se compose de deux parties, Pune dependant 
seulement des coordonnees du Soleil et Pautre de celles de la 
planete, et queF" n’est autre chose que la fonction m\ dun°312 
(Chap. XXIV). D’ou cette consequence que, si l’on se borne aux 
equations (1 bis) et (2 bis), le mouvement du Soleil B par rapport 
au point D et celui de la planete P par rapport au point G sont des 
mouvements kepleriens. 

D’autre part, le mouvement de la Lune par rapport a la Terre 
est celui quia ete etudie dans les Chapitres XXV a XXIX. Nous 
supposerons done qu’on a completement determine ce mouvement 
en appliquant les procedes exposes dans ces Chapitres. 

Les quantites 

dh y'i (8 = 4,5,6,10,11,12) 

s’exprimeront done en fonctions du temps et des douze elements 
(canoniques ou elliptiques, cf. n° 08 ) de Porbite de B autour 
de D et de celle de P autour de G, ou bien, si Ton prefere, en 
fonctions des deux longitudes moyennes de B dans son mouve- 
ment kepl^rien autour de D et de P dans son mouvement kepl^- 
rien autour de G, et des dix autres elements (canoniques ou ellip- 
tiques) des deux orbites. 

Les quantites 

x'h y'i (* = i,a,3) 

pourront s’exprimer en fonctions du temps, des six Elements de 
Porbite elliptique de B autour de D et de six autres constantes 
d’integration, ou bien encore en fonctions : x° de la longitude 
moyenne du Soleil, e’est-a-dire de l’argument ; 2 0 des cinq autres 
elements de Porbite elliptique du Soleil; 3° des trois argu- 
ments t, t*, x 2 ; 4° des trois constantes E*, E 2 , m (ou de trois 
fonctions quelconques de ces trois constantes etdes cinq elements 
de Porbite solaire). 

373. Ainsi nos 18 variables x- Ll y\ se trouventexprim^es enfonc- 
tions de 0 arguments variant proportionnellement au temps, qui sont 
les deux longitudes moyennes de B et de P, et les trois arguments x 7 
x i7 et de i3 constantes d’ integration. 

Quand on a int£gr£ les equations (1 bis ) et (3 bis), on connait 
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les relations qai relient les 18 variables, ces 5 arguments et 
•ces 1 3 constantes. 

Supposons maintenant que nous poussions plus loin l’approxi- 
mation et que nous revenions aux equations (i) et (2). Nous 
pourrons alors definir 18 variables nouvelles qui seraient lides 
aux 18 variables anciennes x 1 et y 3 par les mdmes relations que 
Fdtaient nos 5 arguments et nos i3 constantes quand nous nous 
•contentions des Equations (1 bis) et (2 bis). Ces 18 variables 
pourront dtre regardees comme les elements oscalateurs des 
troisorbites deB autour de D, de Pautour deG, de Aautour de G. 
Seulement ces Elements osculateurs ne seront plus, les uns des 
fonclions lindairesdu temps, les autres des constantes; tout ceque 
nous pouvons dire, c’est qu’a cause de la petitesse des termes 
compldmentaires U 5 et U c , les uns varieront a pea pres propor- 
tionnellement au temps, et les autres tres lentement . 

Operant tout a fait comme au n°79(t. I, Chap. IY), nous allons 
faire un changement de variables en prenant pour variables nou- 
velles ces 18 elements osculateurs; mais, pour l’application de la 
mdthode de Lagrange, il convient de choisir ces variables (que 
nous n’avons pas encore completement definies) de telle fagon que 
la forme canonique des Equations ne soit pas altdree. 


i° Pour les 6 elements du Soleil, nous choisirons les elements 
canoniques 

l- 1 ! X, "Hsi 


ddfrnis au n° 88. La longitude mojenne X n’est autre chose que 
noire argument t 3 ; d’ailleurs y/S? — h* sera de l’ordre de E 3 , et 
pourrait jouer le merne role que E 3 dans I’analjse des Chapitres 
prdcddents. Dans ces conditions, 

(3) x'idyl-h dy\ -+- dy' 6 — \dL — d\ 

est une diffdrentielle exacte, ce qui est la condition pour que les 
equations conservent la forme canonique. 

2 0 Pour les 6 Elements de la planete, nous choisirons dgalement 
les dldments canoniques 

du n° 88. Mais d’ailleurs ces dldments ne joueront aucun rdledans 
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P analyse qui ya suivre; et en ce qui les concerne, les Equations (i) 
nous donneraient simplement les perturbations du mouvement de 
la planele par le systeme Terre-Lune suppose r^duit a un point 
mathematique; ces perturbations out ete determines dans le 
Tome I. 

3° Pour les 6 Elements de la Lune, nous prendrons les 3 argu- 
ments t, z { , et les 3 constantes B, B*, B 2 du Chapitre precedent. 
Si nous rapprochons les fonnules 


(Xy — Y#) du, 


rii 


2*'4r'=2^ x 

'^xdX-d&^^X'dw 

S', = <£,+ n 3 {Xy — Yx), 
u = wt = x 3 , y'i=m\f h 

dX — rffl' =2 B dz— E 3 $> dx„ 

2*' d/~ dQ -= 2 


des n os 318, 320 (Chap. XXIV) et 360 (Chap. XXIX), nous 
pourrons ecrire 

m\ dO ? ^ m\ B* m\ E 3 <J> d-z z — m r t (Xy — Y x) 


Toutes ces formules supposent que les Elements du Soleil, et en 
particular E 3 , sont regardes comme des constantes* Si nous sup- 
posons de plus di z = o, nous yerrons que 1’ expression 

( 4 ) x \ dy\ *+- x r 2 dy\ -4- x\ dy\ — m\ ( B d x -+■ B , h- B 2 ) 

est une different! ell e exacte, en supposant que les six dldments 
solaires (en y comprenant E 3 et t 3 ) soient regardes comme des 
constantes. 

II vaudra mieux d’ailleurs Ecrire la relation pr^cddente sous la 
forme suivante (en divisant par zra'), 

do." =2®' dy—^i B dx+- E»4> dz 3 — (X^— Y*) dz 3 , 


ou, mieux encore, nous remarquerons que t 3 ne joue pas le m£me 
rdle que les autres arguments t et nous mettrons en Evidence le 



CHAPITKE XXX. 


1'iO 

terme en en posant 

B clz = B cl'z -+- B j cl'z i h— Bo dz 2 

ct en Ecrivant par consequent dans les formules precedentes 

an lieu de^Bcfc, ce qui donne 

dQ” *' dy" B dr. — d-. s (B 3 - E s 4> - Xy -4- Yz ). 

Nous rappellerons ensuite la formule 

<j> r 

Zl = Kh-E 3 ^ 
n 2 

du n° 360, et nous poserons 



Nous aurons alors 

1 5) 

cette formule suppose que tous les elements solaires, sauf t 3 , sont 
des constantes. Si nous regardons de plus t 3 comme une constante, 
nous aurons P expression 

(4 bis) 

qui sera une differentielle exacte. 

376. Nous n’avons pas a revenir sur ce qui concerne les Equa- 
tions (i). On integre d’abord les equations approchEes (i bis), qui 
dEfinissent le mouvement kEplErien de la planete par rapport au 
point G, et du SoJeil par rapport au point D; on en dEduira l’in- 
tEgrale des Equations exactes (i) par la mEthode de la variation des 
constantes. Ce n’est pas autre chose que l’Etude des perturbations 
mutuelles de la planete et du systeme Terre-Lune (supposE con- 
centrE en son centre de gravite D), Etude qui a EtE faite complete- 
znent dans le Tome 1. Passons done aux Equations ( 2 ), que nous 
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ecrirons 

, F* d z 

dx\ _ dy'[ _ m' t . 

{b) dt~ dfi ’ dt dx’t 

Nous allons prendre pour variables nouvelles les ^ldments 
canoniques solaires et les six variables t, t 2 ? B, Bj, B 2 . Les 
coordonnees x f , les etQ v vont^tre desfonclions de cessix varia- 
bles t et B, de etdes cinq autres elements solaires que j’appel- 
lerai les y/. 

La formule (5) suppose que ces y; sont constants; si on les 
regarde comme variables, cette formule doit 4tre completee et 1 on 
doit ^crire 

( 5 his ) dQ!' =^P i x ' dv H — ^ B dt — ( G — 4>i ) -4-^ T id^i. 


en posant 


r /= 


dQ" , dr" 

d'(t Zd d'(i 


Pour savoir ce que vont devenir les Equations (6) par ce chan- 
gement de variables, il fain se reporter au theoreme du n" 12. 
Dans ce numero on envisage un systeme d’dquations canoniques 
ou F depend explicitement du temps; on fait un change ment de 
variables, les variables nouvelles x 1 et y f etant fouctions des va- 
riables anciennes x et y et du temps t. On suppose que 




soit une diflKrentielle exacte, en snpposant dt = o, et que 


cc' dy ' — df — W dt 

soit une diflerentielle exacte pour dt ^o ; dans ce cas les equations 
conservent la forme canonique, mais^a fonction F doit £tre rem- 
plac^e par F — W. 

Nous pouvons appliquer ici ce theoreme, car, les Equations (i) 
£tant int^grees, les elements y* etT 3 peuvenl^tre regard^s comme 
des fonctions connues du temps. Nous pouvons done dcrire 

dQ a x'dy" — ^ Brft + W dt, 


(5 ter) 



11 % 
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en posant 



W=(Oi 

— G)~ 

d~ z ' V 

dt 2d 


et les equations (6) deviendront 



(7) 

dB 

dt> 2 

dx _ 

cft>, 

dt ~~ 

dt ’ 

dt 

rfB ’ 

en posant 







m. 

- W. 



II est bon de remarquer que cette analyse ne suppose nullement 
qu’on ait choisi pour les y/ les elements canoniques du Soleil, 
qu’on peut tout aussi bien choisir des fonctions quelconques de 
ces six elements canoniques et en pardculier les elements ellip- 
tiques. 


377. Nous devons parmi les elements solaires faire une distinc- 
tion; les coordonn^es x , y , telles qu’elles ont ete calculees 
dans les Chapitres precedents, dependent deT 3 , anomalie moyenne 
du Soleil, ainsi que de l’excentricite solaire E 3 , et du grand axe 
de Forbite solaire inversement proportionnel a la parallaxe a. 
Elies ne dependent pas des trois autres elements qui sont la longi- 
tude du perihelie terrestre et les angles qui definissent Forienta- 
tion de l’ediptique mobile par rapport aux axes fixes. II en est de 
m£me de X, Y, Z et Q ;/ ; ces fonctions sont egalement indepen- 
dantes de ces trois derniers elements, que j’appellerai les elements 
d } orientation. 

Nous avons dcrit plus haut 

2 ^' dy = ^ » dX -4- c Xy - Y») d*u 

raais en supposant d^t = 0 . En regardant les y* comme des variables 
il convient d’ecrire 


d’ou 


2*' dy = 2* dx H- c Xy - Yx) dx t A, d v , 
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Si y; est un element d’orientation, on aura 


et par consequent 


dor dx 

d*{£ ~ d'{i 


r* = 


d&” ^ dX 

_ — 2d X d'a 


— A./ = — A*. 


Qu’est-ce maintenant que A/? Si les trois Elements d’orienta- 
tion y*, y 2 , y 3 subissent trois accroissements <fy l9 <fy 2 , cfy 3 , tout 
se passera comme si les axes mobiles des .r, y, z subissaient trois 
rotations infiniment petiles, par rapport aux axes fixes des x \ , x'. 2 , 
Si Ton s’arrange pour que les axes mobiles coincident a l’ori- 
gine des temps avec les axes fixes, ces trois rotations s’op^reront 
autour des trois axes et Ton aura 


^37 = TjX XZj 


d; 


Tar. 


378. Si nous supposons que la masse m 40 soit nulle, les y; se 
r^duisent a des constantes ainsi que et l’on a 

4' = “' 

On a d’autre part 

F'-FS = T 1 +U 

d’ofi. finalement 

<J>* = G. 

G est fonction simplement des B et des y*; les Equations canoniques 
se r^duisent a 


dB dG dxi dG 

~dt ~ ~dk. = °’ "5? =V/ = — 5iT<* 

Nous pourrons done ecrire 

— dG — v -h Vj dB\ -4- Vj 



CHAPITRE XXX. 


en regardant les y comme des constantes. Si nous rapprochons de 
la formule du n° 363 


nous en tirerons 


dB “ JB dbi *+■ dv\ - 4 - B 2 dv%, 


H — G = B v -4— B 1 v 1 B 2 v 2 , 


plus une fonction arbitraire des y/ que nous pouvons supposer 
nulle. Nous aurions pu deduire celte m&me formule de la for- 
mule ( 9 ) du n° 363 qui peut s’ecrire 

H = Bv + B I, V j — h“ B 2 Vo B3V3-4” ICV3 
et de la definition de G au n° 375 qui peut s’eerire 

G — B 3 v 3 -f- Kv 3 . 

379. Supposons maintenant que m l0 ne soil pas nul; alors ~ 

ne sera plus nul, mais tres petit; — ^ ne sera plus egal a une con- 
stante, mais nous pourrons poser 

r 

— = 1 -4- s, 

/i 2 dt 

e etant tres petit. Nous aurons d’autre part 


d’ou enfin 




Nous pouvons poser encore 

G = G 0 ■+■ 3G, 

G 0 ^tant ce que devient G quand on y remplace les y* par leurs 
valeurs initiales et etant, en consequence, fonction seulement de B, 
et B 2 , et SG etant tres petit, puis ecrire 


R =~7+( G “*i) S 




<J>2 = G - 4 - R, 

ce qui nous donne finalement les equations 


dB_dK 
dt “ dz 9 


dG_dR 
dB dB 9 
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R etant tres petit; nous pouvons dans les derivees de R remplacer 
les variables par leurs valeurs approchdes ; R et ses derivees pour- 
ront alors £tre regardees comme des fonctions connues du temps. 
Ce seront des fonctions periodiques par rapport aux cinq argu- 
ments 

T, Ti, T 2 , Tj, 

etant l’anomalie moyenne de la planete, lesquels arguments, r&- 
duits a leurs valeurs approchees, sont des fonctions lineaires du 
temps. Les derivees de R se pr^senteront done sous la forme 

2a.cos(cc<-i- ft), 

les A, les a et les (3 etant des constantes. 

On aura alors les B par de simples quadratures; les B etant de- 
termines et les yPayant6t£ prealablement a l’aide des Equations (i), 
il en sera de meme des derivees ^ qui ne dependent que des B et 
des y ; on pourra avoir les t par quadrature. C’est absolument ce 
que nous avons fait dans les Chapitres IV et V. 

380. 11 y a une autre maniere de comprendre l’application de la 
methode de la variation des constantes. Si m ]Q 6tait nul, nous 
aurions certaines relations entre nos variables , y\ (i = i, 2, 3), 
les elements lunaires B et t et les elements solaires t 3 et y /; soient 

(9) y> =/(B, y/,tj) 

ces relations. Si m i0 n’est pas nul, de sorte que les elements de 
l’orbite solaire soient variables, nous pouvons conserver les rela- 
tions (9) comme definitions des elements osculateurs B et t; c’est 
ce que nous avons fait jusqu’ici, mais on peut aussi op^rer autre- 
ment. 

Soient y? les valeurs initiales des y/; soit 

= /l$t -h £ 0 , 

n[l et e 0 £tant les valeurs initiales de n*> et de t 3 ; remplagons alors 
les Equations (9) par les suivantes : 


(9 *w) 


Xi 

ft 


f ( B, *C, yj , ). 
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Ces six Equations (9 bis) definiront les six Aments oscula- 
teurs B, B 4 , B 2 , t, t 4 , tt 2 - Les B et les c definis par les equa- 
tions (9 bis) ne sont pas les m£mes que les B et les t definis par 
les Equations (9), mais ils en different fort peu. 

Remarquons que, les y®, /ig, e 0 etant des conslantes, les elements- 
solaires variables n’interviennent pas dans la nouvelle definition. 

Que devient V equation (5 bis) 


On doit y remplacer les y/, n 2 et r z par y®, /i® et ona alors 




dx% 

n i 


= City 


puisque les y® et e 0 sont des constantes, et il reste 
dor = X 1 dy - B dx - ( G„ - <*>' 5 ) dt. . 


De plus, comme <J> 4 depend des constantes solaires, il faut rem- 
placer <£ 4 (qui est une fonction des x , des y', de t 3 et des y*) par 
ce que devient cette fonction quand on y remplace y** et par y?* 
et t®; je designe par le r^sultat de cette substitution, d’ou 

W = <*>« — G 0 . 


Nous aurons alors les equations 


(7 bis) 


dB _ d<P 2 dx __ d< t> 2 

dt ~~ dx ’ dt~~ dB 


en posant 
d’ou 


$ s = G 0 -t- R, R= — ' +(<!>! — 4>J). 

TiX j 


Remarquons que G 0 ne depend que des B. Nous retombons done 
sur des equations de la forme (8) qui se traitentde lameme fagon. 

M. Newcomb ( Investigation of inequalities in the motion of 
the Moon produced by the action of the planets, Washington, 
Carnegie Institution, juin 1907) emploie un procede mixte; il 
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emploie, en effet, le procede da n° 379 pour les Elements d’orien- 
tation et celui du n° 380 pour les autres elements. 

381. En general, les termes provenant de Paction des planetes 
sont fort petits et ne deviennenl sensibles que par l’effet des petits 
diviseurs si la p^riode est tres longue. Ils ne seront le plus souvent 
appreciates que dans le cas d’une double integration, qui intro- 
duit au denominateur le carre du petit diviseur. Nous aurons done 
la partie la plus importante d’un terme a longue periode en nous 
bornant a Integration des equations suivantes : 

d . B dB dx dQt 

dt ~~ dx 9 dt~~ dB 


Nous n^gligeons ainsi dans le second membre de la seconde 6qua- 
tion le terme — qui ne subirait qu’une simple integration. Si 

nous appelons SB, Sy, St les inegalites dues a un terme donne SR 
de la fonction perturba trice, et si g est ce que devient G quand on 
y remplace les B et les y par leurs valeurs initiales, nous pourrons 
encore ecrire 


, N d SB d 6R 

< 10) -ar=nr’ 


d 5 xj y d z G 

dt ~~ dBi dBj k 




d* G 
dB £ d'(k 


$7 k- 


382. Nous devons distinguer V action directe et 1 'action indi- 
recte d’une planete. Si l’action directe existait seule, tout se passe- 
rait comme si, le Soleil B et la planete P assujettis a decrire des 
orbites kepleriennes, I’un autour de D, l’autre autour de G, la Terre 
et la Lune ^taient soumises a l’attraction de ces deux astres mo- 
biles; si Paction indirecte existait seule, tout se passerait comme 
si, la planete n r existant pas, le Soleil B ^tait assujetti k decrire, 
par rapport a D, l’orbite troublee due a Paction de la planete. 

Dans le cas des Equations du n° 380, Paction directe pro- 

vient du terme et Paction indirecte du terme — <£? Fdans 
ce cas, dans les Equations ( 10 ), on doit remplacer G par G 0 et les 


derivees sont nulles, car G 0 ne depend que des Bj . 

Dans le cas du n° 379, on obtiendra Paction directe a l’aide 


des equations (8), en regardant les y ,■ et^~ comme des constantes 
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et en reduisant Rau terme • On obtiendra Paction indirecte, en 

m l 

definissant Jes variations des y* et de par le moyen des equa- 
tions (i) et en reduisant R aux termes 

(G 

Lorsque la planete troublante est tres voisine du Soleil, Paction 
directe et Paction indirecte sont sensiblement egales etde signes 
contraires. 

En effet, soit G, le centre de gravite de P et B. Dans le cas qui 
nous occupe, le point G! decrira une orbite sensiblement keple- 
rienne autour de D, analogue a celle que decrirait le point B si la 
planete n’existait pas, mais dont le grand axe, pour un meme moyen 
mouvement, se trouve multiplie par 



Si a la limite nous negligeons la distance PB, tout se passe pour 
Paction directe comme si la masse du Soleil etait augmentde 
de m i0 et, en ce qui concerne Paction indirecte, comme si la dis- 
tance BD etait multipliee par 



il y a done compensation en ce qui concerne les inegalites corres- 
pondantes (e’est-a-dire celles qui ne dependent pas de 
P angle PG, D). 

II y a compensation aussi pour les in^galitds proportionnelles a 
la distance PB (et qui contiennent cet angle PG*D comme argu- 
ment), car Pecart BG< multiplie par m h est egal a PecartPG* mul- 
tiplie parm 10 . La compensation n’existe plus pour les inegalites 
proportionnelles aux puissances superieures de PB, mais ces ine- 
galites sont beaucoup plus faibles. 

383. Nous devons distinguer deux sortes d’inegalites planetaires 
periodiques. Celles de la premiere sorte sont celles dont l’argument 
ne depend que de t 8 et de t*. Dans ce cas on a, en se reportant 
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aui equations (io), 
d'oii 

SB = o, 

Ces inegalites sont clone dues presque exclusivement a Paction 
indirecte (je veux dire que les termes provenant de Paction indi- 
recte subissent seuls une double integration), et Paction directe est 
negligeable surtout si la periode est longue. 

Les seuls y* dont depend G sont E 3 et le grand axe solaire a ! . 
Le terme enE 3 est de beaucoup le plus petit. On pent clone eerire 

d fa/ __ __ d*G . , 
dt ~~ d&i da' 0<Z ’ 


d SR 
d'z 


d 

dt 


-2 


d 2 G ^ 
dliidy 


et comme on a d ’autre part, A etant la longitude solaire, 


d 8X ^ 

dt ~ 


^ , 

oa , 

•2 a 


on voit que les inegalites St, St*, 8t 2 sont sensiblement proportion- 
nelles entre elles et a Pinegalite solaire ok. L’inegalite St engendre 
une inegalite a longue periode de la longitude vraie cle la Lune; les 
inegalites St* et St 2 engendrent des inegalites concomitantes a 
courte periode de la longitude vraie. La principale deces inegalites 
de la premiere sorte est engendree par Venus et a pour pe- 
riode 8t 4 — i3t 3 . 


384. Les inegalites de la seconde sorte sont celles dont Pargu- 
ment depend de t, t 4 ou t 2 . Pour celles-la SB* n’est plus mil. An 
contraire, comme Sy* ne contient que des termes independants 
de t, t* ou To, on devra faire 3y* = o, d 7 ou 

dfai _ drG ^ 

~~~ 2d dB £ dB A k ' 

Ces inegalites sont dues surtout a Paction directe. Pour les de- 
terminer, nous ecrirons Uc en negfigeant la paraitaxe sous la 
forme 

U 6 i — 3 cos 2 PDA k n9 
m[m l0 ~ 2 PD* ’ 

P. - II (a). 
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OU 

* 

-Hil — — ■ 77~ ~ - ( XX* -H 3 -+• 1 2 3S- r -+- [2 (D&y-h 12 C& ), . 

m\ 4PI>° 


X = ■+* oc'i — 2^2 3 , "lib = — #S*» 

3 = #2#3, CD=^ / i^s, ^ ? 


et ou X', Hb' 7 3', O', £' sont formes avec les trois projections du 
vecteur PD, c’est-a-dire avec 


ou 


a?',, -hXiP's, a/ 12 4 -Xa?' 6 , 


X = 


JYl\ 


tilt* H- /?Zj -+- ^7 


comme db, i»t>, £, (6 avec a?' , a?' 2 , #' 3 . 

Onvoitque ^ 5 , ••• dependent settlement de tt s et t 4 , 

tandis que Jb, Ht>, ... ne contiennent que des arguments depen- 
dant des coordonndes lunaires, a savoir (si Ton se borne aux termes 
elliptiques) des arguments de la forme 

zj'zz-h/czi pour X, 

2t - 4 - 2'r 3 -h pour -ife et C, 

t -H r C3-H(v'““ I^a+^J P our ^ ®t (0. 


Les principales inegalites de cette sorte sont celle de Hansen, 
due a Paction directe de Venus, quia pour argument 


Zi -f- 1 6^:3 — i8t 4 (periode 23 q ans) 


provetfant du terme en?* de X et du terme en i 6 t 3 — i 8 t 4 de-^j 
et celle de Neison, due a Taction directe de Jupiter, qui a pour 
argument — "h — 3t 4 (pdriode 37 ans), provenant du 

terme en a'r-j-aTs — de tlb et et du terme en 3^4 de pgg 



Je me bornerai a re nv over pour plus de details au Memoire cite 
plus haut de M. Newcomb, ainsi qu’au Memoire de M. Radau 
dans Je Tome XXI des Me moires de V Observatoire et au rdsum<§ 
qu'cn a fait Tisserand dans le Tome III de sa Mecaniqae celeste. 
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38o. Pour ^tudier les accelerations sdculaires , nous devons 
d’abord nous reporter a ce que nous avons dit au Chapitre V, 
n° 105, au sujet de Finvariabilite des grands axes. 

Soil plus g^n^ralement un systeme d ’equations canoniques 


0 ) 


dxi dF dyi dF 

dt dyt* dt dxi* 


F = F 0 -{- pR, 


p. etant tres petit. Je distinguerai deux sortes de variables a?/, que 
j’appellerai les x\ et les x je designerai de raeme parj/'- ety". les 
variables conjugu6es des x\ et des x \ . 

Je supposerai que F 0 depend seulement des x f et est indepen- 
darit des x ", des y f et desj)/'; quant a R, il depend des quatre 
sortes de variables, mais il est periodique par rapport aux y ! et 
auxj/'. Nous supposerons en outre que R depend directement du 
temps; plus pr^cisdment, nous supposerons que R est periodique 
par rapport aux y r : aux y 31 et a un certain nombre d’arguments w 
qui sont des fonctions lineaires connues du temps. 

En premiere approximation, od a 


dx' _ dF 0 _ dx n _ ^F 0 _ 

dt ~ ~dy> ~ °’ dt ~ dy ~ °’ 

x' = const., x" = const., y" = const., 


dt “ dor* 1 



dF 0 
dx 1 


= const. 


Pour la seconde approximation, nous remplacerons dans les 
derivees de R les variables par les valeurs approch^es que nous 
venons de trouver et nous aurons 
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et le second membre sepresentera sous la forme d’une serie trigo- 
nometrique. 

En eflfet, R est une fonction periodique des y' : des j/ et des w ; 
les w sont des fonctions lineaires connues du temps, il en est de 
m^me des premieres valeurs approchees des y’\ quant a celles 
desy ; , ce sont des constantes. 

Si nous supposons qu’il n’y a entre les et les (c’est- 
a-dire entre les mojens mouvements) aucune relation lineaire a 
coefficients entiers, les variations s^culaires ne pourraient pro- 
venir que de ceux des termes de cette serie trigonometrique qui 
sont independants a la fois des y / et des w. Or tous ces termes 
disparaissent quand on differentie R par rapport a y\ . Done 
il n’y a pas de termes sdculaires dans les x \ . 

C’est la la generalisation du theoreme sur l’invariabilitd des 
grands axes. 

Appliquons-le aux equations (8) du Chapitre precedent* G va 
jouer le rdle de F 0 , R celui de uR, les B celui des x\ les t celui 
des y\ et enfin ^3 et celui des w; nousn’aurons pas de variables 
analogues aux x ,! et aux j", tandis que dans les equations (5 bis) 
du n° 93, Chapitre IV, les L, les A, les p et les <0 jouaient res- 
pectivement le role des x\ desy, des x u et desy'; en premiere ap- 
proximation on a 


B = const., 


dz d& 

dt ~~ 


as const. 


En deuxieme approximation on a 


d% dK* 
dt ~~ dz* 


Dans le second membre on remplace les B, les y par leurs valeurs 
■approchees qui sont des constantes, les t, t 3 el par leurs valeurs 
approchees qui sont des fonctions lineaires du temps. 

On obtient ainsi une serie trigonometrique. Les termes de cette 
serie qui sont ind^pendants a la fois des c, de et de T4, et d 7 ou il 
pourrail r^sulter des variations seculaires, disparaitront quand on 
diff^rentiera par rapport a a ou a done les quantiles B, 
B 2 ri? eprouvent aucune variation seculaire. 

386. Le numero pr^c^dent nous apprend que les variations 
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seculaires des B, 

8B, 5Bj. 

sont nulles. (Test sur cette circonstance que s'appuie Brown pour 
determiner les accelerations seculaires 


OV, Svj, Sv 2 

des divers mojens mouvements. Pour cela rappelons la formule du 
n° 363 : 

dH = Bflfv + B[rfv|+n s dv z . 

Si nous regardons les B, H, y, et v 2 comme des fonctions de 



de E, et de E 2 , il viendra 


d H 

= B +B, 

dv\ 


Bj 

d'Jv 

i d'i 

d'j 


~dZ’ 

] dH 

= B, 

dv 1 


b 2 

d'>i 

j dEi 

dEi 


d E,’ 

I dH 

= B, 

dv\ 


b 2 

dv* 

[ dE« 

dE% 

-H 

dE<% 


Si nous negligeons Ej et E” et par consequent B* et B 2 qui sont 
respectivement di visibles par EJ et Ej, il restera 


et par consequent 


dH 

dv 


= B, 


,dH 

° &> 


= o . 


Les B, H et les v sont fonctions non seulement des trois constanles 
lunaires v, Ej et E 2 , mais encore de deux des constantes solaires, 
k savoir le demi-grand axe a' et l’excentricite E 3 . En vertu du 
tb^oreme d’ Adams, H ne contient pas de termes en E^ et Eij, de 
sorte qu’en n^gligeant les puissances superieures de ces quantiles 
on aura 

dH _ dH d* H _ d* H _ 
dEi dE% dv dEi dv dE% 1 

et qu’il reste 


s dH. 

' dv 


_ dm, 

” dv ° V 


dm 

d't d 


SE 3 = o. 


(3) 
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Comme 8E 3 est connu par la theorie des planetes, cette equation 
donnera 8v. Nous trouverons ensuite 


(4) 


8v/= 

av 


cbn 

dEl 


8E 3 . 


Nous devons remarquer, en effet, qu’a ce degrd d’approximation 
on a 


Ei 8Ej = Ei% 8 E 2 — o, 


et, comme B, etanl divisible parEf nous pouvons ^crire B 4 = C,EJ ? 
il vient 

(5) 8B 1 = Ef 0G, G! 3(Ef) = o, 

d’ou, puisque nous negligeons Ef, 

8(E?) = o, 

et de meme 8(Ei;) — o. 

Les equations (3) et (4) determinent les accelerations secu- 
laires ov, Sv { , ov 2 . Mais pour cela il faut se servir des expressions 
de v* et de en fonctions de v et de Ef, ou, ce qui revient au 
m&me, de cedes de c et de g en fonctions de m et de Ef. Il faut, 
d’autre part, connaitre l’expression de H en fonction de v et de E 3 . 

11 nous suffit de rappeler qu’en vertu du th^or&me d’Adams, 
quand on neglige la parallaxe, H n’est autre chose, a un facteur 
constant pres, que le terme constant du developpement trigono- 
m^trique de 


387. Avant de pousser plus loin rapproximation en tenant 
compte de Ef et Ef, commen$ons par remarquer que nous avons 


<6) 


dvt d- v/ 
dv " d>* 


d*'jj 

du <a?E s 


SE 3 , 


d’ou il r^sulte que les Equations (3), (4) et (6) nous donnent en 
premiere approximation 


f dvi 
dv 


Cela va nous permettre de passer a la deuxi&ne approximation. 
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Reprenons les notations du n° 371 et ^crivons 

H = H,+ a Ef-t- 2^EfE|-hcE| = H 0 -+-H 4 , 
vi = Xi ■+• Ef -+• y.'i E|, 

Vg= -+- (j-aEj-h p-2 E|, 

Bi*pE| f B 2 = y e 1* 

Dans la premiere approximation, nous avons r£duit H, v 4 et v 2 a 
leurs premiers termes, H 0 , et X 2 . 

Passons ala deuxi&me approximation. Comrae S(Ef) et 8 (E£) 
sont de l’ordre de Ef et E*, nous voyons que 


qui sont des polynomes du deuxieme degT<* en Ef, E*, S(Ef), 
8 (EJ), sont de l’ordre de EJ ou E* et par suite n^gligeables, de 
sorte qu’on a 


sj dH _ „dHo _ & H 0 
dv ° dv dv 2 


3v -H 


d* Hp 
dv dE 3 


aE s . 


D’autre part, 


dH 

dv 


B -+- Bi 


dvj 

dv 


B 2 


dv 2 
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d’ou, en nous rappelant que les SB sont mils, 


( 7 ) 


d 2 H 


d* H 


dv i 


. dv s 


-^r $v ' + "^7^ oE a = = 


dv dE 3 


dv 


dv 


Gomme dans tous les termes du second membre figure en. 
facteur Ef ou Eij, nous pourrons dans ce second membre rem- 

placer 8 ^— • et S ^ par leurs premieres valeurs approves, de sorte 
que l’equation ( 7 ) nous donnera la nouvelle valeur de 8 v. 


388. Nous avons ensuite, pour les Svi, 

(8) 5v/= ^ Sv -t- SEj + & 6(Ef ) ■+• « ( BJ ). 

II faudra cette fois, dans le caleul de tenir compte des 

termes p./Ef 4 - |xj EJ . 
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ACCELERATIONS SECULAIRES. 


Pour calculer o (EJ ), nous nous servirons de liquation (5) et 
nous y ferons 


Gi— p. 


5G t = as = 


Nous pourrons d’ailleurs, dans cette formule, remplacer 8v par 
sa premiere valeur approchee. On calcuJerait S(E^) de la mdme 
maniere. 

On peut considerer les developpemenls des v/ comme connus, 
et par consequent aussi ceux des A/, jjl/. jjl't . On connait egalement 

celui de et par consequent celui de H et ceux de H 0 , a, b , c. 

Quant a p et a y, iis nous sont donnes par le theoreme d’Adams 
du n° 371 , d ? oii 

IM A ’ ( A ’ 

On remarquera que, dans toute cette analyse, les elements 
d'orientation ne jouent aucun r<Me, d’ou il suit immediatement que 
les deplacements seculaires de l’ecliptique ne peuvent exercer 
aucune influence sur l’acceleration seculaire du mouvement de la 
Lune, ce qui confirme les resultats obtenus autrefois par M. Pui- 
seux. 

II resterait a parler de ce qui concerne les inegalites dues a 
Paplatissement terrestre ; en Fabsence de nouveaux travaux sur ce 
sujet, je me bornerai a renvoyer le lecteur au Tome III de la Meca~ 
nique celeste de Tisserand, pages t 44 et suivantes. 


FIN DE LA DEUXIEMH PARTIE DU TOME II. 
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Donner des ph&nomenes yiectriques une explication m6canique complete, 
r6duisant les lois de la Physique aux principes fondamentaux de la Dyna- 
mique, c’est lit un problkme qui a tente bien des chercheurs. N*est-ce pas 
cependant une question un peu oiseuse et oh nos forces se consnmeraient 
en pure perte? 

Si elle ne comportait qu’une seule solution, la possession de cette solu- 
tion unique, qui serait la v6rit6, ne saurait 6tre pay£e trop cher. Mais il 
n'en est pas ainsi : on arriverait sans doute it mventer un mScanisme 
donnant une imitation plus ou moins parfaite des ph^nombnes &ectrosta- 
tiques et yiectrodynamiques. Mais, si Ton peut en lmaginer un, on pourra 
en imaginer une infinite d’autres. 

11 ne semble pas d’ailleurs qu’aucun d’entre eux s'impose jusqu’ici it 
notre choix par sa simplicity. Dks lors, on ne voit pas bien pourquoi Tun 
d’eux nous Ferait, mieux que les autres, p£n6trer le secret de la nature. li 
en r^sulte que tous eeux que Ton peut proposer ont je ne sais quel 
caraet&re artificiel qui rypugne k la raison..* 

S’il est oiseux de ehercner it se repr£senter dans tous ses details le 
m£canisme des phynomhnes electriques, il est tres important, au contraire, 
de montrer que ces ph6nomhnes obyissent aux lois g^ndrales de la Meca- 
nique. 

Ces lois, en effet, sont independantes du mycanisme particulier auquel 
©lies s’appliquent. Elies doivent se retrouver invariables k travers la 
diversity des apparences. Si les phynomenes electriques y ychappaient, on 
oevralt renoncer a tout espoir crexplication m6canique. S’ils y obdissent. 




la possibility de cette explication est certain©, et Ton n*est arryry que par 
la difficult© de choisir entre toutes les solutions qu© le problem© comport©. 

Mais comment nous assurons-nous de la conformity des lois de rElee- 
trostatique ©t de I’EIectrodynamique avec les principes d© la Dvnamique? 
C’est par une s£rie de comparisons ; quand nous voudrons analyser un 

E henomene electrique, nous prendrons un ou deux phdnomenes mdcaniques 
ien connus et nous chercherons h mettro en Evidence ieur parfait paraliy- 
lisme. Ce parail41isme nous sera ainsi un garant suffisant de la possibility 
d’une explication mecanique. 

L’emploi de F Analyse mathymatique ne servirait qu’4 montrer que ces 
comparaisons ne sont pas seulement de grossiers rapprochements, mais 
qu’eltes se poursuiveni jusque dans les ddtails les plus prycis. Les limites 
ae cet Ouvrage ne me permettront pas d’aller aussi loin, et je devrai me 
borner a une comparaison pour ainsi dire qualitative. 
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